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FOREWORD 


comprehensive Symposium the conformity model prototype be- 
havior has been the minds the Executive Committee the Hydraulics 
Division for number years. The Committee Hydraulic Research (J. 
Stevens, Am. Soc. E., Chairman), appointed April, 1934, has an- 
nounced this subject one its major interests and the present Sympo- 
sium the result the efforts Subcommittee comprising Clarence 
Bardsley, Am. Soc. E., and Joseph Tiffany, Jr., Jun. Am. E., 
with Lane, Am. Soc. E., Chairman. Credit should also 
given Paul Thompson, Assoc. Am. Soc. E., and Herbert Vogel, 
Am. Soc. Colonel Thompson succeeded Major Vogel the com- 
mittee and was turn succeeded Mr. Tiffany. 

Two the papers were first presented the San Diego, Calif., meeting 
the Division July, 1941; one was presented the Denver, Colo., meeting 
the Division July, 1940; the remaining papers have not been presented 
publicly any form. Some papers this subject have been published 
separately and these are included the Bibliography (2) (3) Each author 
outstanding the art model experimentation and the analysis results, 
and has had opportunity determine how nearly the model tests conformed 
the prototype observations. 

The objective this research can stated single question: the 
prototype act the model rather difficult obtain complete 
conformity for two three reasons. For example, working conditions may 
not even approach the load limit that was used the model; and the prototype 
rarely built precisely accord with the model. The purpose the Sym- 
posium assemble the results wide variety studies which actual 
comparisons have been made that hydraulic engineers may judge the 
ciency the model tool design. 


Numerals parentheses, thus: (1), refer corresponding items the Bibliography, which appended 
the Symposium separate unit. 
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HYDRAULIC STRUCTURES 


JUN. AM. Soc. 


Some the general factors involved making model-prototype com- 
parisons hydraulic structures are mentioned this paper. Reasons for 
desiring comparisons, the types made, the difficulties encountered, and the 
measurements and instruments required are discussed briefly. Citing two 
the comparisons made the Bureau Reclamation examples, account 
given small irrigation structure, and the pressure and discharge tests 
the 102-in. river outlets dam western state. 


INTRODUCTION 


Although great progress was made during the decade the use 
hydraulic models aid the design hydraulic structures, there has been 
but little progress comparing the performance model and its prototype. 
considerable interest and importance make such comparison 
order determine the degree with which the model may relied upon for 
quantitative results even the laws similitude have been rigidly followed. 
the field hydraulic structures this particular value when applied 
discharge coefficients, water-surface profiles, erosion river beds, and pres- 
sures boundary surfaces, among others. nearly equal importance the 
fact that reliable prototype measurements will furnish accurate design data 
enabling more economical designs the future. 

The reasons why little progress has been made comparing the hydraulic 
performance model and its prototype may explained lack in- 
terest, particularly after the prototype satisfactory operation; natural 
reliance model theory; and the difficulties making prototype measure- 
ments. Fortunately, considerable interest now being shown for making 
these comparisons and believed that this Symposium not only will assist 
materially helping create more interest but, perhaps, from the exchange 
ideas, the difficulties involved making the comparison will lessened. 

Accordingly, this paper will treat some the recognized difficulties pro- 
curing prototype measurements, together with the measurements and instru- 
ments required. Two types comparisons between model and prototype will 
follow, one qualitative and the other quantitative. The former the easiest 
and the most frequently made. consists describing, for example, the ob- 
served performance some structure for several flows, and then comparing 
these observations with those made similar flows the model. recog- 


Hydr. Testing Laboratory, Bureau Reclamation, Denver, Colo. 
Lt., Army, Waterways Experiment Station, Vicksburg, Miss. 
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nized that this type comparison particular value only those most 
familiar with the structure, not too convincing others, and leads the en- 
gineer nearer actually verifying the quantitative results the model. 
Obviously, quantitative comparison the most valuable since not only 
will permit direct comparison between several quantitative factors the model 
and prototype, but will also begin establish the degrée with which the 
model can relied upon. The hydraulic model reliable high degree, 
because the laws similitude themselves are fundamentally sound when 
used correctly; but importance obtain concrete evidence and thereby 
recognize the influence any digressions made conducting the model ex- 
periments. 


Factors DIFFICULT 


Prototype Measurements.—The difficulty obtaining accurate prototype 
data the principal reason why correlation work has been delayed. Although 
several methods are available for measuring the flow, the difficulty lies not only 
actually being able procure the measurements, but also being able 
obtain data applicable the particular feature the prototype for which 
comparison made. these respects the discharge, velocity, and water- 
surface measurements are the most troublesome. 

Ordinarily the discharge passing structure can measured gaging 
station downstream. During flood flows, however, frequently impossible 
obtain accurate measurement because change river section due 
bed movement, because the peak flow may pass before complete measure- 
ment recorded. Occasionally flood will occur night, whereupon 
automatic gage, installed, must relied upon the measurement lost. 
the other sometimes necessary wait many years before 
structure operates, and even then sufficient range head and discharge may 
not obtained. Frequently the discharge passing structure may derived 
from the simultaneous operation spillway, outlet works, and power house. 
such event usually necessary wait until the flow derived from 
only one source, but that time the range measurements desired may have 
passed. Similarly, confluence streams exists between the structure and 
gaging station, heavy rainfall occurs, discharge measurement adjusted 
applicable the structure may questionable value. 

Perhaps the most trouble encountered when desired obtain ve- 
locity and water-surface measurements, not only for comparison with the model, 
but also for procuring fundamental data. the flow travels slow velocities 
and confined narrow channels, these measurements may made with 
minimum difficulty. the flow travels high velocities, however, the 
measurements are extremely troublesome, regardless the size channel. 
For example, how would one measure accurately the velocity and water sur- 
face the flow the spillway section dam western state the 
dam question here where the velocity excess 100 per sec, and the 
water surface affected air entrainment? impracticable use cur- 
rent meters pitot tubes measure the velocity, point gages obtain 
the water surface. 
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The futility using pitot tubes was realized the dam discussed here, 
more from the trouble maintaining the equipment than from the magnitude 
the quantities measured. During construction three streamlined, 
bronze piers, about in. high, were fastened bronze base plates the 
entrance and the lip the spillway bucket. Each bronze pier supports five 
stainless-steel pitot tubes whose static and dynamic openings are connected 
within the piers hydrostatic pressure cells similar those used for measuring 
gallery higher elevation copper tubing extending through the dam. 
provide check measurements, by-pass formed each cell and connected 
copper tubing manometers the gallery. spite all the precautions 
made, many the copper tubes leading the gallery became plugged and all 
attempts clear them failed. Some the diaphragms the pressure cells 
were ruptured during the attempts clear them. also feared that debris 
the tailwater has plugged coated some the pitot-tube openings and 
may have damaged the entire assembly from impact, which may even more 
severe since the start spillway operation June, 1942. result, 
doubtful that any reliable measurements can taken with this equipment, 
even though carefully planned and all precautions taken assure efficient 
arrangement. was realized that the pitot-tube assembly might inade- 
quate, the salt-injection method will attempted making use electrode 
plates, which were also installed intervals down the spillway section the 
dam. yet suitable methods have been developed measure the water 
surface for the high-velocity flow occurring spillways high dams 
other open channels. 

Difficulties are also encountered when measuring velocities large outlets 
under high heads because the tremendous forces involved, which either 
necessitate extremely rigid instruments prohibit the use instruments en- 
tirely. access also eliminates any measurements, particularly pre- 
vious preparations have not been made for installing the necessary equipment 
during construction. 

The success with which other measurements can made the field, such 
pressures boundary surfaces, vibrations, and air-demand, depends more 
the careful planning the instruments, the ability foresee their work- 
ability, and their installation, than upon the difficulties actually being able 
obtain the measurements. 

Planning and Cost.—Since the success any verification depends upon the 
flow measurements made model and its prototype, compulsory 
formulate plans for the comparison during the laboratory experiments and 
during the design the structure. Such plans should provide for quantitative 
model data covering several flows, not just the maximum, and for the installa- 
tion piezometers, gages, and special instruments the prototype design 
dictated the model tests. result careful planning, the unfortunate 
condition having prototype data one range and model data another 
will avoided, even though might possible regulate the prototype tests 
conform similar ones the model—a condition which rarely exists. 
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The success any planning, moreover, depends upon the circumstances 
under which the laboratory tests are made. Usually tests are performed under 
considerable strain because the limited time available after the preliminary 
design submitted for testing and before the final design prepared for con- 
struction. Discovery faults design during construction also requires 
model experiments made quickly not interfere with the progress 
construction. Because this constant rush, the model data usually 
not cover the desired range, thorough checks the data are not possible, and 
the adequacy the prototype measuring equipment may suffer. 

some flow measurements have been made over period years 
structure never previously tested model, considerable interest 
construct model and secure the data necessary for comparison. This 
reversed procedure has been tried the University Iowa, Iowa City, 
Iowa (see paper Mr. Soucek this Symposium), academic study, and 
was used recently the Bureau Reclamation, but for the primary purpose 
revising hydraulic conditions existing structure. 

The cost taking field measurements will greatly reduced plans pro- 
vide for the installation equipment during construction. Once the proper 
arrangements have been made, the actual cost making the measurements 
will small comparison with the entire cost the structure itself. The 
cost instrumentation, the other hand, relatively high and many 
prohibits any field measurements. 

Properties the Model.—Not all the difficulties making comparisons 
are attributed the prototype tests; fact, most them must-be traced 
directly the model experiments, assuming, course, that the prototype data 
are reasonably accurate. This evident since any comparison actually 
check the model tests, not the prototype. open-channel studies 
small error introduced because usually impossible select scale ratio 
that will permit rigorous similitude between model and prototype roughness, 
the model surfaces being too rough. During experiments different scale 
models sluices the Assuan Dam (4), Egypt, was found that models 
having surface roughness commonly used today agreed with the prototype 
discharge within 2%. When the model surfaces were roughened artificially 
with sand, however, the model discharge decreased 5%. From measurements 
five geometrically similar models spillways, Eisner (5) showed that 
there was definite variation the coefficient discharge with the scale ratio. 
The relative roughness the boundary surfaces was believed the impor- 
tant factor, the effects viscosity being almost negligible. Because the 
uniform variation between the coefficients discharge and the relative rough- 
ness, Eisner predicted the prototype coefficients reasonably well, comparison 
extrapolating the Froude number. 

much greater error introduced model studies closed conduits be- 
cause practically always impossible operate the model Reynolds 
numbers high those obtained the prototype. important realize, 
therefore, that the model friction losses are too large, and that correction 
should made the model data before estimating any corresponding pro- 
totype values. Other errors introduced such fluid forces capillarity and 
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elasticity are negligible the problems usually studied hydraulic labora- 
tories. any event well know the requirements for similitude, how 
apply them, and their effects the results. 

regard the model data, they are usually more accurate than the pro- 
totype measurements, since errors the field are not controlled easily 
they are the laboratory. this regard should determined the 
model and prototype dimensions are geometrically similar, because frequently 


the field structure may not constructed accordance with the dimensions 
used for the model. 


FLOW MEASUREMENTS 


Since necessary take nearly the same measurements the field 
the laboratory, instruments should available for measuring the discharge, 
pressure boundary surfaces, water-surface elevations and profiles, velocity, 
vibration, air-demand, bed movement, and for photography. 
mind the difficulties mentioned obtaining prototype measurements, several 
well-known methods will outlined. Although the instruments used the 
laboratory and the field are basically similar, the prototype equipment must 
designed withstand tremendous forces and convenient for use. With 
little precedent follow, one must rely good judgment for the design the 
instruments and, the same time, able foresee their efficacy. 

flows can usually measured current meters, 
pitot-tube pitometer traverses, weirs, Parshall flumes, floats, and the salt- 
injection method. The particular method used depends the accuracy re- 
quired, the size and type channel, quantity water, and its velocity. 
Accordingly, current meter used most frequently rivers, lined channels, 
Smaller quantities canals are measured with weirs Parshall 
flumes. For high-velocity flow lined channels, the salt-injection method 
used unless the discharge measurement can made some other place where 
more convenient method could adapted. 

Discharges closed conduits are obtained pitot-tube traverses, the 
salt-velocity method, the Gibson method, and venturi meters when in- 
stalled the conduit. These methods are adaptable only plans have been 
made for their use, since has been pointed out that often impossible 
gain access many conduits dams. access impossible, the discharge 
probably could obtained downstream gaging station. 

most these methods calibration the instruments required, and 
debatable whether the calibration some them will the same the 
test site the place calibration. Using rating curve the proto- 
type which has been established from model tests obviously defeats the purpose 
making comparison. 

which are used almost exclusively obtain pres- 
sures boundary surfaces, should installed the prototype positions 
similar those the model. The installation piezometers two the 
102-in. river outlets the dam offered example. Piezometers 
these outlets have been placed the vertical and horizontal diameters. 
Each opening consists high-tension bronze plug drilled with }-in. 
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hole and screwed into the steel plate lining the conduits, with the orifice 
end projecting about in. past the interior surface. After this projection 
had been ground flush with the surface lining, radius was formed 
the orifice end and all burrs were removed. short length }-in. standard 
brass pipe connects each plug steel junction box, where reducing coupling 
used connect each piezometer plug outside diameter 0.035-in. 
wall, soft copper tubing. Junction boxes placed opposite each plug are con- 
nected 2-in. galvanized conduit, through which the soft copper tubing 
collected, and lead six headers recessed inspection gallery below the 
outlets. The headers, one placed above the other, are made 2-in., extra- 
strong brass pipe, with brass air cocks placed rows along the top 
each header make connection with each piezometer line. Each header 
connected high-pressure line one end and bleeder pipe the other 
end, which, turn, connected pressure measuring instrument through 
by-pass. 

similar method could used for measuring pressures concrete sur- 
faces, except the piezometers would attached metal plate set the 
concrete. any event the installation must made during construction. 

measure the pressures, common practice use water mercury 
manometers, and dial gages. Manometers are not always satisfactory because 
their fragility, the fluctuations the mercury columns, and the magnitude 
the pressures frequently encountered. Commercial-type dial gages are usu- 
ally not precise enough. fluid-pressure scale the type used for turbine 
testing proved satisfactory high-head tests the two dams. Pressures 
great 300 per in. can recorded within 0.25 per in. pressures 
100 lb, and within 0.25% any pressure greater than 100 with this 
instrument, which compact and simple construction. When reading 
piezometer pressures, the pressure transmitted copper tube grid mounted 
gimbals weighing beam the base the instrument. The weighing 
piston, which lapped fit the cylinder, restrained from being forced 
upward. The resulting downward reaction the weighing cylinder trans- 
mitted, therefore, the lower weighing beam, and measured upper scale 
beam mounted above the base. The weighing piston rotated small 
electric motor while readings are being made avoid static friction with the 
walls the cylinder. 

Pressure cells used combination with group piezometers and 
oscillograph may used, but the cost such equipment prohibitive 
many cases unless available from previous experiments. This method 
has one advantage that continuous record the pressure fluctuations 
recorded the film sensitized paper the oscillograph. 

making any pressure measurements either the model prototype, 
the care with which the piezometer openings are installed cannot over- 
emphasized. Any burr will affect the reading seriously, and excessive surface 
roughness either upstream downstream from the piezometer opening will 
influence the measurements. essential have high-pressure water 
line for flushing the piezometer lines clear air bubbles, and the tubing from 
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the piezometer the measuring instrument should alined, possible, 
eliminate any bends sags that might trap air. 

Water these measurements are not important the 
discharge and pressure, nevertheless check desirable such design assump- 
tions freeboard and roughness coefficients. connection with the work 
the Society’s Special Committee Hydraulic Research, the problem 
entrainment high-velocity flow requires, part, measurement the water 
surface open channels field structures. comparison the low-velocity 
flow models where the water surfaces are relatively smooth and well-defined, 
there actually definite water surface for the high-velocity flow field 
structures, but only mass spray and aerated water. Efforts point-gage 
such water surface narrow wasteway the Yakima project, Washington, 
were successful measure, but only water surface was actually 
obtained. 

Although point gages may useful for measuring water surfaces narrow 
channels, might more feasible wide channels use system cables 
manipulate heavy target, which could located transits. profile 
the water surface could also obtained suspending plumb bob from the 
training walls, from grid painted the training walls; but such measure- 
ments would affected the wall friction. 

For measuring water-surface elevations reservoir, use has been made 
the familiar staff gage and several commercial instruments. Those based 
the principle long air column transmitting the pressure should 
avoided because the effect temperature change the air column, and 
because possible air leaks the tubing. 

methods available are directly associated with discharge 
measurements, use made current meters, pitot tubes, salt-injection 
method, and floats. Although somewhat unsuccessful Grand Coulee Dam, 
the installation pitot tubes, discussed herein, may prove successful under 
different circumstances. The accuracy required for model-prototype com- 
parisons eliminates the use surface and subsurface floats. 

Vibration.—It frequently desirable measure the vibrations occurring 
prototype structure, not necessarily compare with the model establish 
relation between the two, but more ascertain what effect the vibrations 
may have the structure itself. The flutter nappes discharging over the 
lips drum gates the crest movable weirs dams one example such 
vibration (6). The discharge large valves and regulating gates also 
frequently accompanied vibrations. Accelerometers and pressure cells were 
used with oscillograph measure vibrations due fluttering nappes 
dam another western state. The instruments were designed the principle 
that change resistance occurs with change pressure pile carbon 
disks. The accelerometers measured the accelerations three directions 
the dam, and the pressure cells were used measure the amplitude and fre- 
the inducing forces exerted the dam the flutter the nappes 
spilling over the lip drum gates. The instruments were calibrated noting 
the deflections produced the moving elements the oscillograph under 
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various conditions created shaking table. Less complicated instruments, 
such vibrating reeds, which give only the frequency, have been used else- 
where, the amplitude being determined magnifying the vibration and using 
oscillograph. 

often necessary and desirable measure the quantity 
air required prototype relieve sub-atmospheric pressures under 
nappes, conduits downstream from needle valves and regulating gates, and 
regulating gates and valves themselves. difficult problem deter- 
mine the size air ducts required the prototype model tests, because 
the uncertainty the similitude involved; thus any check the prototype 
will aid considerably establishing the transference relations. measure 
air quantities the prototype (and the model) anemometer, orifice plate, 
pitot tube could installed the air ducts. Although not entirely suitable 
for prototype measurements because the large quantities involved, 
rotameter may used for model tests. 

Bed Movement.—The scour river channels downstream from dams and 
spillways obtained surveys soundings. These measurements, 
however, can used only for qualitative comparison. 

Photography.—Still and motion pictures all flow conditions observed both 
laboratory and field work are invaluable for purposes comparison. Fre- 
quently startling similarity seen between model and prototype photo- 
graph some flow phenomenon. For still pictures, the photographic equip- 
ment should include least standard view camera and miniature camera, 
both with fast lens. The motion-picture camera should equipped for slow- 
motion. The higher the number frames per second, the better, because 
many prototype flows cannot studied unless considerably slowed down 
the motion pictures. Telephoto lens attachments are considered essential for 
both the still and motion-picture cameras. 


EXAMPLES QUALITATIVE COMPARISONS 


relatively small irrigation structures diver- 
sion dams, check drops, wasteways, and turnouts offer opportunities for making 
model-prototype comparisons because their continuous operation various 
discharges during irrigation season. The primary concern these struc- 
tures prevent scouring the canal downstream using hydraulic-jump 
stilling basin wherever possible. Consequently, sufficient make only 
qualitative comparison any erosion that occurs. Most the small irrigation 
structures use today are based designs developed over period years, 
that model studies have not been utilized except for revising some existing 
structure inadequate design. Such condition required laboratory in- 
vestigation Check Drop the Main Canal western project. 

The gradual growth irrigation, the Main Canal 
had increased capacity. maintain proper gradient the 
canal, twenty-three check drops were installed intervals about one and 
one-half miles. These drop structures, which were built during the period 
were similar design, but decreasing size downstream the 
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discharge the canal became less due deliveries. Each structure consists 
rectangular concrete basin divided into bays concrete piers which are 
surmounted steel brackets. placing flashboards between each pair 
brackets, the depth flow between successive drops can regulated. Wing 
walls each end the concrete basin extend into the canal banks right 
angles. The drop water surface these structures varies from in. for 
the maximum in. for the minimum, with maximum discharge 1,300 
per sec the first drop gradually decreasing 520 per sec the 
last one. 

_Not many irrigation seasons had passed before excessive scouring developed 
the canal immediately downstream from the drops, particularly the larger 
ones. All attempts check the erosion failed, until finally the width the 
canal below the drops was increased much 50% some cases. 1938 
the problem eliminating this scour was assigned the hydraulic laboratory 
the Bureau Reclamation, Denver, Colo. 

Model the drops were similar design, Check Drop was 
selected for study, and scale model was constructed. The initial tests 
the laboratory demonstrated that the model would reproduce the flow con- 
ditions and scour existing the field (Figs. and Further analysis 
disclosed that, because the small drop water surface, standing waves de- 
veloped instead more efficient energy dissipator, such the hydraulic 
jump. high surface velocities proceeded downstream the center 
the canal, thereby producing large return eddies along the canal banks which 
eroded them severely. 

satisfactory solution this problem was not readily apparent, for not 
until many designs had been tested was solution found which would eliminate 
the unfavorable conditions completely. This solution consisted adding 
curved training walls the entrance the drop, extending vertical walls and 
concrete floor downstream form stilling pool, and placing deflector 
between these walls directly above the concrete basin force the flow under 
the tailwater. These revisions created more uniform flow distribution 
through the structure, and provided energy dissipation sufficient eliminate 
the high surface velocities which had caused side eddies scour the canal banks. 

Comparison.—After the prototype had been revised, observations its 
performance showed that was acting manner similar that predicted 
the model tests. Figs. 1(c) and 1(d) show the revised design the model and 
prototype. The collection debris the prototype structure evidence 
roller action the stilling pool. This indicative energy dissipation 
throughout the entire depth, instead only the surface noticed Figs. 
and 

During the laboratory studies similar collection debris was noticed, but 
was given little consideration since interfered with photographing the 
action the model. the field, however, has been considerable value 
clearing the canal large quantities debris, conveniently removed the 
drop structure. Three miles downstream Drop (which was also revised), 
there little trash remaining the canal. 
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The results this comparison are considerable satisfaction those 
familiar with the problem, but previously mentioned such correlation 
more convincing those directly concerned. Nevertheless, this study 
offered convincing proof that model tests are excellent tool for the 
hydraulic engineer, and that they can relied upon high degree. 


EXAMPLES QUANTITATIVE COMPARISONS 


The 102-in. outlets were provided three tiers 
each the dam release total discharge 225,000 per sec 
Each outlet has bellmouth entrance followed length straight 


Top Dam, 1311.08 


135'x28' Drum Gate 
10n 


102" Ring Seal Gates 


Outlet Works 
Trashrack 


1036.67 


Paradox Gate 


102" Diameter Outlet Conduit 


conduit, which terminates elbow and cone the downstream face 
the dam. outlet closed service gate and emergency gate, 
partial gate openings never being used with ring-seal paradox gates. Except 
for part the lower tier, each outlet has steel-plate lining. The three 
tiers outlets El. 1136.67, El. 1036.67, and El. 935.76 are designated 
the upper, intermediate, and lower outlets, respectively. Since the outlets 
occur pairs each tier, centers, further designation given east 
(E) and west (W), and the block number which they occur. Fig. for 
example, shows upper outlet 63W discharging down the face the dam. 
During the design the outlets, plans were formulated for installing 
piezometers one intermediate outlet, 51E, and one upper outlet, 51E, 100 
immediately above. Piezometers were not placed the lower outlets because 
they were not completely steel-lined and because detailed model studies were 
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not made them. locations similar those the model, series eleven 
piezometers have been placed the bellmouth entrances along both the invert 
and crown (vertical center line), and the sides (horizontal center line). 
Similarly, the horizontal part the two conduits contains twelve piezometer 
rings the intermediate outlet, and ten the upper outlet; both have ten 
rings the elbow and cone. The installation methods have been described 
and the vertical section intermediate outlet shown Fig. indicates 
the distribution piezometers. 

Model résumé the 1:17 scale model tests the outlets 
the dam has been given another paper tests were concerned 
with the bellmouth entrance, the profile the outlets, the elbow and cone, 
and the deflectors placed the spillway immediately above the exit 
the upper and intermediate outlets, which can seen Figs. and 4(a). 
Measurements the model were concerned primarily with the pressures the 
boundary surfaces, measured with water and mercury manometers, and the dis- 
charge for several heads, obtained from weir measurements. 

Only the model data from the intermediate outlet tests will used for 
comparing with the prototype. Some tests were made the laboratory 
the lower outlets, but piezometers were installed the prototype, and 
since the upper outlets are quite similar those the intermediate tier, 
additional tests were not required. 

Prototype tests were made the intermediate and 
eight the upper outlets. Heads above the exit ranged from 234 


“ay, 
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for the intermediate outlet, and from 133 for the upper outlet. 
addition measuring pressures with the fluid pressure scale, the reservoir 
elevation, water temperature, tailwater elevation, discharge, and the number 
outlets open each tier were recorded. make comparisons with the 


Downstream Face Dam 
Elbow and 


Deflector 
Axis Dam 
y2 °19' —Trough 
(a) VERTICAL SECTION inside Diameter 
36.67 


model the intermediate outlet, was possible use only those tests which 
the adjacent outlet, 51W, and the upper outlet, 51E, were closed. This not 
only conformed with the laboratory arrangement, but also eliminated from 
consideration any secondary effects the pressures. These effects were 
discovered first during separate model investigation the bellmouth, and 
later during the prototype tests when was noticed that pressures the bell- 
mouth were slightly lowered when the adjacent outlet 51W was operating, 
while pressures the elbow and cone were raised when the upper outlet 51E 
was operating. 

Notation.—The letter symbols used this paper are defined where they 
first appear and are assembled for convenience reference under ‘‘Bibliography 
and the end the Symposium. 

Model-Prototype Pressure first comparison given 
Fig. intended show the characteristics the hydraulic gradient. Be- 
cause the field tests could not always made when the reservoir reached 
corresponding elevation, the model and prototype heads shown are not 
parable. Upon first consideration would seem desirable for the heads 
agree, but even they did, agreement between the model and prototype 
pressures would lacking. The reason for this evident from Fig. where 
may seen that the model tests were made much lower Reynolds num- 
bers (R); hence the friction losses 
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were proportionately greater coefficient friction; loss head due 
friction, feet water; average velocity through conduit; length 
gravity). Consequently, when the model pressure heads are scaled the 
prototype the Froude law, they will lower than corresponding proto- 


Values 


Reynolds Number, 


FRICTION 


type values. This fact generally well-known and only mentioned show 
the effect not satisfying the similitude requirement that the model should 
tested the same Reynolds number the prototype. This condition 
incompatible with hydraulic models based Froude’s law, addition the 
fact that the extremely high velocities required comply with the Reynolds 
law are impossible obtain many laboratories. any event the model 
results are the side safety. 

Notwithstanding the lack complete similitude, study Fig. indicates 
certain similarity between the model and prototype gradients. Whereas 
the prototype data are quite consistent, the model data appear erratic 
some the piezometers. 

Pressure Drop make more valid comparison than the 
one Fig. the pressure drop the bellmouth has been analyzed. 
choosing short length conduit, the effect the differences friction 
negligible, and using dimensionless parameters, comparable heads are not 
required. Accordingly, shown and curves were drawn for the 


crown and invert using ordinates and abscissas which 


Lu 
the head above El. 1036.67; the pressure any piezometer the bell- 
mouth referred El. 1036.67 datum; the pressure the reference 
piezometer, number 11, referred the same datum; the distance from 
the face the bellmouth any piezometer; and the distance from the 
face the bellmouth piezometer number 11, shown Fig. 


values express the pressure drop any piezometer percentage 


the total, whereas the distance ratio eliminates any small differences 
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between the location piezometers the model prototype. Two the 
prototype tests were omitted from the diagrams not confuse them with 
more superimposition the points, but all the model values are included. 
The actual pressures the piezometers the bellmouth were all positive 
(above atmospheric) except for low heads, and these did not exceed vacuum 

The prototype curves illustrate that the percentage pressure drop prac- 
tically constant nearly all the piezometers regardless the head. This 
not true for the model tests, even though the average values the model are 
within approximately the prototype. Reasons for the dispersion the 
model data may explained experimental error, different approach con- 
ditions, and unstable flow the model lower Reynolds numbers, although 
above critical Reynolds numbers. 

Experimental errors are probably due the installation the model 
piezometers, even though great care was taken install them properly. 
Nevertheless, any minute obstructions the boundary surfaces and the 
piezometer openings have serious effects. The laboratory experiments, more- 
over, were made under considerable duress enable the design keep ahead 
construction the field. Consequently, the customary close checks 
the data could not made complete they would have been under more 
normal procedure. The different approach conditions refer the absence 
the laboratory trashrack upstream from the bellmouth. However, was 
observed from separate tests the bellmouth that the trashrack effects were 
negligible. 

Evidently the differences friction play small part even though com- 
parison confined the bellmouth alone. estimate can made readily 
the relative roughness either the model prototype surfaces, but since 
the type surface each nearly the same, the relative roughness the 
model would greater. Consequently lower Reynolds numbers, the 
resistance factor for the model will vary considerably more within its range 
tests. can assumed further that the model tests were for rough flow, 
whereas the prototype tests were more nearly smooth flow. Finally, might 
added that transition from laminar sublayer turbulent sublayer will 
develop less rapidly the model, thus influencing the velocity distribution 
and pressures. reasonable assume, therefore, that the model flow 
not sufficiently stable give the consistent data obtained for the prototype. 
Accordingly, believed that the laboratory tests should have been made 
greater velocities (higher Reynolds numbers), disregarding heads based 
the scale ratio the model. Had this been done, the model data analyzed 
Figs. and would have been more consistent and would have approached 
the prototype values more closely, although certainly good agreement 
for the average curves. 

Pressure Drop Elbow and curves have been drawn 
Figs. and show the pressure drop the elbow and cone manner 
similar that described for the bellmouth comparison. eliminate the 
differences friction the elbow, the pressure drop has been expressed 
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datum, the start the elbow and the end the cone, respectively; 
the pressure any piezometer referred the same datum; the 
distance along the center line from piezometer any piezometer; and 
the center-line distance from piezometer piezometer (see Fig. 12). 
Figs. and two prototype tests have been omitted for clarity, but all model 
tests are shown. The heads, however, are referred El. 1021.67 instead 
El. 1036.67, and average center-line pressure drops have been added. The 
reason for negative ordinate values for the pressure drop the crown due 
increase pressure from piezometer 23, thus making the numerator negative 
the ratio Although not shown Figs. and 11, the center-line 
pressure drop falls good average between the curves for the crown and 
invert for both model and prototype. This may seen from the curves given 
Fig. 

Figs. and show the data for both the model and prototype almost 
equally consistent, the model values still being slightly erratic. Reasonable 
agreement, however, lacking. general, the pressure drop indicated from 
the model tests about 30% greater than shown 
from the field measurements. This unsatisfac- 
tory agreement may explained the principal 
reasons applied the bellmouth comparison; 
namely, experimental errors and unstable flow for 
the range Reynolds numbers the model. 

computing the model data for the elbow 
and cone, was evident from the excessive disper- 
sion the points that some the pressure read- 
ings were error. adjust these, curve 

Fic. head versus pressure was plotted for each piezo- 

meter. Since these curves should straight 

line was possible detect which readings were erroneous and adjust 
them the line. Unfortunately, the pressures the start the elbow, 
and the end the cone, are not initially correct are not properly 


adjusted, then the total pressure drop causes error the ratio 


plotted Figs. and 11. 


One the field tests was similarly adjusted where obvious errors were 
noted. checking the piezometer numbers attached each header, was 
found that the erroneous readings were all group located one piezometer 
header. was concluded that the pressure from another piezometer the 
same header was being registered the time this group was being recorded. 
error this nature easily made when many piezometers are required, 
unless particular care taken certain that only one piezometer being 
read time. 

Since the length the elbow and cone nearly six times that the bell- 
mouth, evident that the differences between the model and prototype 
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friction will have still more effect the comparison, even though study 
the losses shows the bend loss centrifugal effect considerably greater. 
The latter seems noticeable the prototype curves Figs. and 11, 
where the points are more than for the bellmouth measurements 
shown Fig. tests other structures, known that slight 
change diameter the cone exit will cause appreciable changes the pres- 
sure gradients through the pen- 

stock, making more 

did not permit 

more complete 

the data for the elbow and 
which might more 
valid method for comparing the 


Downstream Face Dam 


Steel Lining 
model and prototype dis- 


and Are Measured 


Fig. Curve based Along Center Line 


the model tests for 
Curve the model rating curve 
corrected take into account the difference between model and prototype fric- 
tion. make this correction, was first necessary estimate average fric- 
tion for the prototype outlet. This was not readily accomplished 
because lack reliable friction-coefficient data Reynolds numbers 
high was estimated that value 0.0085 for the 
friction factor would reasonably accurate. The model friction coefficients 
and Reynolds numbers were then computed from the laboratory pressure tests 
and plotted Fig. Since the model friction too large, the model heads 
for given discharges would too high. Accordingly, assuming that all losses 
may sealed the prototype except friction, the model heads were reduced 


72 


the value (fm which was obtained from Fig. and 


0.0085 estimated. After the pressure data for the prototype became 
available, the friction coefficients and Reynolds numbers were computed and 
plotted Fig. The average value for was 0.0095 compared 
estimate 0.0085. Two 40-ft reaches were used for the prototype computa- 
tions, starting distance approximately eight diameters downstream from 
the entrance the outlet, corresponding reaches being used for the model 
computations. That this probably not sufficient length develop 
normal velocity distribution seen Fig. where the upper points were 
obtained using the first 40-ft reach, whereas the lower points were obtained 
from the following 40-ft reach. 

The discharges given Curve Fig. 13, which were used for computing 
the prototype friction factors and Reynolds numbers, were developed from 
prototype pressure measurements, treating the bellmouth submerged ori- 
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fice. The average pressure heads measured piezometer ring enabled the 
values computed (see Fig. The four piezometers this 
point are located in. from the entrance and the tangent between the 
bellmouth and transition immediately downstream, shown section Fig. 
Using the area the bellmouth this point, and assuming coefficient 
discharge, the nozzle formula 


was applied. Similar results would obtained using and the area the 
start the bellmouth, but with 0.60, applying the same area and 
coefficient with the head the bellmouth reduced the losses the bell- 
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Fic. 13.—CoMPARISON DISCHARGES 


mouth instead These computations are based the premise that the 
bellmouth fits the jet from sharp-edged orifice placed the plane the 
upstream face the dam. 

Curve Fig. 13, was obtained from discharge measurements river 
gaging station one- mile downstream from the impossible 
pass the flow the river through only one outlet, the discharge appli- 
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cable single outlet had computed indirect First, 
rating curve was obtained for one lower outlet before the water surface 
the reservoir reached the intermediate outlets. This was accomplished 
plotting effective head (reservoir elevation minus tailwater elevation the 
gaging station) against the discharge for one lower outlet (total measured 
discharge divided the number lower outlets operating). soon the 
intermediate outlets began operate, the total river discharge increased ac- 
cordingly, but knowing the effective head the lower outlets and the number 
operating, was possible obtain the discharge per outlet from their rating 
curve, and thus the total discharge contributed the lower outlets. Sub- 
tracting this value from the measured flow the river, and dividing the 
number intermediate outlets operating, gave the discharge for one outlet 

This method assumes that each outlet will discharge the same amount, 
and that the jets from the intermediate outlets will not affect the effective 
head the lower outlets. Although these assumptions are not strictly cor- 
rect, believed that Curve Fig. 13, and certainly Curve would agree 
within any rating curve obtained measuring the discharge only 
one outlet. 

comparing the four rating curves over range head from 230 
ft, the uncorrected model discharge (Curve from 10% 21% less 
than the prototype discharge computed the bellmouth (Curve 3), whereas 
the corrected model discharge (Curve only less. Similarly, 
the uncorrected model discharge (Curve from 11% 27% less than the 
measured prototype discharge (Curve 4), whereas the corrected model dis- 
charge (Curve only 15% less. 

Summary.—The agreement between the model and prototype this quan- 
titative comparison reasonably good except one instance. Considering 
experimental errors and lack similitude relative Reynolds’ number, the 
pressure drop the model bellmouth agrees with the prototype within 5%; 
the pressure drop the elbow and cone the model, the other hand, 
greater than the prototype; and, finally, the model estimate the pro- 
totype discharge is, average, low. Even though complete accord 
could not obtained, the scale model estimates are the safe side and 
the field measurements show that excellent prototype design has been 
evolved. follows, therefore, although complete similitude was lacking, that 
the model could relied upon prove not only the inadequacy several 
proposed designs, but also develop the satisfactory design finally adopted. 


CONCLUSIONS 


make better model-prototype comparisons, evident that plans 
should prepared for procuring field measurements during the model tests 
and while the prototype being designed. This not only enables sufficient 
range model data obtained, but also makes certain that the in- 
struments required the prototype will located properly and incorporated 
part the design the structure. Even these precautions may not 
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yield the desired results because the difficulties actually obtaining proto- 
type flow measurements. Thus, the reasons for few quantitative compari- 
sons are traced the necessity waiting years before sufficient water 
available, the difficulties measuring large discharges and high velocities, 
the fact that the measured discharge derived from several sources in- 
the one used comparison. Finally, the model experiments 
themselves introduce objectionable factors influencing comparisons, the chief 
one being the inability test the model conduits the correct range 
Reynolds numbers. any event important build hydraulic models 
large scale ratio possible. 

The qualitative and quantitative examples model-prototype comparisons 
given this paper reaffirm the high degree which the model can relied 
upon. Although the avowed purpose any comparison, particularly the 
quantitative type, should affirm the theory hydraulic models 
modify the transference relations necessary, cannot accomplished from 
only few comparisons. present, the emphasis should laid the 
limitations and accuracy quantitative estimates based model tests. 
this regard, equally important show both inferior and good model- 
prototype comparisons. The lessons learned therefrom will certainly demon- 
strate the adequacy present-day model research. 
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DRY-DOCK SUCTION CHAMBER 


Other papers this Symposium describe studies which the analyses 
model and prototype data could based upon direct quantitative comparisons. 
This paper will devoted study wherein, due the lack certain 
measurements the prototype, the comparison model and prototype per- 
formance was approached from inductive standpoint. The problem under 
investigation was the unsatisfactory operation pumps used drain large 
dry dock; the medium study was 1:10-scale model constructed the 
Waterways Experiment Station. 


The dry dock, together with its culvert system and suction chamber, 
shown the general plan drawing Fig. passages 
the suction chamber from adjacent dry docks are not shown.) The suction 
chamber shown greater detail Fig. 15, and the arrangement pump 
intakes (1, and and sluice gates (3, and indicated. This arrange- 
ment, the high rate discharge (130,000 gal per min maximum per pump), 
and the limited dimensions the suction chamber, gave rise unstable and 
sometimes violent flow conditions, which were manifested the pumps 
vibration, loss capacity, and symptoms cavitation. (All such objection- 
able features will covered the general term ‘‘rough Test 
runs the prototype yielded notes and observations the character 
pump performance, but did not provide any numerical data evidence 
the character flow within the suction chamber. 

constructing model for the investigation, was essential that the proper 
distribution flow through the three sluice gates assured. Accordingly, 
the entire system side culverts and inlets was reproduced the model. 
was not feasible reproduce the prototype pumps means miniature 
pumps homologous series; however, the proper metered discharges were 
induced through the intakes ordinary centrifugal pumps. For purposes 
visual observation, the suction chamber and pump intakes were made trans- 
parent pyralin. 

gain some idea flow conditions within the suction chamber, com- 
pressed air was injected into the model under operating conditions similar 
those known give unsatisfactory pump operation the prototype. From 
observations the air bubbles these preliminary tests was learned (a) that. 
flow entering the suction chamber from the sluice gates created initial rotation 
about one another the pump intakes (depending upon the particular 
conditions operation), and (b) that initial vortices formed were accentuated 
censors.) 
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the water was drawn upward into the bellmouth intakes. Under such con- 
ditions was reasonable expect that the pressure the axis any intake 
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affected (that is, the center the vortex) would frequently fall below the 
suction lift capability the pump. was concluded that the rapid and vio- 
lent fluctuations pressure attending vortex action within the intakes were 
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responsible for alternate breaking and restoring suction head the 
pumps, resulting the aforementioned symptoms. 

support this conclusion, and provide basis for the comparison 
model and prototype performance, was desired that record made the 
amplitudes and frequencies pressure fluctuations within the pump intakes. 
was believed that such data would permit the establishment criterion 
for pump operation, whereby rough pump operation would revealed 
record showing fluctuations great amplitude and frequency, and whereby 
smooth operation would indicated slight variations pressure. From 
such criterion, the model reproduction prototype performance could 
confirmed, and the effectiveness any corrective measures under trial could 
evaluated. 

obtain record pressures within the suction bells, the instrument 
shown Fig. was devised. The operation this instrument very simple: 
Pressure impulses originating the pitot-tube opening are transmitted through 


Leads 


the confined liquid the tubing the face the pressure cell; the resulting 
deflections this diaphragm are recorded electrically oscillograph give 
continuous time-record the instantaneous values pressure. 

Using this instrument the model was run under each operating condition— 
combination pumps and sluice gates operation—for which observations 
had been made the prototype. Oscillograph records taken during few 
such tests are plotted Fig. 17; and for each the test conditions the observa- 
tions made under comparable conditions the prototype are recorded 
Table (In Fig. 17, and all other figures, dimensions are expressed terms 


prototype equivalents.) comparison model and prototype observa- 

tions may seen that they are consistent all cases. Furthermore, the 
criterion for rough and smooth pump performance definitely established. 
From such model-prototype comparison was safe assume that predic- 
tions the effectiveness corrective measures could relied upon. 
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was considered that the unsatisfactory performance the pumps would 
improved the vortices the intakes could broken up. For 


(Discharge 130,000 Gal per Min) 


Sluice gates 


Pump op- 
open 


erating Corresponding prototype observation 


Extremely bad operation 

Generally unsatisfactory operation 

With more than one pump operating Sluice Gate the farthest 
unit could not operate; generally unsatisfactory operation 

Under initial conditions Pump operated satisfactorily 

Operation very rough 


When Pump was shut down, power input Pump decreased. 
When Sluice Gate was opened power input Pump decreased again and its operation was good. 


Initially (see Figs. 17(g) and 17(h)) 


umps and were each discharging 130,000 gal per min with Sluice 
Gate open. Then (see Fig. 


was shut down and Sluice Gate was opened. 


pose, vertical vanes the design shown Fig. were installed the model 
pump intakes. order that flow longitudinal direction unobstructed, 
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transverse vanes were not extended the floor the model suction chamber. 
the prototype, however, all vanes were extended the floor. The results 
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model test run with these vanes place are shown Fig. 19. com- 
paring these results with those Fig. 17, the improvement evident, and 
smooth operation the pumps indicated. Modifications the basic design 
vanes were also tested, but improvement over the basic design was found. 

Accordingly, vanes the basic design were adopted for installation the 
prototype. final step the analysis, the model predictions were thor- 
oughly substantiated the prototype; reported that operation the 
pumps, following installation vertical vanes the intakes, has been entirely 
satisfactory. complete model-prototype confirmation would include the 
measurement pressure fluctuations the prototype but thus far such mea- 
surements have not been practicable. 
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INTERESTING LOCK CONSTRUCTION 


JAMES HARTIGAN,® Esq. 


lock and dam, one the major projects constructed eastern 
river the Tennessee Valley Authority (TVA), was designed 1934. 
model the proposed lock was built and tested the Engineer Sub- 
Office, Iowa City, Iowa, for the purpose aiding the hydraulic design 
the filling and emptying system. the Tennessee development con- 
tinued, other locks were designed and models were used study special 
problems and check the final plans. 

The question the reliability hydraulic lock models simulating condi- 
tions which would exist the full-size structures developed with the increased 
use this method design. Accordingly, steps were initiated test model 
existing lock which exact geometric similarity could obtained. 
The project chosen for the correlation was lock, the first project completed 
the river the TVA. model the lock was built December, 1936, 
and comparisons were made between model tests and limited amount 
prototype test data available that time. The results these comparisons 
were satisfactory, but the prototype data available were not adequate 
warrant drawing definite conclusions. 

The division engineer the Ohio River Division initiated steps for testing 
number prototype locks two eastern rivers verify the model tests 
and study, under actual operating conditions, the hydraulics various 
modifications the port and culvert system. Tests were made four locks 
one river, and another eastern river lock. Model tests had been made 
only three locks, and this paper will confined discussion tests 
made these three locks and their models. 


The locks question are all equipped with the same general type filling 
and emptying system. Longitudinal culverts located the side-walls the 
lock open into the upper pool through gallery intakes and into the lower pool 
through similar discharge sections. Short lateral ports connect the culverts 
with the lock chamber near the floor. The flow each culvert controlled 
valves upstream from the first and downstream from the last lateral port 
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the chamber. filling the lock chamber the downstream valves remain 
closed, the upstream valves are opened prescribed rate, and water from the 
upper pool enters the culverts through the intake galleries and flows into the 
lock chamber through the short lateral ports, the flow water continuing 
until the water levels the chamber and upper pool are equalized. empty 
the lock the upper valves are closed, the lower valves opened, and water the 
chamber discharges into the lower pool until their levels coincide. 

Fig. shows the essential details the hydraulic system one lock. 
The intake galleries and culverts are designed effect uniform accelera- 
tion the flow possible from the face the lock wall the valve section. 
The ports (see Fig. 21) were streamlined with well-rounded entrance corners 


HORIZONTAL 
SECTION 


VERTICAL SECTION (b) DISCHARGE PORTS 
(a) FILLING AND EMPTYING PORTS 


eliminate separation the jets from the leeward walls. The bores were 
tapered effect greater discharge capacity and reduce turbulence 
the chamber expanding the jets and decreasing their velocity. The dis- 
charge section was also designed with view obtaining uniform distribution 
discharge the lower lock entrance and minimize turbulence the 
immediate area. 

Another lock has system essentially identical, only enlarged scale 
approximate proportion the relative dimensions the structures. 

The hydraulic system still another lock typical the designs used 
inland waterways locks prior the use models analyzing the hydraulics 
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involved their operation. Deviation from perpendicular lines effect 
better hydraulic efficiency was practically avoided. Simplified construction 
was given primary consideration the expense efficiency throughout 
lifetime operation. 


THE TESTS 


The parallel tests that were made corresponding model and prototype 
locks will described general terms order eliminate the necessity 
discussing the experimental work performed each individual structure. 
Each type test considered was performed least one pair similar locks 
and most cases all three pairs considered this paper: 


(a) The stages near both ends the lock chamber were observed with 
respect time filling and emptying operation. From these data were 
plotted filling emptying curves, rate-of-rise rate-of-fall curves, and 
water-surface slope curves, all with respect time. Comparisons model 
and prototype were made the basis total time required complete the 
operation, the coefficients discharge for the entire hydraulic system, the 
maximum value the rate change stage the lock chamber, and the rate 
decrease the flow after the maximum had been reached indicated 
the slopes the rate rise fall curves. 

(b) Velocities were observed the ports the intake section, the mani- 
fold section, and the discharge section, and from these data comparisons the 
distribution discharge the respective sections were made. 

(c) Pressure measurements observed the ends the culvert transition 
sections were used conjunction with volumetric discharge measurements 
the lock chamber determine discharge calibrations and venturi coefficients. 


The primary use for hydraulic models provide data for design. 
such, the functions model are generally completed before the prototype 
built, and not infrequently the final design modified from that used the 
model, due contingencies which had not arisen until the final stages the 
design were reached and which were considered too insignificant warrant 
the expense checking the model. result, very few prototypes are 
built exactly like the models from which their designs were developed. Such 
differences were found the locks covered this study, but wherever possible, 
corrections will applied dimensionless parameters used validate the 
desired comparisons. 


CoMPARISON PROTOTYPE AND 


Intake primary consideration the design the intake sections 
one lock was uniform distribution flow the ports. With flow 
one direction only, tapered culvert section could used. That the design 
successfully approximated uniform distribution the four ports one side 
the culvert, and that agreement existed between model and prototype dis- 
tributions, are indicated the graph Fig. 22. The discharge the model 
intake ports varied uniformly from 21.5% 28.5% and the corresponding 
prototype ports from 22% 27%. 
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Transition rating equation for culvert transition section, 
expressed terms culvert discharge and pressure head differential between 
two sections based the Bernoulli theorem; thus: 


The coefficient discharge for the transition section meter equal the ratio 
the observed the theoretical rating. comparison prototype and 


Percentage Total Flow 
uo 


Port Number 


model locks made, shown Table the basis the performance 
their culvert transition sections flow meters. The lack exact geometric 
similarity the transition sections model and prototype accounts for the 
small difference the constants 
the theoretical rating equa- 
tions; however corresponding 
equations and coefficients for 


model and prototype are vir- 
Coefficient 
Observed Theoretical Lock Stage 


tervals during filling operations 

the model and prototype are 
shown Fig. 23. The points defining the filling curves are averages two 
stage measurements, one near each end the lock chamber, observed each 
time interval. 

The filling curves for the model and prototype locks indicate considerable 
difference the total time operation and the form the filling curves. 
There were two factors which materially affected the results these tests— 
one difference construction the lock-chamber ports (Fig. 24), and the 
other difference the time opening the culvert valves. the model 
the total cross-sectional area the lock-chamber ports the smallest section 
was 125.0 per culvert, whereas the prototype the corresponding areas 
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were 145.0 the river wall and 147.8 the land wall. The culvert 
valves the model were opened time corresponding min the proto- 
type, whereas the filling test shown the right valve the prototype required 
min sec open, and the left valve min sec. third factor which 


Rate Rise, Feet per Minute 


Filling Curve 

Water Surface Slope 

Lift 55.0 61.0 


Minutes and Seconds 12:00 10:10 


Time of Opening Right Valve 
Time of Left Valve 


5:04 2:00 


a 


Water Surface Elevation, Feet 


ny 


Time 


CHARACTERISTICS 


also affected the filling time was difference the lifts. the model lift 
was simulated, whereas the prototype was the maximum differ- 
ence water levels existing during the period the tests. 

the two filling curves shown Fig. were matched the instant when 
the corresponding locks were full, noted that they would follow approxi- 
mately the same path for the period after the valves were completely open. 
The larger throat area the prototype lock-chamber ports, which tended 
increase the efficiency the prototype over that the model, appears have 
been offset for the most part the larger and more efficient divergence angle 
the model ports. 

The maximum rate rise fall the water level lock occurs shortly 
after the culvert valves are fully open. Because the difference time 
valve operation, the rate-of-rise and rate-of-fall curves for the model and proto- 
type one the locks differ considerably. The maximum rate rise the 
model was 11.1 per min min sec, and the prototype reached 
9.0 per min min sec. manner similar the filling curves, the 
rate-of-rise curves for model and prototype follow nearly parallel paths after 
the peak rate has been reached; the average slope the curve for the proto- 
type 1.40 per min?, while that for the model 1.24 per 
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The differences construction details and valve operation, noted the 
discussion the filling tests, apply also with respect the emptying tests. 
However, the larger port area and smaller lift existing the prototype, tending 
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reduce the emptying time, were practic- 
ally offset the longer valve period. The 
emptying time the prototype was min 
sec, whereas the model was min 
sec. The dissimilarities construction 


details the two structures make 


sons the filling and emptying charac- 
teristics unsatisfactory. 

Lock-Chamber tests have 
indicated that longitudinal slope the 
water surface the lock, such was cre- 
ated surges the energy change from 
kinetic potential accompanying the de- 
crease velocity flow from one end 
the lock the other, the primary factor 
that influences movement vessel during 
lockage. Improper operation the valves 
and uneven distribution discharge the 
ports side culvert system will promote 
surging the lock; and, once started, surg- 
ing will continue throughout 
part the filling period. important, 
therefore, that the filling controls oper- 
water the beginning the filling period, 
and that the system designed dis- 
tribute the lockage inflow evenly among the 


lock-chamber ports. 

From the port velocity tests made the model and prototype locks, 
distribution-of-flow graphs for the ports the right river wall for the 
filling operation were plotted shown Fig. 25(a). The flow distribu- 
tions the model and prototype differed considerably, due largely the dis- 
similarity the ports. earlier laboratory experiments was found that 
the distribution discharge from the ports during filling operation was 
influenced large extent the value the cumulative cross-sectional area 
the lock-chamber ports the smallest section. When the port area exceeded 
certain value, the shape the hydraulic gradient was changed that during 
part all the filling period water would flow from the lock chamber into 
the culvert through number the ports the upstream end the manifold. 
The prototype port area exceeds that the model 16% the river wall 
and 18.2% the land wall. Excess area the prototype equivalent 
that three ports. The effect the excess port area the prototype 
indicated Fig. 25(a). Discharge through ports and was reversed, flowing 
from the lock chamber into the culvert during the entire period the filling 
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operation, and, although the flow port was not measured, undoubtedly 
was negative also. Therefore, the distribution patterns the model and 
prototype manifolds cannot correlated. 


Port Number; River Wall 


AVERAGE 


Description 


per Second: 


Discharge refers flow into, and out of, lock chamber through the ports the river wall. Proto- 
type. Model. 
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very close similarity the flow distribution the model and the 
prototype during emptying operations shown Fig. 25(b). The distribu- 
tion flow through the ports was only slightly affected variations head, 
and the discrepancies port dimensions the prototype had less effect the 
emptying than the filling operation, undoubtedly because the flow the 
venturi-shaped ports was converging instead diverging. 

Velocity Discharge discharge sections one the locks, 
with the exception slight differences few dimensions, were constructed 
like the model. Flow conditions during the first part the emptying opera- 
tion differed the two structures, but after the valves were opened completely, 
distribution the flow from the ports was found quite similar. 


Port Number; Land Wall 


Fig. shown the distribution flow from the discharge ports the 
land wall the model for head ft. The curve for the prototype shows 
average distribution for four heads (25 ft, ft, ft, and ft), all which 
occurred after the valves were open and flow conditions were relatively stable. 

These curves for model and prototype are similar having the same 
general trend larger discharges the downstream ports. The increase 
the discharge for each successive port downstream not orderly the 
discharge section (which has culvert variable area) the lock- 
chamber section which has uniformly spaced ports and constant area culvert. 
The percentages culvert flow carried the eight upstream and eight down- 
stream ports the discharge section were 44.4% and 55.6% for both model and 
prototype, variations between model and prototye having balanced each 
half the section. 

Lock coefficients for lock hydraulic systems presented 
this paper were computed for the period after the valves were fully opened and 
were based the minimum culvert area. the prototype lock system the 
section having the least area through which the water passes, both filling 
and emptying, the 12-ft 12-ft culvert below the transition section. 
the model the section minimum area was the throat the lock-chamber 
ports, this area being 250 compared with the minimum culvert area 
288 ft. The culvert area the prototype was the same the model, but 
the cumulative port area the minimum section was 292.8 ft. 
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The so-called lock coefficient which bulk coefficient representing 
the efficiency the lock, equal the ratio the observed discharge 
the theoretical discharge computed the plain orifice equation using in- 
stantaneous heads: 


and cross-section area two culverts. The filling and emptying 
cients determined for the prototype and model locks are shown Fig. 
and Table 


TABLE 


No. Valve cham- Lock- 
ber Throat cham- 
end ber 
end 
(1) (2) (3) (4) (6) (7) (8) (9) (10) 
(a) 
Lock 
right 
2nd Lock 
Model 1.00 7.20 12.6 0.148 24,300 160 240 240 
1.07 7.08 12.4 0.273 24,300 160 240 240 
Model 1.00 6.15 12.5 0.117 24,300 128 102.8 211.8 
Prototype 1.12 5.50 12.7 0.118 24,300 128 102.8 280.8 
(b) OPERATION 
Lock 
Model 13.48 8.6 0.045 288 250 640 
5.12 13.30 7.3 0.030 288 292.8 640 
Lock 
Model 1.00 8.03 11.7 0.059 24,300 160 240 240 
1.07 7.37 11.8 0.106 24,300 160 240 240 
Model 1.00 7.60 9.3 0.029 24,300 128 102.8 211.8 
Prototype 1.12 7.60 8.8 0.052 24,300 128 102.8 280.8 


Maximum rate rise fall. Maximum water surface slope. Distribution based equal number 


The average coefficient for filling the prototype lock based the culvert 
area, which the smallest section the course flow, was 0.91 shown 
the legend. the model the coefficient based the culvert area was 0.87. 
However, the smallest area the course flow the model was the throat 
the ports, and when the model coefficient correlated the prototype 
value applying the ratio culvert port area, 1.15, increased 
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1.00, which not agreement with the prototype value. would appear 
that the venturi type flow the ports during filling makes comparison the 
lock coefficients unsatisfactory because the divergence effect the model 
and prototype was not the same. 

The average coefficients for the emptying operation, based the least 
culvert port area, are for the prototype 0.73 and for the model 0.64 1.15, 
0.74, which are agreement within 2%. The flow through the lock- 
chamber ports during emptying operation not similar that flaring 
tube, and the decrease the efficiency the model indicated the lock 
coefficient directly proportional the reduction the area the controlling 
section. 


AND PROTOTYPE 


Intake Ports Lock-Chamber Discharge Ports 


Scale 


Lock Transi- ratio 
hydrauli 
Up- Down- Up- Down- Up- Down- 
stream stream stream stream stream stream 
half half half half half half 
(11) (12) (13) (14) (15) (16) (17) (18) (19) 
OPERATION 
0.87 0.93 45.5 54.5 40.4 59.6 
0.91 0.92 47.0 53.0 17.3 82.7 
0.62 31.3 68.7 
0.88 0.93 45.5 54.5 
46.5 53.5 
0.95 0.95 46.0 54.0 40.8 59.2 
OPERATION 

0.64 42.1 57.9 44.4 55.6 
0.73 40.5 59.5 44.4 55.6 
0.56 30.2 69.8 
0.58 35.0 65.0 
0.63 44.3 55.7 44.4 55.6 
0.63 46.5 53.5 43.6 56.4 


ports upstream and downstream. Model test Norris Laboratory. 


the orifice equation, used computing lock-chamber discharge, 


the head, for filling, the upper pool elevation minus the lock-chamber 
elevation, and, for emptying, the lock-chamber elevation minus the lower pool 
elevation. This head simply the hydrostatic head and not the complete 
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effective head acting the column water the culverts. Considering 
the filling operation, the lock-chamber water surface does not stop when 
reaches upper pool elevation but rises above it. Therefore, the instant when 
the lock full the effective head not zero but has some positive value. 
When the dynamic head factor included the equation for discharge the 
resulting coefficients for decreasing heads become smaller instead larger. 
Thus, the prototype tests, shown Fig. 27, the ascending value the 
coefficient the head decreases can attributed the fact that the effective 
head was not used the computation, whereas the model tests the decrease 


> fea] @ 
Water Surface Elevation, Feet Above Mean Sea Level 


Rate Rise, Feet per Minute 


Filling Curve 
Rate of Rise oy 
Water Surface Slope 

Upper Pool (Ft) 553.9 

Lower Pool El (Ft) 503.9 
Lift 50.0 
Time of Operation, in 

Minutes and Seconds 7:05 
Time of Valve Opening 1:04 


Time Minutes 


CHARACTERISTICS 


the coefficients for low heads may indicate that the hydrostatic head more 

nearly equal the effective head the model than the prototype. 
Coefficients based the effective head (that is, static plus dynamic) and 

coefficients based static head alone have the same adverse characteristic 


560 
555 
4 
W 
LEGEND 
505 
tir 


October, 1942 


MODEL-PROTOTYPE CONFORMITY 


not being entirely constant throughout the lock operation. 


1315 


the other hand, 


the effective head not easily measured the static head. view the 
difficulties establishing the exact value the effective head and the ques- 
tionable advantage, the comparison model and prototype, the coefficient 
thus obtained over that based the static head, all lock coefficients presented 
this paper are based the static head alone. 


Lock Stage measured with respect time during filling 
operations the Wheeler model and prototype lock chambers are shown 
Fig. 28. The points which define the curves are averages two stage mea- 


surements, one near each end the lock chamber. 


the lock-chamber ports are shown 
Fig. 29. 

the prototype, was necessary 
make the tests with existing eleva- 
tions 553.9 for the upper pool and 
503.9 for the lower pool; tests had 
been made previously the model 
with elevations 556.0 and 506.0, 
respectively. each case the head 
was ft, but the pool levels used 
the model were 2.1 higher than 
the prototype. Fig. the filling 
curve the model shown ob- 
served and also adjusted agree with 
the lower pool elevation the proto- 
type. The prototype lock filled 
min sec, whereas the model filled 
min 12sec. Maximum rates rise 
were 12.4 and 12.6 per min, and 
the slopes the rate rise curves were 
1.90 and 1.94 per min?, respectively. 

The emptying curve the model 
was similarly observed and adjusted 
agree with the upper pool elevation 
the prototype. Emptying time 
the prototype was min sec, and 
the model was min sec. 


Typical sections through 


560.0 


Top Lock Wall 
: 
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Center 
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3' o" 


Maximum rates fall the water surface 


the chambers were 11.8 and 11.7 per min, respectively, and corresponding 
slopes the rate fall curves were 1.9 and 1.6 per min?. 

The valve opening periods were nearly identical the prototype and model, 
and the crests the water surface change curves occurred about the same 


time the two structures. 


the filling tests the model was slightly more 


efficient than the prototype, whereas the emptying tests the reverse was true. 
Velocity Lock-Chamber Ports.—Distribution flow the lock-chamber 


ports for the left river wall the prototype and for the right wall the 
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model shown Fig. 30(a) for the filling operation. The curves for the 
model were based velocity tests with steady flow through the system, whereas 
the curves for the prototype were based velocity measurements made during 
filling operations. The distribution flow observed the lock-chamber 
ports the model good agreement with that the prototype (see Table 
3A). Inflow into the lock through the upstream five ports the prototype 
was 32.4% the total, and the model 31.3%. The corresponding values 
for the downstream five ports were 67.6% and 68.7%, respectively. 


30) 


AVERAGE AVERAGE 


Description 


Cubic Feet per 
econda: 


Discharge refers flow into, and out of, lock chamber through ports river wall. Prototype. 


head. Model 


good agreement also noted the distributions discharge the lock- 
chamber ports emptying operations one model and prototype locks which 
are shown Fig. 

Lock and emptying coefficients for one the lock 
hydraulic system the prototype and model are shown Fig. 31. These 
coefficients are based observed discharges computed from the average rate- 
of-rise and rate-of-fall curves and theoretical discharge computed from stage 
records (see Eq. average filling coefficients for heads above are 
0.61 and 0.62 for the prototype and model, respectively. Corresponding 
emptying coefficients are 0.58 and 0.56. The lock coefficient varies the 
square root the slope the rate-of-rise rate-of-fall curves. For that part 
the filling operation for which coefficients were computed, the model had the 
steeper rate-of-rise curve and the larger filling coefficient, and the prototype had 
the steeper rate-of-fall curve and the larger emptying coefficient. 

one lock, the pattern the coefficients another indicated that 
there was less difference between the static and effective heads the model 
than the prototype. 


Intake the model tests one lock was found that the 
distribution flow through the intake ports was fairly uniform and that 
was affected only slightly variations the head. The intakes used 
one model were also used another model, and further tests were made 
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the latter model determine the distribution flow through intake ports. 
prototype tests flow distribution one lock was found that the dis- 
tribution flow the four ports one side the intake section varied from 
20.5% 28% and that the pattern flow distribution was similar that 
for the intake ports another lock. 

Transition prototype and model transition sections are cor- 
related using the dimensionless coefficient discharge because correspond- 
ing dimensions the culvert which pressure measurements were taken were 
not similar. Coefficients for the prototype and model transition section meter 
and 0.93, respectively—are good agreement. 

Lock Stage with respect time the lock chambers (see 
Fig. 21) during filling operations the model and prototype are shown 
Fig. time required fill the lock chamber with 39.2-ft head the 
prototype min sec) was sec, about 11%, less than the time required 
fill the lock chamber the model. The rate-of-rise curves indicate that 
the water surface rose faster the prototype than the model during the 
time the valves were opening. believed that this difference due largely 
personal factor introduced the manual operation the valves the 
model. After the valves are open the rate-of-rise curves are better agree- 
ment, having slopes 2.77 and 2.57 per min? for prototype and model, 
respectively. The maximum rates rise were 12.7 and 12.5 per min the 
prototype and model, respectively. 

Emptying curves are shown Fig. The same period time, min 
sec, was required empty the model and prototype lock chambers with head 
39.0 ft, and the emptying curves were practically identical. Maximum rates 
fall the prototype and model were 8.8 and 9.3 per min, respectively, 
and corresponding slopes the descending rate-of-fall curves were 1.22 and 
1.27 per min?. 

Lock-Chamber Ports.—The distribution discharge the 
ports observed filling operations the prototype, the Iowa City model, 
and model the lock tested the TVA Hydraulic Laboratory Norris, 
Tenn., are shown Fig. 34. 

Velocities measured the lock-chamber ports the models with steady 
flow are shown the distribution flow curves percentages the total 
flow. Velocities measured the lock-chamber ports the prototype lock 
during filling operations were plotted against time, and from these graphs the 
port velocity distribution for given static head was obtained, using the time 
interval which the given static head occurred. 

The distribution flow during the filling operation was better the models 
than the prototype. Percentages total flow carried the seven ports 
the upstream half the chamber and the seven ports the downstream 
half the chamber are: For one model (average four heads), 47% and 
53%; for another model (25-ft head), 46.5% and 53.5%; and for the proto- 

type (average six heads), 40.8% and 59.2%. Patterns flow distribu- 
tion the two models are very similar, and they agree more nearly with that 


the prototype 35-ft head than with the average prototype distribution 
after the valves are open. 
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Water Surface Slope ---- 
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The distribution flow the lock-chamber ports during the emptying 
operation the model and prototype was better agreement than the dis- 
tribution during filling. Percentages total flow carried the seven ports 


Port Numbers 


OnE AVERAGE 
Description 


Cubic Feet 
per Second: 


Discharge refers the flow into the lock chamber through the ports the land wall. The curve 
for the 35-ft head plotted separate line Fig. 34(a), the remaining six this group being averaged. 
Fig. all seven are averaged. Prototype. Model. 
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the upstream half, and the downstream half the chamber are, for one 
model (25-ft head), 44.3% and 55.7%, and for the prototype, 46.5% and 
53.5%. The models had the more even distribution during the filling opera- 
tion, and the prototype during the emptying operation. 

Discharge preliminary tests were made the distribution 
flow the discharge ports one the models. The discharge sections 
one model were used also another model, and since the preliminary 
tests agree with similar data from one the models, further tests were 
considered necessary. 

From prototype port velocity tests one lock, the distribution flow 
the sixteen discharge ports the land wall were computed for seven heads 
ranging from ft. Distribution flow among the sixteen discharge 
ports remained approximately the same for the seven different heads, and 
the culvert flow carried the eight upstream and eight downstream ports 
was 43.6% and 56.4%, respectively, based the average distribution for 
the seven heads. 

Lock coefficients for the lock hydraulic systems the 
prototype and model, for heads greater than ft, are: 


Observations on: Filling Emptying 
0.95 


For the emptying operation the coefficients are good agreement, having the 
same average value for prototype and model, but for the filling operation 
the prototype coefficients average higher than those for the model. 

The prototype ports and the model ports for which data are presented 
this paper have the same dimensions the throat, but the prototype ports 
are larger the lock-chamber end. Consequently, the divergence angle for 
the prototype, 11° 44’, larger than that for the model port, 50’, and the 
ports are not geometrically similar. This greater divergence angle undoubt- 
edly factor contributing toward the higher coefficients and faster rate 
filling observed the prototype. 

Although the area the culverts was used computing coefficients for 
the lock hydraulic system, the minimum area the prototype and the model 
the cumulative area the throats the lock-chamber ports. If, instead 
the culvert area, the minimum port area was used, the average coefficient 
for emptying would 0.78. For emptying one lock the average coefficients 
based least port culvert area were 0.73 and 0.74 for prototype and model, 
respectively. 

The coefficients for prototype lock based static head gradually in- 
crease the head decreases, and for the model they gradually decrease the 
head decreases, until just before the end the lock operation where values 
the coefficient increase sharply. 
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SUMMARY AND CONCLUSION 


Table the results model and prototype studies three locks are 
summarized. includes operation data and areas critical sections the 
lock hydraulic systems. 

Tests one lock offered the best opportunity for comparing model and 
prototype performance because the corresponding structures were geomet- 
rically similar and lock operation data were agreement. Moreover, 
will noted that this lock had the most consistent agreement between 
model and prototype the three locks tested, the only material difference 
the performance the corresponding structures being the longitudinal 
slope the water surface for which readings the prototype lock were larger 
than the model. 

For model and prototype locks another there was agreement lock 
operation data, and the two structures were geometrically similar except 
the lock end the lock-chamber ports. more efficient filling operation 
was observed the prototype than the model, the former having the greater 
divergence angle the lock-chamber ports. Furthermore, during the filling 
operation, the less efficient model had the more even distribution flow 
the lock-chamber ports. interest note the performance Gunters- 
ville model and prototype locks during the emptying operation. There was 
agreement operation time, lock coefficient, and the distribution flow the 
lock-chamber ports irrespective the difference lock-chamber port diver- 
gence angles. The throat the lock-chamber port, for which the cumulative 
area was the minimum the system, appeared effective the control 
section the system during the emptying operation. 

There was less agreement the performance one model and proto- 
type locks because, addition dissimilarity the throat area the lock- 
chamber ports, lift and valve time the two structures were not agreement. 
However, the performance this model and prototype during the emptying 
operation, similar that other locks, showed agreement the lock co- 
efficients based the minimum areas the respective systems and the 
distribution flow the lock-chamber ports. 

When geometric similarity does not exist the important elements two 
structures not expected that their operations should conform 
the basic laws similitude. However, where geometric similarity did ob- 
tain, and where analogous functions could correlated means the basic 
laws, this study substantiates the theory upon which the use hydraulic 
models founded and, believed, tends justify the confidence engineers 
have placed upon this method analyzing difficult problems design and 
operation within the past few years. 
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DAM OUTLET WORKS WESTERN STATE 


The principal hydraulic model tests the dam outlet works (34), which 
included flow through the intake tower, penstock turbine, and outlet-valve 
manifolds, were made scale), extremely smooth, sheetmetal 
model which had friction coefficient 0.021 Reynolds number 
80,000 and used maximum per sec water. Performance and 
efficiency tests the 30-ft diameter prototype, utilizing maximum discharge 
20,000 per sec with measured velocities per sec and 
Reynolds number from made possible obtain data from 
which various hydraulic features the model and prototype could com- 
pared. These included loss coefficients for the trashrack, bends, branches, 
valves, and straight-pipe friction. full the data was presented 
before the hydraulic section the Annual Convention the Society 
July, 1940 This paper summary such the results may 
generally useful hydraulic design, together with suggestions concerning the 
technique field measurements which, found practicable, may make possible 
better comparisons between models and prototypes future tests. 


The letters Figs. and indicate points which measurements were 
taken. Distances are given units pipe diameters, that the 25-ft section 
being converted into its 30-ft hydraulic equivalent. Three series compre- 
hensive tests were run from 1937 1939 with heads varying from 213 337 ft. 
During the tests the turbines were shut down, varying discharges being obtained 
manipulating the 84-in. outlet valves. Pressures were measured with 
calibrated Bourdon gages and differential mercury manometers, utilizing the 
dead-ended penstocks for piezometric connections. Tests showed that such 
large connections taking off angle resulted the gages registering 
excess pressure 0.04 

Discharges were obtained impact-pressure tube traverses right angles 
one another 96-in. section pipe located immediately downstream from 
reducer and front valve point Utilizing these measurements 
and the simultaneously measured heads the other valves, the total discharge 


for the entire battery valves was computed. The coefficient the pitot 
tube the formula, 


Civ. Engr., Bureau Reclamation, Denver, Colo. 
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obtained from still-water rating, was corrected for wall proximity effects 
and for the effects due curved streamlines caused the valve and the 
obstruction the streamlined struts used support the instrument (see 
Fig. 37). These latter corrections were determined the photoelastic labora- 


Water Surface 
1157 


050 


60° 
TURBINES 


tory the Bureau Reclamation, using the electric-analogy method (see 
Fig. 38). Fig. shows the results found the laboratory, the close corre- 
spondence between the laboratory result and that obtained the rating 
channel for obstruction being interest. Since offered 


Lake Water 
Surface 


Lost Head 


Distance from (Equivalent 30-Ft Diameter) 


the struts variable percentage each circular lamina the pipe, 
separate coefficients were applied each individual reading the pitot 
traverse. 

Plots the velocity and pressure distribution the cross section indicated 
that upstream from the pitot tube the friction 22.5 diameters pipe plus 
15% reducer did not introduce enough resistance normalize the upset 
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velocity distribution due branching, the velocity contours showing the 
crescent-like form (typical flow below bend) superimposed the upstream 
cone pressure caused the downstream valve. The latter caused differ- 


Coefficient 


Percentage Obstruction 


ences.in pressure water, nearly 0.5 between the pipe wall 
and the center the pipe. Worth considering valve design the fact that 
the energy flow was also unsymmetrical, the half the pipe away from the 
center curvature the branch carry- 
ing 10% 15% more energy than 
the inside, with the possibility that pres- 
sure unbalance existed within the valve. 

the whole, the pitot-tube gag- 
ings were reasonably self-consistent. 
estimate their errors and com- 
parison with current-meter gagings 
the river indicated standard error 
for the largest discharges the -order 
2%. 

The interpretation ‘of losses pen- 
stock, with fittings but relatively few 
diameters apart, depends upon abil- 
ity separate fitting losses from those 
due straight pipe. pipe entrance 
fitting tends result contracted 
jet high-velocity water which changes 
velocity distribution, the 
length required being dependent the 
roughness the pipe wall and Reynolds’ number knowledge the 
energy grade based only the readings wall piezometers, the profile 
the energy grade characteristically will found show apparent increase 
energy the first diameters following fitting, pass through point 
inflection and then will become concave upward, approaching asymptotically 
the friction gradient. The model test, moreover, indicated the possibility 
that nodes pressure about 0.03 located diameters apart possibly 


Center Line 


38.—Tray 


0.92 
Field Observation 
0.96 
From Rating Channel 
~ 
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caused the expansion and contraction the high-velocity jet, were super- 
imposed the profile immediately downstream from the branch fitting. 
This possibility requires more investigation before can taken fact. 


ution 
wo 


Horizontal Velocity 


1.02 


Distance from Center Line 


0.5 W=25D 


The foregoing phenomena make almost impossible power installation 
find run pipe which can certain that only straight-pipe friction 
reflected differences the readings wall piezometers. For the present 
experiments, the most reliable result estimates straight-pipe friction were 
obtained utilizing the measurements made the 13-ft turbine penstocks 
connection with the Gibson (36) efficiency tests the turbines. The 
measuring section Fig. included diameters pipe, located diam- 
eters from the branch. The results, Fig. 40, gave value 
0.0207. Utilizing the model tests, which showed that such location 
about 0.03 the measured drop should ascribed the upstream diversion 
loss, the friction coefficient for 13-ft diameter straight pipe was estimated 
and Nikuradse’s sand roughness factor 0.00134. This latter factor 
can interpreted meaning that the bitumastic-covered surface the 13-ft 
pipe was rough hydraulically had been coated with uniform sand 
grains 0.1 in. diameter. According the field observers, broom marks 
and joints were readily visible the bitumen, 20-ft straightedge showed 
openings below it, and flattened bolt heads projected in. from the welded 
interior the pipe. Friction coefficients for other sizes pipe were computed 
from these measurements, assuming them vary inversely the value 
for the 30-ft pipe being 0.0134. These unexpectedly high values can 
compared with the results obtained from concrete pipe comparable values 
plotted the graph Fig. 40. Also the same graph are plotted over- 
all losses per diameter for each the sections measured—A etc. 
(see Fig. 35). Since the values are approximately constant 0.03 per 
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diameter, this figure can used over-all check for hydraulic designs 
which the hydraulic losses due friction and fittings are separately estimated. 


The general results the test 


are best shown the form graphs. 


Fig. shows, terms velocity head, profiles the accumulated total 
energy loss and the total hydraulic loss ascribable fittings for both the 


model and prototype. 


Fig. 40(a) comparison the drop the intake 


tower and the valve discharge coefficients plotted against For the former 
smooth curve has been passed through the 1939 and 1938 measurements, based 


Junction Loss, 


ischarge 


Fie. Loss 


velocities determined the pitot tube. 


gagings are shown for comparison. 
manifold plotted against the perce 


The results from the river 
Fig. shows branch losses the valve 
ntage total flow diverted the branch: 


Curve Reynolds’ number 
(Fig. 41) Structure (in millions) 
Bureau Reclamation model............... 0.05 
Thoma model (interpolated 0.10 
Bureau Reclamation prototype............ 40.00 


Tr 
1.0 
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Table compares the various loss coefficients for the entire structure orig- 
inally estimated, found from the model, and estimated from the prototype 
tests. 

Quantitative correspondence between model and prototype can expected 
only where Reynolds’ criteria can satisfied and where absolute roughness 
the model can scaled down accordance with model laws. This was not 
possible the Boulder model which, while smooth possible, was materially 
rougher than demanded similitude. Referring the fitting loss Fig. 36, 
can seen that the prototype point diameters from the entrance, 
there was indicated loss about one half that measured the model and 
point diameters downstream, value 85% large. will noted, 
moreover, that point the accumulated fitting loss the prototype 
shown only one half that observed diameters upstream immediately 
downstream from the intake. will realized once that this decrease 
accumulated loss downstream measured wall piezometers due the 
fact that the profiles are only nominal energy grades and include corrections 
for actual velocity distributions. appears probable from these tests that 
the prototype, with relatively less friction than the model, requires greater 
length, measured diameters, regain its pressure head following the con- 
traction below the intake. The required length may analogous the 
distance required produce normal distribution below intake smooth 
pipe, for which Goldstein (37) gives the formula 


Applied the present tests, four times the length required the model, 
measured diameters, would required produce normal velocity distribu- 
tion the prototype. Whatever the reason, believed that the over-all 
measured prototype fitting loss, 0.85 (which includes entrance, two bends, 


TABLE 4.—Comparison Loss 


Loss CoEFFICIENT 
30-Fr PIPE 


No. Description 


Prelimi- 
Proto- 
nary Model 
design type 


Intake drop from the lake within the 0.81 
Estimated entrance loss diameters from intake................. 1.10 0.68 0.63 
Estimated vertical and horizontal bend 0.40 0.08 
Estimated losses for eight 12-in. tie rods and four 13-ft branch 
Measured total entrance and fitting losses diameters.......... 1.75 1.00 0.85 
Branch loss with 100% diversion, measured diameters from branch 0.26 0.14 


the disturbances caused openings for four branches, together with their 
eight tie rods, and reducer), provides satisfactory check the 
model results. Table this over-all loss has been broken down into its 
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components and compared with the standard loss coefficient originally used for 
the preliminary design the outlet works and with those later obtained from 
the model. Although the breakdowns are only estimates based con- 
sideration all the data, these fitting coefficients should helpful designers 
until better values become available. estimated that the values have 
standard error the order 15%. 

The entrance loss into the tower through its radial ports, 0.63 given 
from the foregoing breakdown, sometimes confused with the piezometric 
drop which occurs from the lake the water surface within the tower. The 
latter (as shown the bulletin (34) the model tests) involves reac- 
tion head due the water turning downward within the tower and not 
absolute measure entrance losses. This coefficient tower drop, how- 
ever, obtained from various gagings, forms good means comparing 
performance model and prototype. The results are shown Fig. 40. 
Although first sight the variations the different tests may seem large, 
the 12% range the experimental values believed comparable the pre- 
cision the majority prototype tests. The increase the coefficient with 
increasing indicated the curve averaging the 1939 and 1938 prototype 
results, may, perhaps, partly accounted for the formation vortices 
the intake which were observed the 1938 measurements. that time the 
head the intake was only from four six times the mean velocity head. 
Since the velocity heads for some the stream filaments passing around certain 
unavoidable sharp edges the intake, measured terms the mean velocity, 
may have been even larger multiple the mean head with the corresponding 
absolute pressures approaching zero, this condition which should taken 
into account designers, the vacuum may cause cavitation and the en- 
training vortices large quantities air. Without complete exploration 

The model showed loss through the screens the prototype 
the loss was less than one tenth this value. This undoubtedly was due 
the fact that velocity through the model screens was below critical 
value. model should not reproduce from the prototype details within which 
subcritical velocities may exist during the model tests. 

considerable interest designers the loss due closely connected 
90° vertical and 40° horizontal bend. Although the computed loss coefficient 
for the prototype 0.08 may have considerable percentage error, the total 
measured loss for all the fittings was such that the bend loss, regardless 
assumption, must have been small. the writer’s knowledge this one 
the few prototype tests that has been made bends. 

The junction loss for the branch from the manifold shown Fig. 
Curve for the model interpolated from very complete series tests 
described the bulletin (34). Curve was likewise interpolated from 
Thoma’s well-known tests (38). Curve was obtained from the results 
commercial laboratory connection with the design the turbine 
manifold. will noted that the observations defining its averaging curve 
are widely scattered. Any one who has examined critically the technique 
which fitting losses must obtained, set forth the bulletin (34), will 
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appreciate how difficult obtain reliable results and will know that 
only the exceptional test can the absolute values the results accepted. 
enough state that coefficients fitting loss are not fixed quantities 
depending merely They tend decrease with increasing smoothness 
and with increase Reynolds’ number, and greater the greater the dis- 
tance below the fitting which they are measured. Also, noted 
that with multiple similar branches from manifold, the loss will different 
each successive offtake due change the value the percentage the 
main flow diverted. 

The needle valve discharge coefficients are shown Fig. 40. will 
noted that the model with increasing the coefficient increased, whereas 
the prototype decreased. The first case undoubtedly reflects the decrease 
friction through the valve with increase the second case analysis 
the velocity distribution due curved streamline within the valve indi- 
cated the probability that certain points vacuum was approached high 
values Once this point was reached, the effective area the waterway 
would reduced, cavitation would commence, and the discharge coefficient 
would decrease. The performance the prototype was line with these 
predictions. 

believed that the whole the use model for predicting quanti- 
tative results the outlet works was pre-eminently successful, nearly halv- 
ing the error design that would have existed without its use. the 
other hand, much more experimental data and development theory are 
needed predictions losses within 10% are desired the prototype. 
ticularly required are methods for forecasting pipe friction terms quanti- 
tative pipe roughness and for estimating fitting losses various distances 
downstream the stretch before velocity distribution has become normal. 
Finally, must not forgotten that there way which factor 
safety can applied the economic layout the flow line hydroelectric 
development. overestimate hydraulic losses leads the overdimension- 
ing the conduit with excess interest charges perpetuity, whereas under- 
estimate leads continuous waste power excess friction losses. 
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SPILLWAY COEFFICIENTS 
EASTERN DAM 


The determination discharge coefficients for the dam spillway crest 
gates described this paper. Actual coefficients calculated for the gates 
from records their performance during the flood January and February, 
1937, are compared with the coefficients predicted from tests made 
scale model the Bureau Reclamation. This comparison indicates 
that the coefficients predicted the model average 3.9% larger than those 
actually measured the dam during this period. 


GENERAL CONDITIONS 


The dam was built part program navigation, flood control, and 
power development. controls the discharge 2,912 miles land 
drained two rivers. 

concrete gravity structure having total length 1,860 and 
maximum height, the roadway the top El. 1061, 265 The spill- 
way crest El. 1020 and 300 long, divided into three bays each 100 
long. Three hydraulically operated steel drum gates are installed the con- 
crete spillway crest, and the crests the gates may raised any elevation 
between 1020 and 1034 retain flood peaks. Eight outlet conduits permit 
the discharge water when the reservoir level below the spillway crest. 
The maximum discharge these conduits with reservoir level El. 1020 
approximately 37,000 per sec. The power house contains two generating 
units, each having capacity 56,000 kva. 

Fig. shows general plan and elevation the dam and section through 
the spillway. Fig. gives the dimensions the spillway crest detail. 

Flood January, 1937.—Unusually heavy rains resulted flood stages. 
During this period the entire flow rivers above the dam was stored lake. 
This necessitated raising the spillway drum gates hold water maximum 
about above the concrete spillway crest. Late January was seen 
that moderate discharges the dam would not cause further rises the 
river, and accordingly water was released from the lake provide storage 
space for additional runoff. Much the water was released lowering the 
crest gates and allowing discharge over the spillway. Advantage was 
taken this opportunity obtain the necessary measurements for determin- 
ing the discharge coefficients the spillway for heads the gates great 
9.8 and for discharges great 30,000 per sec. 

(Deleted censors.) 
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Data AVAILABLE 


Gate Operating program water release the dam re- 
quired the discharge held approximately constant for considerable 


S 
~ 


20-Ft Diam. 


Tailrace 


Switchyard 
Spillway 


DOWNSTREAM ELEVATION 


Center Line Road 1061.0 


Crest 1020.0 


Center Line 
865.0 


SECTION THROUGH SPILLWAY AND SLUICE 


time. was necessary change the elevation the crest gates, 
since the water level the reservoir was changing. Accordingly, the gates 
were reset about every two hours. complete record the operation the 
spillway crest gates, giving the time which changes were made the gate 
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settings and the elevations which the gates were set, was available from the 
records kept the power house. Records reservoir water-surface elevations 
were also available. 


Spillway Crest 


Discharge gaging station has been established the 
Geological Survey point 0.9 mile below the dam. standard recorder 
was installed, and continuous record for the period spillway discharge 
was available. The rating curve well defined since the control solid 
rock diagonal reef about 150 below the gage. number discharge 
measurements were made during the flood period. Only one these 
much from mean rating curve, and most the measurements are less 
than from the curve. 


CALCULATION DISCHARGE COEFFICIENTS 


Notation.—The letter symbols used this paper are defined where they 
first appear and are assembled for convenience reference under 
and the end the Symposium. 

Correlating the computing the discharge coefficients, was 
necessary that the discharge over the spillway should the same that the 
gaging station and that the river stage the gaging station should steady. 
previously mentioned, the gates were adjusted about every two hours, 
causing fluctuations the discharge. most cases 2-hr period was sufficient 
allow the stage the gage become steady, since the changes discharge 
were usually small. When observations were used which were only two hours 
apart, the elevations the gate crest and the reservoir level were taken 
immediately before the change—that is, the end the 2-hr period. The 
gage height the gaging station was also taken the end the period. 
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When large changes discharge occurred, greater time was required for the 
river stage become steady. When appeared from the recorder chart that 
steady condition had not been established the gage, the data were not used. 
data were used determine discharge coefficients there was any doubt 
their value, either because unsteady flow for any other reason. Use 
river discharge the gaging station represent the discharge over the spillway 
the same time justified because the care rejecting data which did not 
represent steady flow and because the relatively small storage between the 
dam and the gaging station. 

Table Cols. inclusive, lists the data regarding reservoir and 
gate-crest elevations obtained from the power-house records. Col. Table 
gives the gage height, feet, the gaging station, corrected from the con- 
tinuous chart. 

Calculation rating curve for the gaging station was con- 
structed from the discharge measurements plotting the discharges against 
gage heights and drawing curve through them give greatest weight 
those measurements that were made during the period observation. 

There practically inflow the river between the dam and the gaging 
station. The gaging station records indicate that, just before spillway dis- 
charge began, the combined inflow and leakage through the dam was about 
per sec. Col. 10, Table gives discharges taken from the rating curve, 
less per sec for all discharges less than 10,000 per sec. For dis- 
charges greater than 10,000 per sec, the correction was neglected since 
amounted less than 0.5%. The discharges Col. 10, determined de- 
scribed, are considered the spillway discharges for the corresponding 
reservoir and gate-crest elevations Cols. 

Calculation Discharge Coefficient discharge coefficient, was 
computed means the customary weir equation, 


which: discharge, cubic feet per second; discharge coefficient; 
length weir crest, feet (for one dam, 300 ft); and head weir 
crest, feet. 

Cols. and Table show that the crest gates were not all the 
same elevation but varied somewhat, giving different heads the gates. The 
discharge coefficient varies slightly with both the elevation the gate crest 
and the head, but the variation not perceptible for the differences 
The maximum variation head between gates 5.7% for Observation No. 
The results the first six observations are quite widely scattered, and, these 
are omitted from consideration, the maximum variation head 1.5% 
Observation Variation the discharge coefficient due this change 
head may neglected. slight error also introduced the average head 
used the weir formula, but with maximum head difference 1.5%, 
this error also may safely neglected. Col. Table gives the average 
head the three gates for each observation. Values calculated from 
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TABLE 
(Calculation Prototype Coefficients 


Elevation? Crest 
Observa- 
tion No. Time 
Gate Gate Gate 


February 
8:50 a.m. 30.49 24.44 24.44 24.52 24.47 
5:00 p.m. 30.37 24.34 24.34 24.42 24.37 
10:15 p.m. 30.32 24.24 24.24 24.32 24.27 
February 


8:45 p.m. 
February 
2:40 a.m. 


p.m. 


February 
3:15 a.m. 29.03 21.32 21.32 21.40 21.35 
12:00 28.77 21.12 21.12 21.20 21.15 
2:00 p.m. 28.68 21.02 21.02 21.10 21.05 
4:30 p.m. 28.61 20.00 20.00 20.00 20.00 
February 
12:05 a.m. 28.37 20.71 20.71 


Add 1,000 Cols. inclusive, obtain actual elevations. 


(1) (2) (3) (4) (6) 
January 
11:05 a.m. 30.35 29.28 29.28 29.34 29.30 
1:10 p.m. 30.41 28.68 28.68 28.74 28.70 
3:10 p.m. 30.45 28.27 28.27 28.33 28.29 
5:00 p.m. 30.47 27.86 27.86 27.93 27.88 
February 
8:30 a.m. 30.70 27.56 27.56 27.63 27.58 
; 6 11:05 a.m. 30.73 27.05 27.05 27.12 27.07 
7:00 p.m. 30.78 26.10 26.10 26.17 26.12 
9:00 p.m. 30.78 25.74 25.74 25.82 25.77 
11:00 p.m. 30.78 25.54 25.54 25.62 25.57 
February 
1:00 a.m. 30.77 25.20 25.20 25.27 25.22 
3:00 a.m. 30.77 24.90 24.90 24.97 24.92 
5:00 a.m. 30.75 24.55 24.55 24.62 24.57 
17 8:25 a.m. 30.17 24.04 24.04 24.12 24.0 
18 12:45 p.m. 30.08 23.74 23.74 23.82 23.77 
eS 19 1:45 p.m. 30.03 23.54 23.54 23.62 23.57 
2:45 p.m. 30.01 23.13 23.13 23.21 23.16 
29.98 22.53 22.53 22.61 22.56 
29.86 22.33 22.33 22.41 22.36 
29.73 22.12 22.12 22.21 22.15 
10:00 a.m. 29.54 21.92 21.92 22.00 21.95 
1:00 p.m. 29.45 21.82 21.82 21.90 21.85 
26 00 29.49 21.72 21.72 21.80 21.75 
27 - 29.13 21.52 21.52 21.60 21.55 1 
f 33 1:30 a.m. 28.33 20.61 20.61 20.69 20.64 : 
4:30 a.m. 28.23 20.51 20.51 20.60 20.54 
8:00 a.m. 28.16 20.41 20.41 20.50 20.44 
10:00 a.m. 28.17 20.00 20.00 20.00 
2:05 p.m. 28.10 20.11 20.11 
4:10 p.m. 28.07 20.11 20.11 
7:25 p.m. 28.05 20.31 20.31 20.40 20.34 
10:10 p.m. 27.95 20.21 20.21 20.30 20.24 
February 
6:00 a.m. 31.05 22.63 22.63 22.71 22.66 
9:00 a.m. 31.09 22.28 22.28 22.36 22.31 
12:00 31.10 21.72 21.72 21.80 21.75 
12:00 mid. 30.96 21.12 21.12 21.20 21.15 
February 
4:30 a.m. 30.86 21.02 21.02 21.10 21.05 
8:00 a.m. 30.74 20.91 20.91 21.00 20.94 
10:00 a.m. 30.67 20.81 20.81 20.89 20.84 
2:00 p.m. 30.58 20.71 20.71 20.79 20.74 
7:00 p.m. 30.43 20.61 20.61 20.69 20.64 
9:30 p.m. 30.35 20.51 20.51 20.60 20.54 
February 
2:30 a.m. 30.16 20.31 20.31 20.40 20.34 
5:00 a.m. 30.07 20.21 20.21 20.30 20.24 
10:50 a.m. 29.82 20.11 20.11 20.20 20.14 
1:45 p.m. 29.72 20.00 20.00 20.00 20.00 
2:15 p.m. 29.70 20.00 20.00 20.00 20.00 
3:15 p.m. 29.65 20.00 20.00 20.00 20.00 
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189.0 


Total, Nos. and 57, 189.0; average, 


ee: 


Observa- 
3.56 3.82 1.506 1 y 
4.17 3.77 1.510 2 
3.96 3.74 1.512 
4.16 3.70 1.516 
3.37 1,519 5 
3.36 19823 6 
3.40 133 7 
3.43 138 8 
3.39 9 
3.36 
3.41 149 1l 
3.31 
3.18 13 
71 18 
75 19 
20 
91 21 
93 22 
196 23 
98 24 
99 25 
26 
108 30 
513 35 
42 
43 
500 
508 47 
q 510 48 ; 
512 50 
514 52 
15 53 
617 54 
17 55 
617 56 
q 17 57 
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the formula, 


are given Col. 11. 


CoMPARISON WITH 


Determination Model Discharge coefficients for 
the dam spillway were measured the Bureau Reclamation 
scale model (39). was found that because the shape the drumgate 
crest, which changed with the elevation the gate, Eq. was neither 
constant nor function the head. The Bureau found that the tests made 
the model could represented the equation, 


which and are constants which vary only with the elevation the gate 
crest. Values and for each elevation the drum-gate crest are given 
(39a) Fig. 44. 


1034 


Elevation Gate Crest, Feet 


3.8 
Coefficient Discharge, 
1.46 1.48 1.50 1.52 1.54 1.56 1.58 1.60 1.62 


3.6 


compare the model tests with the results observations the proto- 
type, was desirable that coefficients discharge expressed the same 
terms for each. The data obtained the prototype not cover sufficient 
range permit the determination and However, values can 
computed from the Bureau’s tests. the two expressions for Eqs. 
and 10, 


Percentage Error: Model Coefficiente 


Papers 
b 1 1.5 ( ) 7 
1030 
1026 
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From this, 


Values and from Fig. are given Cols. and Table and the 
has been headed C,, order distinguish the coefficients computed the 
basis the model tests from those determined the prototype. Col. 15, 
Table compares the two values showing the percentage which 

few these values should not used for drawing general conclusions. 
The first six observations are for relatively small discharges and low heads 
the crest for which errors observation cause large percentages variation. 
Observations Nos. give discharge coefficients for single gate. They 
are probably correct but should not compared with the coefficients deter- 
mined for three gates, since the discharge over single gate subject much 
greater end contractions. these nine observations are eliminated, the aver- 
age the forty-eight remaining values +3.9%, indicating that the model 
tests predicted discharges which average 3.9% greater than those actually 
obtained. The reason for this not known, although possible that the 
model coefficients are not accurately known this range. generally 
believed that, other things being equal, small-scale models will give discharge 
coefficients somewhat less than the prototype. There reason discard 
any prototype observations other than those referred to, all questionable 
observations were omitted before calculations were made. believed that 
the values determined for the prototype are reasonably accurate. 


Percentage Error; Model Coefficients 


Discharge, Thousands per Sec 


Fig. shows the percentage difference the forty-eight observations con- 
sidered, plotted against the discharge. There definite tendency for the 
difference decrease with increasing discharge. this connection, 
interesting the tests performed the Bureau Reclamation, 
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only five actually fell within the range operating conditions represented 
the observations the prototype. other words, the model discharge 
coefficients for this range have been determined largely extrapolation. 
the region more model data approached, the difference becomes less. 
possible that the difference may decrease still further for higher discharges. 
Observations for discharges the range 60,000 100,000 per see will 
necessary before more definite comparison the model and its prototype 
can made. 


The comparison shows that model tests may relied predict proto- 
type discharge coefficients with accuracy least equal that stream- 
flow measurements. The average discrepancy 3.9% not serious and 
certainly less than might expected without the benefit model tests. 
quite possible that this difference might have been much less the measure- 
ments the prototype had been made range discharges that were more 
adequately covered model tests. This seems indicated the trend 
Fig. 45. 

Accurate quantitative comparisons model and prototype discharge 
coefficients are scarce because the model tests are generally made the region 
maximum discharge which seldom, ever, occurs the prototype. Many 
more comparisons could made the model were tested not only for the 
maximum discharge but also for the smaller discharges that occur more often. 
suggested that future tests spillway coefficients made include the 
range frequent prototype operation that will not necessary wait 
for the maximum flood obtain verification the model studies. 
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PRESSURE HEADS WESTERN DAM 


SYNOPSIS 


The main hydraulic features pertaining the design western dam 
were studied models. Among these features were the hydraulic character- 
istics the 50-ft 50-ft crest gates with water flowing beneath them (40), 
and the hydraulic effect departures the profile the constructed dam from 
the design profile. During the latter studies the pressures the model crest 
were determined for various gate openings. 

series piezometer taps, corresponding those the model, were placed. 
means these taps the pressure the crest the dam for various gate 
openings was observed. addition, the surface profile the prototype nappe 
was determined. 

This paper presents the results the pressure-head profiles the under- 
face the overflowing nappe both model and prototype. When the gate 
was the downstream position, the pressure profiles model and prototype 
checked fairly well; when the gate was the upstream position, however, 
anomalies developed that may serve the subject for some interesting 
speculations. 


replica 10-ft length the crest and gate was made scale 
piers there were two gate grooves; one the gate seal was upstream from 
the center line the crest and the other was downstream. 
simulate these positions the entire model crest block could shifted either 
position. The gate was held stay rods that moved parallel itself, 
and the rods were long that this movement virtually simulated the prototype. 
The effect the gate grooves was absent the model, watertightness being 
secured stanching strips rubber. 


The model tests were made primarily determine the effect slight 
departures the constructed crest from the design crest, and establish the 
tolerances that might permitted the crest and ogee profile. These experi- 
ments consisted essentially observing pressures the contact surface between 


(Deleted censors.) 
(Deleted censors.) 


| 
| 
, 
| | 
y ane 
. 
| 
7 
a 
7 
q 
q 
7 


and the three types distortions used the experiments. 


Platform Scales 


Recirculating 
Flume 


Crest Showing Gate 
Upstream Position 


A 


> 


d 


Crest Showing Gate 
Downstream Position 


Mud Sills 


addition tests the model, pressure the crest was observed 
for large range gate openings the middle one three bays 
model. The gate was the downstream position only, and the piezometer 
taps were 5-ft intervals downstream from the center line the crest. Since 


the prototype gate openings were and less, comparisons that value 
only could made. 


Prototype tests were made Bay No. the completed dam October, 
1939. The crest profile this bay was without appreciable departure from the 
design crest. The piezometer taps were spaced conform those the 
model except that the prototype, Piezometer No. upstream from 
the axis, was substituted for model Piezometer No. 9.58 upstream from 
the axis. The protoptype also had, addition those the model, Piezom- 
eters Nos. and located 27.5 and 42.5 downstream from the axis. 

The piezometer taps were made steel plate in. wide set flush with 
the concrete and bent the shape the design crest. From each tap 
iron pipe led through stopcock manifold inspection gallery the 
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the concrete crest and the overflowing nappe. Three distorted crests were 
studied for gate openings 0.25, 0.5, 1.0, 2.0, and 3.0 with the gate both 
the upstream and the downstream positions. Fig. shows the design crest 
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dam. Bourdon pressure gage, reading positive pressures above 
atmospheric, was connected the manifold El. —28.89 ft, that the 
pressure each piezometer could read turn. 


(a) PERFECT CREST PROFILE 
TESTS AND G-10 


Elevation, Feet Above Mean Sea Level 


Distance, Feet 


1.0" DISTORTED CREST PROFILE 
TESTS G-4 AND 


0123 


Scale Feet 


(c) 1.5" DISTORTED CREST PROFILE 
TESTS AND G-9 


(d) 1.0" AND DISTORTED CREST PROFILE 
TESTS AND G-11 


Two series tests were made, one with the gate the downstream groove 
—its normal position—and the other the upstream groove—an emergency 
position. addition observing the pressures the underface the nappe, 
the surface profile the nappe was obtained for distance 24.5 down- 


Axis Dam 
5 
24.00 
—1.0" = 10 . 9 
24.00 
+1.0" 
4q 
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stream from the gate seal when its downstream position and 
downstream from the gate seal when its upstream position. 
The surface profile for 

TABLE ELEVATION length near the crest 


(PROTOTYPE) was also obtained for 

when the crest 
Opening, Feet: was under head 47.8 ft. 
0.25 0.5 1.0 2.0 3.0 the flow through 


any one bay, but tests 
the model indicate that 
under that head this 50-ft 


23.2 23.4 23.9 24.9 25.9 26.9 


—9.0 
The surface profiles (see 
+ 8.0 -.e- | 23.6 |] 23.8 | 24.3 | 25.3 | 26.3 | 27.3 ; Se or 
+32.6 gantry crane which could 


moved perpendicular 
the axis the dam. There 
much spray from the gate groove the piers (see Fig. 48) beneath the gate 
which, however, did not affect the position the steel plate; but when the corner 
the plate touched the solid nappe the plate was kicked violently downstream. 

During these tests the tailwater level was that the crest was not 

submerged. 
CoMPARATIVE RESULTS 

The results the tests models and prototype are presented Table 
other data for 2-ft gate opening being plotted Fig. 49. Observations 
the scale model, shown Fig. 49, were also made for gate openings 
and ft. For gate openings 0.25 and 0.5 the comparison was 
between the prototype and model only (see Table 7). 

The pressure heads both model and prototype are fair agreement for 
the gate the downstream position. The maximum negative pressure head 
observed the prototype was —1.2 for gate opening Piezometer 
No. just below the gate seal. The greatest negative heads either model 
had about the same values approximately the same localities. 

For the gate the upstream position there are some interesting departures 
the prototype pressure heads from those the model. Table 7(b) shows the 
comparative pressure heads (see also Fig. 49(b), for 2-ft opening). Unfor- 
tunately, prototype Piezometers Nos. and were clogged. Every effort 
free them with compressed air failed. 

the protype pressure heads for Piezometer No. could have been taken 
for gate openings 0.25 and 0.5 ft, one might assume that conformity could 
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have been demonstrated, but one dares not any farther with that assump- 
tion. For gate openings and ft, the anomalies between prototype 
piezometer heads Nos. and and the corresponding model piezometer heads 
Nos. and became increasingly obvious. 

drawing the line prototype pressure heads between Piezometers Nos. 
and for the five gate openings, one might have the curve 
conform the line obtained the model, but without the model results the 
prototype line would not have been drawn. Moreover, have drawn 
otherwise than indicated the example (Fig. 49), particularly for the gate 
openings and ft, would have made rather curve. 


The explanation for the discordant results undoubtedly lies the fact 
that with the gate the upstream position there tendency for the jet 
spring free from the crest. This tendency creates zone subatmospheric 
pressure just beyond the axis crest the spillway. The maximum observed 
negative pressure head the prototype this zone —5.0 Piezometer 
No. (Table item 26), whereas the projected pressure head Piezometer 
No. (Fig. 49) might have reached —6.0 ft. 

These negative pressure heads appear reach maximum about 
beyond the crest and gradually diminish positive pressures about 
beyond the crest (see Table 7(b), 4-ft opening). This effect doubtless would 
become more marked and reach maximum some larger gate opening, and 
thereafter gradually diminish, evidenced the last column Table 
where all pressure heads are definitely positive. For the prototype this 
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TABLE PRESSURE (IN FEET*) SPILLWAY 


AND PROTOTYPE 


Distance Pie- 
No. from axis Structure 
(ft) 


Prototype 17.8 9.8 5.8 4.7 7.3 9.5 3.7 
Model 0.9 1.3 1.8 2.7 5.4 
Prototype 0.5 1.4 2.5 4.9 6.0 
Model 0.7 0.7 0.9 1.9 4.0 
Model 1.5 1.6 1.6 1.4 1.3 
20.0 
Model 1.3 1.6 2.1 2.5 2.6 
Prototype —0.2 0.6 0.9 1.5 0.9 8.0 
Prototype 0.1 0.4 1.3 1.3 1.3 9.0 
Model 16¢ 0.4 0.4 0.4 0.5 0.5 
50.0 
Model 15¢ 1.0 1.2 1.2 
(b) CENTER OF SEAL 9 Fr UpstrREAM FROM AxIS 
Model 1.2 2.1 2.9 2.5 2.4 
Model 0.4 0.9 1.5 1.4 
0.2 0.3 0.7 1.9 2.1 
Prototype 0.0 0.3 0.3 0.9 8.0 
Model 15¢ 0.3 0.3 0.3 0.4 0.4 


Pressures are given feet water the spillway floor feet water equals one 
Corresponding prototype piezometer clogged. 


See Fig. 49. 


GaTE OPENING, FEET: 


0.25 0.5 1.0 2.0 


(a) CENTER SEAL FROM AXIS 


3.0 4.0 open 
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perfectly normal behavior, and the line pressure heads believed represent 
the true state affairs. Why, though, did not the model experiments disclose 
the same tendency? All model pressure heads with the gate the upstream 
position are the positive side. 

show these anomalies better, Table has been prepared. gives the 
difference pressure heads between prototype and model downstream from 


TABLE 8.—DIFFERENCE BETWEEN PRESSURE HEADS FOR THE GATE 
THE 
(Prototype Value Minus Value Determined from Model) 


Gate 
opening 


0.0 7.5 8.5 9.0 11.0 20.0 0.0 7.5 8.5 9.0 11.0 20.0 


(a) G-4 (c) G-6 
(6) Test G-5 Test G-7 
0.50 —-1.2 | —0.6 0.1 | —0.1 | —0.4 —0.5 | —1.2 0.0 0.1 | —0.3 | —0.6 | —0.5 


Distance from axis dam piezometer, feet. 


the crest for the gate the upstream position. Nearly all differences are 
negative; that is, the prototype pressure heads are lower than the model 
pressure heads, and these differences increase the gate openings become 
larger. 

Fig. shows that model tests G-6 were made with undistorted crest, 
while tests G-4, G-5, and G-7 were with distorted crests. The depressions 
the concrete crest, however, begin the downstream gate seal, which 
below the axis except for G-7 for which crest was raised between the gate seal 
and the axis. Anomalous pressure heads shown Table 7(b) appear upstream 
these crest distortions; hence, these distortions cannot assigned the 
cause the discrepancies. 

The negative pressures resulting from the nappe adhering the underface 
the gate below the seal (40) obviously canndt assigned cause the 
departure the model from prototype behavior for two principal reasons: 
First, the anomalies occur far downstream from the gate and, second, the 
underface both the model and the protype gate was perforated admit air. 
The perforations the prototype gate consist three rows staggered 
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holes, spaced the equivalent continuous slot parallel the 
rubber seal strip and communicating with the atmosphere the downstream 
side the gate. The first row holes in. downstream from the center 
line the seal strip. The perforations through the underface the 1:5 
model gate simulated those the prototype nearly was practicable. 


CONCLUSIONS 


The writers can offer valid explanation for the anomalous results 
between model and prototype with the gate the upstream position. 

Under such high velocities—45 per sec for the prototype and per 
sec for the model—the difficulty obtaining accurate pressure heads through 
piezometers was fully recognized and every precaution was taken eliminate 
such sources error. departures from conformity are not this nature. 
The model consistently refused show negative pressures where such pressures 
should have obtained just below the crest. 
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METER MEASUREMENTS DAM DISCHARGE 


comparison between series current-meter measurements discharge 
the dam presented this paper. Prototype heads covered the measure- 
ments ranged from 1.7 6.3 ft, and discharges ranged from less than 2,000 
nearly 16,000 per sec. Eleven the twenty-three discharge measure- 
ments were made while the dam was submerged. Submergences ranged from 
less than 12% nearly 90%. 

Model discharges averaged approximately less than comparable 
discharges obtained for the prototype. The comparison would appear some- 
what more favorable made terms head. means certain 
that the discrepancy chargeable the model. 

Both model and prototype exhibit considerably smaller diminution 
discharge capacity, due submergence, than was anticipated. The model 
tests demonstrate the inadequacy single-valued relations between sub- 
mergence and reduction discharge capacity. Observations the behavior 
the nappe are presented. 


INTRODUCTION 


Submerged flow, defined the condition that exists when the tailwater 
dam elevation greater than that the crest, often occurs during 
floods low-head dams. Dams designed remain submerged all times 
are used increase depths lakes and rivers. Despite the considerable 
importance the problem, published data submerged flow over dams are 
exceedingly rare. The writer knows previous model-prototype com- 
parison involving this phenomenon. These facts motivated the investigation 
described this paper. 

PROTOTYPE 


Measurement the dam five cross sections revealed that was built 
close agreement with construction drawings. The river bottom upstream from 
the dam, although somewhat variable transversely, can represented fairly 
well the broken line Fig. 50. The dam founded rock which, 
the toe, has been excavated and scoured maximum depth approximately 
lower than the crest. 

The spillway 273.50 long, composed main section 269.50 long 
and section 4.00 long apparently intended for use fishway. These 
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sections are separated wall 1.48 wide that has its top 1.41 higher 
than the spillway crest. profile, the top the dividing wall roughly 
parallels the spillway. The wall subject overflow when the head the 
dam exceeds 1.41 ft. additional length 6.75 subject overflow 
when the head exceeds 3.10 the spillway proper. This additional length 
flat-crest section, 2.00 wide, one end the dam. 


Radius 


Direction 
Flow 


Maximum Height 


Cross 


comparing discharges model and prototype, the measured prototype 
discharge was reduced the estimated flow over the odd sections—namely, 
the dividing wall and the flat-crest. Owing the short length overflow 
for the wall section and the low head the high flat-crest sections, these 
corrections were small, never exceeding the measured discharge. Since 
any reasonable weir coefficient used making the adjustment would give 
substantially the same results, the computation not explained detail, but 
both measured and corrected discharges are given. 

The laboratory head-gate structures act non-overflow section the 
dam long. The gate pier adjacent the spillway well rounded, and 
contraction the nappe can observed the laboratory end the dam. 
The power house forms short non-overflow section the other end the 
dam. The intake crib, trashracks which are parallel the river bank, 
tends reduce end contraction here also. Since even the Francis correction 
for full contraction would have reduced the effective length the spillway 
only about 0.2% the maximum head observed, adjustment for end 
effects was considered necessary. 

All data pertaining the prototype are based field measurements 
made during the course the study. 


— 


OBSERVATIONS PROTOTYPE 


The head the dam was read staff gage mounted retaining wall 
the laboratory bank the river. The gage only upstream from 
line projected along the spillway crest, but its sheltered position, trans- 
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versely from the overflow section, appears keep the readings from being 
affected drop-down. Some concern was felt over this apparently unfavor- 
able location, but simultaneous observations the water surface with 
surveyor’s level the gage and points clearly beyond any possible drop-down 
effect indicated that the location was satisfactory. The situation seems some- 
what analogous the measurement head for contracted weir the side 
There seems reason why the head could not measured 
nearly the plane the weir, well some distance upstream, provided 
that the distance from the point measurement the edge the weir were 
sufficiently great. water approached the head-gate structures with high 
velocity, conceivable that some up’’ might occur which would 
increase the gage reading. Observations indicated the absence any such 
phenomenon the present case. Because bank friction and the alinement 
the river bank upstream from the gage area, the approach water toward 
this vicinity much slower than out the center the channel. 

The tailwater elevation was measured automatic water stage recorder 
350 downstream from the dam. There controlling riffle about 1,200 
downstream which produces flat slope this reach low stages. During 
the high discharges considered this study, ice jam which extended for 
several miles downstream from the dam exerted similar effect. The ice jam 
built the dam number occasions but was never within 500 
the dam during any the measurements. These data downstream condi- 
tions are cited because desired show that the tailwater level would 
have been affected very little changing the measuring point, provided that 
the point was clearly beyond the immediate effects the dam, was the case. 

Discharge measurements were made from cable station 470 downstream 
from the dam. Price current meters good condition were used. Twenty 
the measurements were taken the 0.2 and 0.8 depth method. The 
remaining three were based surface velocity alone. For these, reduction 
coefficient 0.9 was applied determine the mean velocity the vertical. 


(Some the measurements made during the course the study were 


making the station record during the period ice effect.) With due 
regard for the difficulty making accurate measurements under adverse condi- 
tions, the measurements are considered reliable. number attempted 
measurements had abandoned because floating ice. 


The model was scale replica the cross section the prototype 
but did not reproduce the structure regards length, end effects, bottom 
irregularities. The model was made stretching galvanized iron skin 
plate tightly over welded steel frame. The surface presented the flowing 
water was very smooth. 

The model was operated glass-sided steel flume long, 2.5 deep, 
and 2.56 wide, glass glass. The model was bolted the floor the 
flume from its downstream end. End seals were pitch and marine 
putty, smoothed the profile the model. Bottom seals were portland 
cement mortar. The silt against the upstream side the prototype was 
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simulated inclined plane wood. Wooden baffles upstream from the 
model were used produce uniform approach conditions. 
Water was supplied the flume constant rates from the laboratory 

circulating system and was measured sharp-crested, fully suppressed and 
aerated, rectangular weir 2.51 long and 1.29 high. insure accurate 
rating, the weir was calibrated place. This was considered especially 
necessary view the high discharges which taxed the capacity the weir. 
larger weir would have been preferable. For the calibration, by-pass 
flume leading 20,000-lb capacity weighing tank was installed. ettempt 
was made keep the condition the weir crest the same during the model 
tests during the calibration. 

Prior the tests, was determined that the friction slope the flume 
the depths and discharges which were used would not exceed approxi- 
mately 0.0001. Accordingly, the effect friction the flume could 
neglected. This made unnecessary attain exact correspondence points 
which heads were measured model and prototype, again provided that 
measurements were beyond the immediate effects the dam. The head was 
measured 3.1 upstream from the back the dam and the tailwater elevation 
11.3 downstream from the same point. These points were determined 
experimentally being out the drop-down curve and beyond the rising and 
irregular water surface downstream from the dam, respectively. Both water 
levels were observed with point gages the rack and pinion type mounted 
along the center line the flume. 

The discharge over the model was controlled with gate valve 10-in. 
supply line leading the weir tank. The tailwater elevation was regulated 
with hinged flap gate the downstream end the flume. 


The letter symbols used this paper are defined where they first appear 
and are assembled for convenience reference under “Bibliography and 
Notation” the end the Symposium. 


Fig. shows the relation observed for unsubmerged free discharge 
between the coefficient and the energy head for the model the formula: 


this weir formula: the time rate discharge, cubic feet per second; 
coefficient, dimensionally the square root acceleration; the 
length weir dam crest, feet; the measured head dam, referred 
the elevation the crest, feet; average velocity flow; and ac- 


celeration due gravity. The head corresponding the average velocity 
2 


Eq. used characterize the flow whether free submerged. The 
the discharge which occurs when the flow submerged that which 
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would occur under the same energy head for free flow conditions used 
measure the diminution the discharge capacity caused submergence. 

The results model tests for submerged conditions are shown Fig. 52. 
Following usual practice, the diminution discharge capacity due sub- 


Weir Formula (Eq. 13) 


Coefficient 


0.1 0.2 0.3 0.5 0.6 
Energy Head, Feet 


mergence indicated curves which the submerged discharge, expressed 
percentage the unsubmerged discharge which would obtain under the 
same energy head, plotted against the submergence. For submergences 
lower than 45%, the relation between discharge and energy head was not 
affected the submergence the curves are not carried below 45% sub- 
mergence. The total model discharge and the discharge per foot crest 
length are both given these curves, the latter item being more significant 
since the length the model bears only incidental relation the length 
the prototype. 

Fig. 53, the curves Fig. are compared larger scale. inter- 
esting relation exists between the curves for the five higher discharges. Given 
any two the curves, one curve and well-defined point another, the 
remaining three four curves could fixed linear interpolation extra- 
polation terms the discharge along the scale relative discharges. 
haps slightly better fit the plotted points Fig. could have been 
obtained the curves had been drawn without regard for this property. 
seen that the curve for the lowest discharge does not correspond the 
other five. The energy head for this discharge less than 0.1 the model. 
The matter great interest connection with model-prototype relations 
low heads because there some question the behavior the prototype 
for corresponding discharge. regretted that this question could not 
thoroughly investigated. 

When water flows over submerged spillway, the nappe may either plunge 
downward along the face the spillway (as when the flow unsubmerged) 
may flow away from the dam near the surface the tailwater. 
The flowing nappe condition might further subdivided according the 
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appearance the surface. Sometimes the surface has the appearance 
hydraulic jump. other times there series smooth standing waves. 
These conditions are considerable interest, but these tests the only 
distinction made was that between the plunging and flowing nappe conditions. 


100 


Ratio Submerged Free Discharge, Percentage 
uo 


Submergence, Percentage 


The plunging nappe associated with low submergence, the flowing nappe 
with high submergence. For constant discharge, the initial state tends 
maintain itself that there range submergences within which the 
nappe may either plunge flow, depending the direction from which the 
indeterminate zone was entered. 

Fig. shows the results tests made determine the heads and tailwater 
elevations which the changes from one nappe condition the other occurred. 
The observational error was unavoidably high making this determination. 
The results are self-explanatory. The lines drawn are not intended represent 
the points but merely give better idea the submergences which these 
changes occurred. 
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Special tests were conducted determine whether not the behavior 
the nappe had considered studying the effect submergence upon the 
discharge. This was done manipulating the tailwater create 
alternately plunging and flowing nappes for the same submergence and dis- 


Downstream Head, in Feet 


Change from Flowing 
Nappe 


Change from Plunging 
Nappe Flowing 


0.4 
Upstream Head, Feet 


THE NAPPE OBSERVED THE MODEL 


charge. The nappe condition had apparent effect the relation between 
submergence and discharge. 


CoMPARISON MODEL AND PROTOTYPE 


Table comparison discharges measured for the prototype with 
those indicated the model for corresponding conditions. The prototype 
measurements are arranged ascending order the discharge for free flow 
observations and ascending order the submergence for the submerged 
observations. The comparison could made various ways, the method 
adopted being one which gives the result terms discharge. 

previously explained, the length the prototype spillway 273.50 ft. 


model 22.8 long would have been necessary simulate the 


prototype regards length. The model discharges per foot length, there- 


fore, must multiplied 22.8 length adjustment before being “scaled 
up” Froude’s law. 


0.1 
0.2 0.3 
100 
ppe 
the 
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Table 9(a) gives prototype data which appear require explanation. 
Since model and prototype data are separated the table, subscripts are 
necessary differentiate between them. Table 9(b) gives steps the de- 
termination the spillway discharge indicated the model. Only the 


TABLE MODEL AND PROTOTYPE 


E (E Co- Dis- Flow, indi- Di ff 


1,920 | 1,920} 1.70 | 0.01 | 1.71 |.... 0.143 | 0.0541 | 2.86 | .... 10.155] 1,760] — 83 
1,960 1,960 | 1.74 | 0.01 | 1.75 |.... 0.146 | 0.0558 | 2.87 | .... |0.160] 1,820] — 7.1 
2,190 | 2,190 | 1.77 | 0.01 | 1.78 |.... 0.148 | 0.0569 | 2.88 ... {0.164 1,860 | —15.0 
2,260 | 2,260 | 1.87 | 0.01 | 1.88 ].... 0.157 | 0.0622 | 2.90 | .... |0.180] 2,050} — 9.2 
2,410 | 2,410 | 1.96 | 0.01 | 1.97 ].... 0.164 | 0.0664 | 2.92 | .... |0.194] 2,210) — 83 
2,700 | 2,700 | 3.13 | 0.01 | 2.14 |.... 0.178 | 0.0751 | 2.95 | .... |0.222| 2,53 — 6.3 
2,910 | 2,910 | 2.18 | 0.02 | 2.20].... 0.183 | 0.0783 | 2.96 | .... |0.232| 2,650] — 89 
3,470 | 3,460 | 2.52 | 0.02 | 2.54].... 0.212 | 0.0976 | 3.02 .... 10.295] 3,360) — 2.9 
3,840 | 3,830 | 2.64 | 0.03 | 2.67].... 0.223 | 0.105 3.04 | .... |0.320|, 3,640) — 5.0 
4,500 | 4,490 | 2.96 | 0.04 | 3.00 ].... 0.250 | 0.125 3.08 | .... |0.885] 4,380] — 2.4 
6,180 | 6,160 | 3.57 | 0.06 | 3.63 ].... 0.303 | 0.167 3.17 {0.529 6,020] — 2.3 
6,580 | 6,560 | 3.69 | 0.07 | 3.76 |....] .... 0.313 | 0.175 3.18 ... |0.557| 6,340 | — 3.4 
3,450 | 3,440 | 2.39 | 0.02 | 2.41 |0.28| 11.7 0.201 | 0.090 3.00 | 1.000 |0.270 | 3,070 | —10.7 
9,710 | 9,650 | 4.59 | 0.12 | 4.71 |1.92] 40.8 0.393 | 0.246 3.30 | 1.000 |0.814 |- 9,270 | — 3.9 
12,600 {12,500 | 5.44 | 0.18 | 5.62 |2.37| 43.6 0.468 | 0.320 3.40 | 1.000 |1.089 | 12,400 | — 0.7 
12,6002 |12,500 | 5.19 | 0.19 | 5.38 |2.45| 47.2 0.448 | 0.300 3.37 | 0.998 |1.008 | 11,500 | — 7.9¢ 
9,740 | 9,680 | 4.69 | 0.12 | 4.81 12.57] 54.8 0.401 | 0.254 3.31 | 0.992 |0.834| 9,500 | — 1.9 
3,480 | 3,470 | 2.65 | 0.02 | 2.67 |1.52] 57.4 0.223 | 0.105 3.04 | 0.993 |0.318 | 3,620] + 4.3 
15,9002 |15,800 | 6.03 | 0.26 | 6.29 |3.62| 60.0 0.524 | 0.379 3.47 | 0.980 |1.289 | 14,700 | — 6.9¢ 
14,7502 |14,650 | 6.28 | 0.21 | 6.49 13.84] 61.2 0.541 | 0.398 3.49 | 0.980 |1.359 | 15,500 | + 5.9¢ 
6,600 | 6,580 | 3.65 | 0.07 | 3.72 |3.03| 83.0 0.310 | 0.173 3.18 | 0.965 |0.530 | 6,030 | — 8.4 
6,930 | 6,890 | 3.98 | 0.07 | 4.05 |3.48| 87.5 0.338 | 0.196 3.22 | 0.930 |0.589 | 6,700 | — 2.7 
7,380 | 7,340 | 4.06 | 0.07 | 4.13 |3.62| 89.2 0.344 | 0.202 3.22 | 0.912 |0.592| 6,740] — 81 


Measurement surface velocity with 0.9 reduction; other measurements, 0.2 and 0.8 depth. 


energy head and the submergence this part the table are based proto- 
type observations. brief explanation the computation procedure seems 
desirable. The energy head for the model obtained dividing the corre- 
sponding prototype item 12. Thesubmergence, being numerically 
equal for model and prototype. 

For free flow, when the submergence less than 45% (below which 
value the model indicated that submergence had effect), the discharge 
per foot length the model, computed from the formula, 


_Q_ 1.5 


the coefficient being taken from the curve Fig. 51. 
When the submergence exceeds 45%, the discharge computed Eq. 
must reduced the curves Fig. 53. The model shows, however, that 


_ | 
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the ratio submerged free discharge the same energy head (hereinafter 
called simply the function not only the submergence 
but also the discharge, which yet unknown. Accordingly, the final 
value the discharge ratio cannot found until the discharge known and 
vice versa. solution can found few easy approximations. 


simple procedure assume that the discharge ratio has the value 
100%. first trial value can now found, and second trial value 
the discharge ratio can read from Fig. interpolation. This gives 
second trial value The process repeated until remains unchanged 
further approximation, which will after very few trials. Obviously, the 
procedure can shortened close estimate the discharge 
ratio made initially. trial computations are shown Table but the 
correctness the results can verified without approximation. easier 
method calculation could derived large number discharges were 
determined. 


Rating Curve for Unsubmerged 
from Model Tests 


Prototype Measurement Unaffected 
Submergence; Comparable with Curve 

Prototype Measurement Affected 
Submergence; Comparable with Segment 


Energy Head, in Feet 


Thousands Cubic Feet per Second 


Cols. and 14, Table are comparative data. The spillway discharge 
indicated the model the product three factors—the discharge per 
foot length the model, the length adjustment previously explained, 
and the five-halves power the scale ratio. Col. 14, Table gives the 
percentage difference between the measured net spillway discharges and the 


A- 
Submerged Rating Segment 
Obtained from Model 
ally 
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spillway discharges indicated the model for corresponding energy heads 
and submergences. The difference considered negative when the indicated 
discharge less than the measured discharge. 

Fig. the performance the model compared graphically with 
prototype measurements. The continuous curve energy-head-discharge 
relation for free flow over the spillway, established from tests. 
Observations the prototype when unsubmerged, when the submergence 
did not exceed 45%, are comparable directly with this curve. 

avoid confusing the diagram, continuous curves were not drawn for 
submerged flow. ‘For each submergence which measurement was made, 
short segment the submerged rating indicated the model tests was 
drawn the corresponding energy head. Since submergence reduces the 
discharge for given energy head, the submerged rating segments all lie above 
and the left the curve for free flow. should noticed that the three 
measurements which show the greatest discrepancies are based only surface 
velocity. close agreement can expected for these three measurements. 


CoMMENT 


Effect Submergence the the best known tests 
submerged flow over rounded-crest dam are those performed for the 
Board Engineers Deep Waterways Cornell University Ithaca, 
(41)(42), 1899. Results these tests are frequently quoted. Draw- 
ings the dam which the tests were made are quite generally available, 
illustration included this paper. 

part the relation between submergence and the ratio submerged 
free discharge obtained another dam shown for comparison Fig. 53. 
seen that the effect submergence obtained the present study was only 
small fraction that obtained the dam question. The effect sub- 
mergence the model second dam was not appreciable until the sub- 
mergence exceeded 45%, whereas the first dam was affected much lower 
submergences. 

The comparison presented merely with the idea showing that the use 
submerged flow data can nothing but rough approximation. 
The use the data one dam estimate the effect submergence the 
other would have caused greater errors than would have occurred entirely 
neglecting the effect submergence. 

Analysis the Effect data quoted the foregoing 
section plot smooth curve. This fact appears have given rise 
prevalent belief that, general, the effect submergence the discharge 
over dam can represented single-valued curve relating submergence 
and the discharge ratio. Attempts plot submerged spillway data this 
manner invariably lead wide scatter observation points. This, turn, 
leads the opinion that submerged flow over dams far more unstable 
phenomenon than really the case. The submerged flow tests the first 
dam were made constant head 6.6 ft. made the tailwater elevation 
and the submergence functions the discharge, which had decreased 
the tailwater was raised keep the measured head constant. Study 
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Fig. shows that, set observations were plotted, which the sub- 
mergence was functionally related the discharge, single curve rather than 
family curves would express the results adequately. 

stable control exists downstream from dam, the submergence will 
indeed related the discharge this special way. model tests are made 
for such case, the tailwater must regulated represent the prototype 
condition single-valued relation between submergence and its effect 
the discharge desired. This procedure not recommended. unwilling 
vary the tailwater independently the discharge matter curiosity, 
the experimenter should consider the possibility that the control might change 
and invalidate the results these results were restricted scope. The ice 
effect which caused the unusual submergence observed this study 
pertinent example. 

Measurement Tailwater Society’s Manual Engineering 
Practice No. 11, entitled ‘‘Letter Symbols and Glossary for 
defines submergence follows: 


“The ratio the tail-water elevation the head-water elevation, 
when both are higher than the crest, the overflow crest the structure 
being the datum reference. The distances stream down stream 
from the crest which head-water and tail-water elevations are measured 
are important, but have not been standardized.” 


connection with measurement headwater elevation, there seems 
reason for departing from the procedure commonly used for unsubmerged 
weirs and dams—namely, the measurement point sufficiently distant from 
the overflow safely out the drop-down curve. 

Cox, Am. Soc. (45), measured tailwater elevation very near 
thedam. this point, especially when the nappe plunges, there decided 
depression the water surface. for the flowing nappe condition, the 
water level near the dam lower elevation than farther downstream. 
The reason for negative slope this reach readjustment velocity 
distribution which has effect that similar, regards energy transfer, 
the recovery pressure head which occurs the diffusing cone venturi 
meter. The investigation Professor Cox dealt primarily with the use 
the submerged weir measuring device. For this purpose, the use 
tailwater measuring point very near the dam apparently satisfactory. 

For most purposes, however, believed that tailwater measuring point 
clearly beyond the rising water surface decidedly preferable. Professor Cox 
found necessary consider the behavior the nappe and gave different 
formula for each condition. This complication that should avoided 
possible. recording gages are used, the behavior the nappe cannot 
determined definitely from gage readings alone the submergence the 
critical zone. times during observations the second dam, plunging 
existed over part the spillway length and flowing nappe over the 

The model tests indicated that the tailwater level measured 
sufficient distance from the dam, the behavior the nappe need not con- 
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sidered. This just another way stating that the over-all loss head 
flow over the model did not depend any measurable extent the behavior 
the nappe. 

seems best, also, consider the rising water surface downstream from 
the dam part the phenomenon submerged flow. Otherwise, the 
magnitude the rise water level must estimated fitting results 
model tests stage-discharge relation downstream from the dam site. 
satisfactory method making this adjustment available. probable 
that model tests indicate its magnitude more accurately than can 
computed. 

Comparison Model and the whole, the agreement between 
model and prototype considered satisfactory. Not much weight can 
given the three measurements made surface velocity. the coefficient 
0.85 had been used for these measurements instead 0.90, two them would 
agree very well with the model. coefficient 0.85 well within the range 
possibility, but the originally estimated value 0.90 has been retained. 
may seen from Fig. 55, the agreement between model and prototype except 
for the three surface velocity measurements about good that between 
the discharge measurements themselves. 

For estimating flood discharges from water level observations dams, 
would appear that the use small-scale sectional models must regarded 
being satisfactory the current meter—from the standpoint efficiency, 
far more so. case would estimate the backwater due the dam, 
based the model tests, have been error more than few inches. 

For all but two the twenty-three measurements, the model indicated 
discharge smaller than that measured the prototype. what extent this 
difference chargeable disproportionately large viscosity, surface tension, 
friction, and end effects not known. Overregistration the current meters 
may have accounted for some the difference. 

previously suggested, the nonconformity the curve for the lowest 
discharge Fig. the general pattern for higher discharges special 
interest. seems possible that the model fails for discharge low this, 
but cannot disproved substantiated from data now available. 


The conclusions presented are necessarily tentative. However, they are 
considered constitute logical interpretation the results model tests 
and measurements the dam. 


(a) The effect submergence the discharge affected the shape 
the spillway and perhaps other factors. Data present available are 
insufficient for generalization the effect submergence terms spillway 
shape. The use data cannot 

general, not possible express the effect submergence upon 
the discharge single-valued relation between the submergence and the 
discharge ratio. the range covered these tests, the effect submergence 
greater for high discharges. 
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(c) the tailwater elevation measured sufficient distance from the 
dam, not necessary consider the behavior the nappe determining 
the effect submergence upon the discharge. For this and other reasons, 
the tailwater elevation should measured beyond the immediate effects 
the dam. 

(d) Reasonably good agreement may expected between spillway models 
and prototypes the sizes used these tests under submerged well 
free flow conditions. 


other observations might presented secondary incidental 
conclusions, but these appear have been covered the extent that their 
importance justifies the body the paper. 
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DEVELOPMENT MISSISSIPPI 
RIVER CHANNEL 


The Waterways Experiment Station recently received from the 
Engineer Office, St. Louis (Mo.) District, field data that served confirm the 
results obtained from hydraulic model study (53) (54) conducted the Experi- 
ment Station during 1936 and 1937 reach the Mississippi River. 
careful analysis these data has been made and comparison the action 
the prototype with the indications the model furnishes the basis for 
this paper. 


THE PROBLEM 


The model study involved the development effective plan for improv- 
ing the navigation channel the Mississippi River. certain point the 
Mississippi River enters rock gorge, the bluffs which confine the flow 
the river preceding section. the latter point the rock bluffs the 
left bank the Mississippi River recede, allowing the river widen and 
divide its flow. The natural tendency the river follow the chute 
channel. However, because the more satisfactory navigation conditions 
that would prevail the river could made develop the crossing the 
right bluff line, two pile dikes were constructed 1939 from the left bank. 
These dikes effected considerable improvement the navigation channel, but 
did not eliminate the need for periodic maintenance dredging through this 
crossing. The material dredged from the crossing was spoiled front of, and 
off the end of, one dike and the chute behind this dike attempt 
confine the river flow the crossing. Samples obtained after the completion 
the model study indicated that the river-bed material the vicinity was 
composed sand, gravel, and coarse rock fragments. 


The scale ratios, model prototype, were follows: Horizontal dimensions, 
600; and vertical dimensions, 150. The model was the movable bed 
type, with banks and overbank topography molded concrete and the movable 
channel section molded crushed coal. addition the slope scale 
induced the difference scale ratios, the model was given slight additional 
slope sufficient provide adequate movement the bed material used. The 
scales used were within the range which has been found that velocity distri- 
butions and current directions can reproduced accurately models wide, 
shallow rivers such the Mississippi. 


(Deleted censors.) 
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After thorough adjustment and verification the model, which operat- 
ing conditions were developed such that the movement bed material the 
model was consistent with that the prototype, tests plans were begun 
The testing procedure consisted installing the elements proposed im- 
provement plan the model, the bed which had been molded conform 
the 1936 prototype conditions, and operating the model accordance with 
the stages modified form the 1935 hydrograph (Fig. 56). The 1935 


Fia. AND PROTOTYPE 


hydrograph included within its range adequate distribution critical 
high stages well normal and low stages. Therefore, since was believed 
embrace all conditions flow likely occur the reach question the 
future, this hydrograph was selected for use the model study. Both the 
dredging proposed the improvement plan and the necessary maintenance 
dredging were simulated the model. 

was thought advisable begin testing the improvement plan 
stage comparable that expected the prototype the time completion 
the prototype improvement. was believed that this completion would 
occur November, 1937. Consequently, operation the model was started 
with the 15-ft stage November (modified 1935 hydrograph). Although 
dredging the initial pilot channel the prototype was not actually completed 


until January 30, 1938, the effect the difference time was not considered 
important. 


AND PROTOTYPE IMPROVEMENT PLANS 


The plan included the shortening one dike and the addition 1,500-ft 
trail dike downstream from this dike, dredge cut 250 wide depth 
(low-water reference plane) along the right bank island, use 
the dredge spoil chiefly build the left bank the dredge cut, and riprap 
and mattress protection for the left bank the dredge cut and for the trail dike. 

the prototype, the outer 600 one dike was removed, the 1,500-ft 
trail dike was constructed, and revetments were placed along the left bank 
the new channel during the period from August, 1937, January, 1938. The 
pilot cut for the re-alined channel was made during several dredging periods— 
namely, September 28, 1937, January 30, 1938; March May 28, 1938; 
and August October 31, 1938. This cut was more than 400 wide the 
prototype contrast the 250-ft wide pilot cut the model. under- 


1937 1938 1939 
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stood that the dredge spoil (about 3,400,000 yd) was disposed generally 
accordance with plans followed the model test. addition the works 
scheduled the improvement plan, prior the August, 1939, prototype sur- 
vey, the dike the right bank was extended about 500 across the old 
channel. This extension was intended aid the filling the old channel. 


general, the prototype reproduced accurately the results indicated 
the model study. should stated that some divergence action 
expected, because, the time the August, 1939, prototype survey, the dura- 
tion the high stages (above ft) the prototype was less than half the 
duration like stages the model (Fig. 56). Available data indicate 
greatly divergent tendencies present between model and prototype. The 
district engineer reports that November, 1939, the maximum available 
depth through the old channel was only ft, whereas project depth existed 
through the new channel, showing that the prototype conforming sensibly 
the indications the model. Careful examination prototype and model 
data reveals the following interesting comparisons: 


(a) Until the time maintenance dredging January, 1939, the re-alined 
channel the prototype exhibited persistent tendency This 
agrees with the model results, since the end Run the model channel 
had shoaled controlling depth below the low-water reference plane. 

(b) Since the prototype maintenance dredging January, 1939, the scour 
the re-alined channel indicated the model has occurred also the proto- 
type, although somewhat slower rate. This difference rate scour 
may explained the fact that the plan developed the model was based 
the assumption that the river bed this vicinity was composed average 
river sand containing only small amount larger material, whereas actual 
spoil samples (obtained when the actual improvement was undertaken) were 
found contain 40% rock fragments and large gravel. 

(c) tendency scour was noted the model the upper end the 
old channel along the right bank, and tests indicated that this channel would 
continue fill after each high stage. Thus far, alternate scouring and filling 
has also occurred the prototype, with scouring predominating. 
noted that development the re-alined prototype channel not yet complete. 
believed probable, however, that, once this new channel has fully de- 
veloped, the old channel will become filled substantially predicted the 
model. 


The comparison model and prototype behavior the reach the 
Mississippi River involved will supplemented from time time addi- 
tional prototype data are made available. 
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HYDRAULIC CONDITIONS APPROACH 
LOCK AND DAM 


The new lock and dam here considered located the Monongahela 
River. Because local conditions, the structure, originally designed, 
was subjected adverse currents that had effect the maneuverability 
river vessels the upper approach the locks. moderately high stages, 
the main currents, after sweeping around two comparatively sharp bends up- 
stream, followed the lock side the river, increased velocity they reached 
the upper approach, and then swung sharply across the lock entrance toward 
the center the stream before discharging over the spillway section. Exami- 
nation the site and the current conditions disclosed that the flow charac- 
teristics might have been due several factors, combination any 
these—namely, (a) the effect the upstream bends velocity and current 
distribution; (b) the effect the hydrography the vicinity the approach 
velocities and currents; and (c) the effect contraction the stream 
cross section due the lock and the existence large bar along the shore 
opposite the lock chamber. 


Top Bank 


Bar Removal 


RIVER Guard Wall Extension 


Guide Wall 
Level Pool 737.5 


(Crest 737.5) 


Spur Dike 


60' 


The existence the foregoing condition prompted investigation 
made means hydraulic model tests study the causes the current 
phenomena the upper approach and determine the most suitable means 


improving the entrance conditions. the basis these studies, necessary 
(Deleted censors.) 
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revisions were made the upper approach the lock and the stream 
channel above the dam. These involved extension the guard wall, the 
construction spur dike upstream from the lock, and the removal 
portion the bar opposite the lock chamber. 

determine the accuracy and the fidelity with which the model simulated 
the prototype conditions, certain field tests were conducted. The results 
these investigations showed that, generally, hydraulic similarity between the 
model and prototype existed. far can ascertained from the experience 
navigators traversing the upper approach the lock chamber, the flow 
conditions have been considerably improved since the modifications the 
upper pool were effected. 


DESCRIPTION THE STRUCTURE 


The dam about 550 long from the left abutment the lock, and 
the concrete gravity type. uncontrolled, ogee spillway section, 
with crest El. 737.5 (mean provided. 

accommodate the heavy river traffic this vicinity, two locks having 
lift 10.6 are provided along the right bank the stream. The land lock 
720 ft, and the river lock 360 ft. The tops the guide and 
guard walls are El. 749.0. altered produce suitable entrance condi- 
tions the locks from above, the middle lock-wall and the guide and guard 
walls extend approximately 100 ft, 760 ft, and 855 ft, respectively, upstream 
from the upper gates the land lock. addition, spur dike provided 
point 1,600 upstream from the end the guide wall. The dike normal 
the flow the river and extends riverward from the shore point line 
with the river face the guide wall. consists shell steel sheet piling, 
backfilled for strength and stability. The general details the structure, and 
the location the component parts, are shown Fig. 57. 


Approximately 2.4 miles the river were reproduced the model, the 
major part the reach being upstream from the structure develop ac- 
curately the flow characteristics the upper approach. The model was un- 
distorted, scale 200, and was the fixed-bed type built according 
standard laboratory practice. All necessary details that might affect flow 
phenomena were reproduced faithfully scale, including the lock and spillway 
and the topography the immediate vicinity. The assumed channel roughness 
was reproduced accurately was physically practicable. 

All basic tests for the design the required improvements were made with 
the upper pool 0.5 below the top the lock-walls. The flow corresponding 
this condition was assumed represent the greatest discharge under which 
the navigation facilities could used. The tests consisted analyses the 
resulting velocities and currents, and their effects tows entering the locks 
from above. The investigations involved determination the flow condi- 
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Spur Dike 


Top Bank 


Center Line 


(a) DIRECTION AND VELOCITY, MODEL AND PROTOTYPE 


General River Current Genera 


General River Current 


(c) DIRECTION CURRENTS, MODEL 
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General River Current 


Guard Wall Extended 


> 


Guide Wall 


(Approx.) 


RIVER MODEL 


loci 
per Sec per Sec 


1230' Crest Dam (Approx.) 


General River Current Guard Wall Extended 540' 


Guide Wall 
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SURFACE CURRENTS UPPER APPROACH PROTOTYPE 
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tions for the original design, and the separate and combined effects due 
extension the guide and guard walls, dredging the upper lock approach, 
dredging the bar opposite the lock chamber, and deflection the flow 
means dikes various proportions and located various positions. 

The data obtained from the model were three types: 


(a) Notes visual inspection the flow conditions; 

records confetti movements and the paths traversed 
model tows; and 

(c) Actual measurements the surface velocities the critical areas. 


the analysis the test data for the several proposed designs, the photo- 
graphs disclosed the alinement, direction, and velocity the currents the 
alinement, direction, and length the streaks. addition, the photographic 
record revealed the effect these currents the travel model tows that 
were allowed glide toward the lock without artificial maneuvering. For 
true comparison the data the various designs, the latter procedure was 
standardized the use tow-projecting device that imparted each tow 
the same initial velocity and direction. Samples the photographic data for 
the original design are shown Figs. and surface velocity measure- 
ments were made the region where they were believed have the greatest 
bearing the maneuverability tows. This was accomplished spanning 
the critical region with piece plate glass and marking thereon the location 
floating paper disks regular intervals time. 

provide data for comparison with the prototype results, and thereby 
furnish means verifying the model, special tests were run the model for 
both the original and altered conditions the prototype structure. These 
data were the form photographs the current movements and analysis 
the surface velocities the upper approach, the procedure during these 
tests being the same the foregoing. 


The details the structure, both the model and the prototype, were the 
same, that exact comparison could made for verification, except for 
the effects surface tension and channel roughness—factors which were im- 
possible exact simulation the model. Field tests were made both before 
and after the alterations had been made the prototype. 

The field data were necessarily confined measurements current direc- 
tions and surface velocities. data tow movements, affected the 
flow the approach, were obtained, since the artificial guiding the tows 
the prototype could way simulated tow movements the model. 
The current directions and velocities were obtained tracing the paths 
floats the river and recording their locations regular intervals time. 
This was accomplished spotting the floats simultaneously with three 
transits located strategic points. Corrections were applied the 
measurements for wind velocity noting the wind effect floats the lock 
chamber. 
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VERIFICATION RESULTS AND CONCLUSIONS 


comparison the data obtained the model and the field, the form 
current alinement and velocity vectors, presented Figs. and 61. For 
the original design the approach (Fig. 60), the operating conditions the 
model and prototype were identical, corresponding head 5.6 the 
crest the dam and discharge 30,000 per sec (headwater elevation, 
743.1). Operation for the revised design the upper approach was not iden- 
tical the model and prototype. The head the crest the dam for this 
condition was 6.6 the prototype against 6.0 the model, correspond- 
ing flows 40,000 and and 35,000 per sec, respectively (headwater 
elevations were 744.1 and 743.5 ft, respectively). The model tests had been 
conducted before opportunity had presented itself for obtaining field 
measurements approximating the conditions tested the model for the altered 
design, that identical operation was not possible. However, the differences 
head and discharge are not believed sufficient magnitude invalidate 
the results the verification. 

The degree similarity between the model and prototype results apparent 
from examination the data presented. Exact simulation the model 
the prototype phenomena was not attained; nor was expected. Generally, 
the model data reveal close approximations for the magnitude the prototype 
velocities, for the alinement the currents, and for the nature and location 
the quiescent pool formed downstream from the spur dike. comparison 
the data also shows considerable reduction the velocities flow across the 
lock entrance and the head the guard wall, this condition existing both 
the model and prototype. 

the course the model and field tests, was noted that slight surges 
the flow occurred, especially regions where attempt was made alter 
the flow characteristics, where obstruction flow existed. Because 
these surges, the resulting path flow through any definite point for successive 
intervals time was not identical, although the general pattern flow for any 
area was approximately the same. slight variation the data should 
expected, therefore, both the model and prototype, and comparison the 
data should made the basis general values and flow patterns rather 
than the specific analysis particular flow line. Generally, therefore, 
the velocities, currents, and eddies the model and prototype were similar. 


CONCLUSION 


can concluded that satisfactory verification the model tests for the 
entrance conditions the upper approach the new lock and dam re- 
sulted. Subsequent experience the maneuvering tows the reach the 
stream question has indicated devoid any serious hazards and 
satisfactory every respect, which further attests the success the plan 
alteration proposed result the model tests. many other in- 
stances, model analysis this case provided direct approach problem 
that could not solved purely analytical methods. 
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II. 


The following symbols, used the Symposium, conform essentially with 
“Letter Symbols for prepared Committee the American 
Standards Association with Society representation and approved the Asso- 
ciation 


area cross section, subscript denoting the section location: 
area lock chamber; 
length weir crest; 
discharge coefficient: 
coefficient determined model tests; 
diameter; 
friction coefficient, subscripts and referring model and proto- 
type, respectively; 
acceleration due gravity; 
head above given datum plane: measured tailwater eleva- 
tion, referred the elevation the crest; 
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pressure heads selected sections indicated appropriate sub- 
scripts: 
loss head due friction; 
pressure head differential between two sections; 
head, velocity; 
pressure any piezometer the bellmouth; 
discharge coefficient defined Eq. 11; also constant 
observed constant the rating equation (Table for 
flow meter; 
Scobey’s friction coefficient; 
theoretical constant the rating equation (Table for 
flow meter; 
lengths, such center-line distances from the bellmouth, with sub- 
scripts denote section locations; the distance required pro- 
duce normal distribution below intake smooth pipe: 


7 


constant exponent defined Eq. 11; 


pressure, pounds per square foot; pressure head, 


water; 

time rate discharge, subscripts and denoting model and proto- 
type, respectively; 

discharge per unit width; 

Reynolds’ number; 

radius; 

average velocity given cross section; 


lock-chamber rate rise fall; 
Weber number; 
specific weight water pounds per cubic foot; pressure 
head, feet water; 
central angle. 
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PAPERS 


DISTRIBUTION GRAPHS SUSPENDED- 
MATTER CONCENTRATION 


Damages benefits resulting from sedimentation various points 
drainage basin may due different types sediment, possibly coming from 
different sources. For instance, flood and drainage damage caused stream- 
channel aggradation appears primarily bed-load problem, involving the 
deposition relatively coarse material. the other hand, fine material 
responsible for much the deposition irrigation and drainage canals and 
the increased cost water purification many public water-supply systems. 
Both fine and coarse material contribute the deposits reservoirs, but the 
finer grade sediment normally the chief source damage, because 
constitutes the greater portion the load. Because different types damage 
are caused different proportions the two classes sediment, 
important sediment-load investigations that, not only the total load the 
stream determined, but also the relative proportions fine and coarse 
material. These are the subjects discussed the paper. 


GENERAL 


With the exception dissolved matter, the total solids load passing 
particular cross section stream may classified either bed-material 
load wash load, with definite grain size representing the division between 
the two classes (3).2 The bed-material load, for brevity, 
that part the total solids composed the relatively coarse material (coarser 
than the division grain size), and moves rolling sliding along the bed, 
long jumps between points contact with the bed while temporary 
suspension, rate related the stream discharge. The wash load, except 
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for negligible percentage, moves suspension rate that bears direct 
relation discharge and composed the relatively fine material (finer than 
the dividing grain size). The wash load constitutes much the greater part 
the material suspended-load sample taken flood stage, and may 
regarded sediment relatively permanent suspension. 

The particular grain size that divides the sediment load into these two classes 
varies different streams and, generally, throughout the length particular 
stream, although has definite value particular reach For 
example, the dividing grain size the Colorado River near Imperial Dam 
0.01 (fine silt) (10), whereas the Enoree River near Greenville, C., 
the size 0.351 (medium sand) (3). The dividing grain size can 
determined from mechanical analysis curves samples taken from the stream 
bed immediately after flood, and the finest grain size existing appreciable 
quantities. 

The bed material (coarse material) stream system ordinarily changes 
slowly composition over long period years, all, and the quantity 
moved mainly response discharge. The quantity wash material 
(fine material), however, may vary greatly from year year and from season 
season, because the changing conditions supply that are independent 
conditions flow. The most important variables these completely inter- 
related conditions supply are vegetal cover, tillage methods, stage soil 
erosion, and rainfall characteristics. 

Observations made the Soil Conservation Service (SCS) sediment-load 
laboratory near Greenville (2) have shown the relative magnitudes the two 
types sediment load representative Piedmont stream. the total 
sediment load passing the station the Enoree River during its two years 
operation, more than 90% was found consist wash material. 

Although this ratio between the annual amounts the two classes sedi- 
ment load not expected the same for all streams, even constant for 
the Enoree River, probably represents fairly well the general order mag- 
nitude for streams this size the humid sections the United States where 
the amount wash material subject erosional attack essentially unlimited. 


Vertical Distribution—A general theory the vertical distribution 
suspended sediment turbulent stream water has been established 
result extensive investigations the transportation sediment, both 
the United States and Europe. Based upon this theory, the total suspended- 
sediment load transported vertical section can estimated from deter- 
mination the concentration and the particle size distribution the sediment 
single point the vertical, and the hydraulic characteristics the stream. 
actual practice, however, the procedure more simple because necessary 
apply the theory only sediment coarser than 0.075 (200-mesh), 
since material finer than this size appears uniformly distributed through- 
out the section most streams, and single accurate sample taken any 
point the stream gives the concentration applicable the entire cross 
section. 
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Sampling and Analysis.—By accurate sampling and grade-size analysis 
the suspended matter stream, possible compute the total load 
coarse material—material sizes found appreciable quantity bed deposits 
temporary suspension; and, making sufficiently large number 
observations over wide range stage, relationship between coarse load 
and discharge can established for particular location. Estimates the 


suspended load coarse material, therefore, can made from this 


established relationship, and from record the discharge, without the neces- 
sity further sampling. The load fine material, the other hand, can 
estimated accurately merely from dip sample. However, bears 
direct relationship the discharge, and the past continuous sampling has 
been necessary estimate the total volume fine material carried stream 
during flood. 

Observation Periods.—In small flashy streams almost all the annual load 
wash material carried during few days out the year (6), and the 
content wash material may vary many 100% within the period hour 
during the rising and falling stage the stream. the present time, 
method wide applicability, except frequent sampling, has been developed 
for estimating computing, with even approximate reliability, the total sus- 
pended load carried any small stream. true that some large rivers 
approximate relationship (1) has been developed between discharge and 
suspended load—presumably including both fine and coarse material—and, 
from this, estimates total sediment transportation have been made. How- 
ever, where such relationships have been developed, the period observation, 
general, has extended over only one two years when conditions supply 
apparently were approximately constant. Continued observations over 
longer period time usually have resulted increased scattering data 
(9). The discussion follow concerned primarily with the transportation 
the fine material small streams humid climates and gives general de- 
scription some the sources and factors controlling the supply material. 

Sources Sediment.—The load fine material stream appears have 
its main source, not the stream bed, but other parts the drainage basin. 
possible, course, that small amount fine material any particular 
flood may come from the stream bed, where was deposited pools during 
low stage, where became trapped between larger bed particles, and was 
released when the bed material started move during higher stage. 
general, however, the major sources sediment are cultivated lands, road 
cuts, roadside ditches, gullies, and caving stream banks. The amount 
material that may enter stream from these sources during particular storm 
depends upon such complex, interrelated, and, part, changing factors 
topography, soils, shape and size the drainage basin, land use, vegetal cover, 
and surface and underground storage, well upon rainfall intensity and 
duration. any drainage system, however, certain these factors may 
considered constant, whereas others have seasonal variation. 

Results Past investigations soil loss and runoff 
from relatively small agricultural plots give general information the prin- 
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ciples underlying the supply erosional debris. These studies have shown, 
general, that the sediment concentration runoff water from plot usually 
increases maximum within very short time, and then decreases the 
duration application water continues, approaching constant minimum 
value asymptotically. The high initial sediment concentration appears 
the direct result rain impact. However, rain packing and adhesion from 
wetting and the swelling colloids, the rainfall continues, well the 
variation soil porosity with depth, appear the main factors tending 
reduce concentration. The variation the infiltration capacity with duration 
rainfall also has been observed follow law similar that sediment 
concentration. 

Erosion rain impact and overland flow applies 
more specifically the land surface and the exposed side slopes road cuts, 
ditches, and gullies, but once water concentrated relatively large, well- 
defined channels, the principles stream erosion more nearly apply. Stream 
erosion probably closely related discharge, because, rule, the greater 
the discharge, the greater will the velocity, and hence the capacity dislodge 
and transport particles. any event, however, and despite the complicated 
mechanics which material placed suspension and transported the 
outlet the ordinary small plot, the concentration fine material appears, 
general, vary with the surface runoff. The same reasoning should apply 
natural stream; although, obviously, the problem more complicated, 
because rainfall distribution and areas sediment production, which are 
sense system small plots, vary considerably their characteristics 
throughout the drainage basin. 

Hydrographs and Distribution Graphs.—It has been observed generally that, 
following rainfall during which surface runoff occurs, the discharge this 
surface runoff into stream channels simultaneously over drainage basin 
results pronounced rises stream levels, followed periods decline. 
The concentrations suspended matter streams also have been observed 
rise and decline similar manner (Fig. 1). 

The plotted graph stream flow, showing the rises and declines runoff, 
has been used many investigators basis for the phenomenon 
surface runoff. Sherman, Am. Soc. (8), conceived the idea 
that, for particular drainage area, the surface runoff from rainfall occurring 
within the same time interval will produce hydrographs, the ordinates which 
will vary with the amount the surface runoff. Mr. Sherman has given the 
term the hydrograph surface runoff that results from rain 
falling within unit time, day hour. the unit hydrograph 
surface runoff modified show the proportional relations its ordinates 
the percentage total surface runoff, the resulting graph termed 
tribution The unit-hydrograph method based upon the hypothesis 
that particular drainage basin surface runoff from rainfall occurring 
unit time will produce hydrographs approximately equal bases, and the 
ordinates will vary with the intensity the net rainfall. 

From the consideration, suggested plot studies, that surface 
and suspended-matter concentration appear related the same factors, 
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investigation was made part the SCS sediment-load studies deter- 
mine any such relationship actually existed natural stream. Data from 
the East Fork Deep River, near High Point, were used for this study, and 


Concentration, in Parts per Million by Weight 
Discharge, per Sec 
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upon examining the hydrographs and corresponding graphs suspended- 
matter concentration for all isolated stream rises was found that there was 
“distribution graph’’ concentration analogous the distribution graph 
used analyzing stream flow discussed previously herein, and that could 
used advantage studying the suspended load streams. 


DEVELOPMENT DISTRIBUTION GRAPH 


The data used this study were obtained the East Fork Deep River 
during extensive suspended-matter sampling program conducted the 
United States Geological Survey cooperation with the SCS and the North 
Carolina Agricultural Experiment Station. Collection water samples for 
the determination sediment concentration was made once day during 
periods low water, and more frequent intervals during periods storm 
runoff. Fig. illustrates typical example stream rise, and shows the 
usual frequency sampling considered necessary obtain sufficient informa- 
tion for the determination variations content suspended matter. Other 
basic data secured included the stream-discharge and precipitation records 
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from three standard and three recording gages located and about the 13.9- 
sq-mile drainage basin. 

The basic data from the gaging station are sufficiently complete give 
valuable general information the movement suspended matter typical 
stream the Piedmont. The concentration suspended matter, recorded 
and used the following analysis, includes both fine and coarse material. 
few mechanical analyses samples taken flood stages, however, showed 
that the amount coarse material suspension was insignificant compared 
the fine material that comprised most the wash load, and could dis- 
regarded studies where the movement fine material was primary 
interest. 

should noted here that, the movement the coarse material 
primary importance, careful sampling and special equipment are necessary 
obtain accurate estimate the total load coarse material. This was 
evident the Enoree River studies, where was found that the coarse material 
suspension may appreciable percentage the total load coarse 
material, even though small comparison with the load fine material. 
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Fic. RELATION BETWEEN DISCHARGE AND THE THAT 
PRECEDES THE DISCHARGE 


Fig. illustrates the variation concentration during typical isolated 
rise, and shows the general conditions under which suspended-matter content 
increases very rapidly from the beginning the rise until about the peak stage. 
The maximum content either occurs somewhat advance of, coincides 
with, the maximum discharge the stream, depending upon the stream dis- 
charge (Fig. 2). The intervening period time, between the occurrence 
the peak concentration and the occurrence the peak discharge, termed 
the 
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After the peak concentration has been reached, the content suspended 
matter decreases with the falling stage the stream, but the decrease usually 
less rapid rate than the increase during the rising stage. general, the 
recession curve concentration has long evidenced the muddy 
appearance stream for considerable time after discharge has returned 
approximately normal flow. the discharge and suspended-matter concen- 
tration continue decrease, point reached where the sediment load becomes 
negligible comparison with that transported during the rise. small 
stream this condition usually occurs only relatively few hours after the crest 
has passed. This illustrated the load graph Fig. which shows that 
only negligible load was transported after 3:00 a.m. April 1936. 

distribution graph calculated follows: For particular stream rise 
which sufficient number water samples have been taken define the 
variation concentration throughout the rise, curves suspended-matter 
concentration, discharge, and suspended load are plotted single diagram 
asin Fig.1. segregate the concentration due wholly the rise question, 
horizontal base line drawn ordinate, which appears, inspection, 
the value the concentration prior the rise. (The use recession 
curve, stream-flow studies, adds little the accuracy the results.) 
Having separated the base concentration from the observed values, the unit 
graph concentration represented the ordinates between the base line 
and the actual graph concentration. The time base the concentration 
graph considered arbitrarily extend from the beginning rise (8:30 p.m., 
Fig. the time when the concentration has receded the point beyond 
which the load can considered negligible (3:00 a.m., Fig. 1). The base 
the unit graph divided arbitrarily into time intervals, and the average con- 
centration for each interval then determined. These values are totaled and 
the percentage the total occurring during each the periods computed. 
Time intervals min and were found suitable for representing accu- 
rately the shape the graphs for rises occurring the East Fork Deep River. 

Distribution graphs were prepared for all the isolated rises that occurred 
the East Fork Deep River during the four years record from April, 1934, 
June, 1938. After reducing these various graphs constant time interval, 
was found that the shape the graphs appeared function the period 
stream rise, and somewhat the intensity the rainfall. The period 
rise this case defined the time intervening between the beginning 
surface runoff and the occurrence the peak discharge. Fig. shows the dis- 
tribution for ten storms which the period rise was less. 

These various graphs were prepared from rises that occurred during all 
seasons the year, and, consequently, represent various conditions land use, 
vegetal-cover, soil-moisture, and rainfall characteristics. Rises with the 
highest runoff and soil loss generally occurred during April and September. 
September the losses appeared the result torrential rains, whereas 
April high initial soil-moisture content was the most important factor. Two 
rises occurred during July result high-intensity rains, but the runoff and 
subsequent erosion were not severe for April and September, probably 
because the dry and porous condition the soil, but perhaps also because 
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the total rainfall was less and the evaporation and transpiration probably were 
greater. 

The average distribution graph the superimposed graphs Fig. 
tabulated and plotted Fig. together with the average distribution graphs 


GRAPHS SUSPENDED-MATTER CONCENTRATION FOR TEN 


for other periods rise, the percentages all cases being expressed for time 
interval hr. The graph for period rise less characteristic 
high-intensity rains short duration, when rain impact probably respon- 
sible for placing most the sediment into suspension. The graphs for the 
longer periods rise are characteristic long duration, low-intensity rains, 
where rain impact probably secondary importance placing sediment 
suspension. 

Rises resulting from uniformly distributed rainfall should used, where 
possible, preparing set distribution graphs. Although rainfall records 
are not necessary the application the graphs, such information 
importance judging the relative reliability the computations affected 
rainfall distribution. 
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One-hour intervals 
1.0 4.7 1.5 2.0 0.5 0.4 
9.2 8.7 8.0 4.5 2.4 3.1 
10.0 15.0 8.9 5.0 
11.5 10.0 14.5 9.2 
7.1 7.4 9.2 9.3 9.5 7.7 
5.7 6.4 7.5 8.4 
4.8 5.3 5.6 6.0 6.7 
4.0 4.6 4.0 5.0 6.2 6.1 
3.4 3.7 2.8 4.1 5.4 5.6 
3.1 1.8 3.4 4.7 5.1 
1.0 2.8 4.0 4.4 
1.9 2.4 0.5 2.3 3.4 3.9 
1.6 2.0 2.9 3.5 
1.4 1.9 1.7 2.6 3.2 
1.2 1.7 1.4 2.2 2.8 
1.5 1.2 1.8 2.5 
0.9 1.4 1.0 1.6 2.3 
0.8 1.3 0.9 1.4 3.1 
1.0 1.4 
0.9 1.2 
1.0 


Peak percentages are shown parentheses. 
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APPLICATION THE DISTRIBUTION GRAPH 

illustrate the application the distribution graph suspended-matter 
concentration the calculation suspended load, data were used from rise 
that occurred the East Fork Deep River November 
shows the observed discharge and suspended-matter concentration resulting 
from rain the characteristics illustrated the chart shown Fig. 5(a). 
The record from only one the recording rain gages shown. This 
tipping-bucket recording gage maintained the United States Weather 
Bureau the Greensboro, C., airport. 

The use the distribution graph, developed, requires that the stream 
hydrograph and one measurement suspended-matter concentration during 
rise available. the hydrograph composed several peaks, observed 
sediment concentration for each peak required. desirable that the 
sample for concentration taken at, shortly after, the peak discharge, be- 
cause, during the rising stage, the concentration often fluctuates erratically; 
whereas, during the falling stage, the rate change suspended-matter con- 
tent decreases uniformly, and less rapid rate. 

The hydrograph shown Fig. 5(b) consists two peaks, and therefore two 
observed concentrations are used applying the distribution graph. With 
known value concentration and selected distribution graph oriented with 
respect the peak discharge according Fig. computed distribution 
concentration can made for each peak. Where the effect two peaks 
overlaps, the individual curves are calculated separately, and then added 
graphically give the graph suspended-matter concentration for the com- 
plete hydrograph. The suspended load then easily calculated from the dis- 
charge and the calculated concentration. the load calculated for several 
points throughout rise, complete graph the load may constructed. The 
total tonnage sediment represented the area under the load 
desired, this load then may converted volume selection volume- 
weight ratio. 

the first rise, shown Fig. the period rise approximately hr. 
With peak discharge per sec, the time-of-lead zero according 
Fig. and the peak concentration therefore assumed coincide with the 
peak discharge. The 4-hr distribution graph, shown Col. the tabulation 
Fig. used and plotted Fig. 5(c), with the peak coinciding with the time 
peak discharge. happens that the observed concentration 860 ppm, 
selected for use the calculations, was taken the instant peak discharge; 
consequently, the concentration which the curves calculated and ob- 
served concentration will arbitrarily made coincide. multiplying this 
peak concentration 860 ppm the ratio the ordinate percentages for each 
time interval the distribution graph the peak percentage, calculated 
curve (curve Fig. 5(c)) concentration obtained for the first rise. may 
noted that the recession part this curve extends over period time 
sufficiently long contribute the total amount sediment suspension 
during the second rise; consequently, distribute the concentration the 
second rise, segregation the concentration resulting from the first rise will 
required. 


Oct 


Concentration 
IBeste nar Rainfall Depth 


Concentration 


Suspended Load 


7 


October, 1942 SUSPENDED-MATTER CONCENTRATION 


Concentration 
(Parts per Million) 


Concentration 
(Parts per Million) 


Suspended Load 


Rainfall Depth 


(Tons per Hour) 


1600 
LEGEND 
1400 Observed Concentration 
@ Observed Concentration.Used in 
1200 +120 Distributing the Calculated Concentration 
800 


N 


Calculated Concentration 
(2nd Rise) 


Load=91 


Noon A.M. 


GRAPHS SUSPENDED-MATTER CONCENTRATION 


1393 
Ww, \ 
Discharge 
800] Calculated Concentration 
Rise) 
= 
400}- 
a 


1394 SUSPENDED-MATTER CONCENTRATION Papers 


From the known period rise and the peak discharge, the 2-hr distribution 
graph (Col. Fig. selected and oriented (Fig. 5(c)), that the peak 
percentage coincides with the peak discharge, because here again Fig. shows 
zero time-of-lead for this discharge. The observed concentration used 
calculating the curve concentration this second rise segregated from the 
first rise subtracting the calculated concentration the first rise for the 
corresponding time; that is, concentration 390 ppm (read from the 
lated concentration curve for the first rise) subtracted from the total observed 
concentration 1,450 ppm give segregated concentration 1,060 ppm. 
The distribution-graph ordinate for this corresponding time 13.9%. 
multiplying the segregated concentration 1,060 ppm the ratio the distri- 
bution-graph percentages (Col. Fig. and the percentage 13.9, curve 
concentration (curve Fig. 5(c)) obtained for the second rise. adding 
graphically the individual curves for each rise (curves and Fig. 5(c)) 
calculated curve suspended-matter concentration obtained for the complete 
hydrograph. 

Although these calculated and observed curves concentration were made 
arbitrarily coincide two points, visual examination shows that agreement 
between the two can considered satisfactory. Close agreement further 
confirmed the load curves Fig. 5(d) where the calculated load for the flood 
period agrees within with the observed load. Note that the suspended 
load transported during periods low water negligible comparison with 
that transported during floods; hence, estimates load are made only for 
flood periods, sufficiently precise estimates will obtained for most practical 
problems. 

The primary advantage distribution graphs the fact that they permit 
estimate the suspended load stream made with only single obser- 
vation suspended-matter concentration during rise. This should 
particular value localities where sufficient funds are not available for ob- 
taining personnel for sampling throughout the flood period. Fewer samples for 
analysis also decrease the laboratory costs. This method, course, subject 
all limitations applicable the use the “‘unit hydrologic studies; 
that is, sufficient records particular location stream must available, 
some reasonable assumptions must made obtain reliable set distri- 
bution graphs. 

The application this method compound hydrograph the East Fork 
Deep River shown Fig. For this flood the time-of-lead, the period 
rise, and the observed concentration used distributing the concentration are 
shown for each rise. The actual recorded peak discharge was used orienting 
the position the distribution graph with respect the peak dischdrge the 
use instances, however, might more desirable use the 
peak discharge after segregation from the antecedent flow. general, however, 
appears that considerable leeway permissible the selection distri- 
bution graph and time-of-lead, without greatly influencing the estimate the 
suspended load. noted that the calculated and observed loads for 
the flood period July 29-31, 1936 (Fig. 6), agree within 10%, accuracy 
which, for most practical purposes, may considered satisfactory. 
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Items that may particular interest are the average infiltration capacity 
and the ratio surface runoff rainfall for the stream rise, shown Fig. 
the graph rainfall indicated the portion the rainfall that was assumed 
contribute the segregated surface runoff used the calculation both 
infiltration capacity, and the ratio surface runoff rainfall. The average 
infiltration capacity for the watershed has been calculated the method de- 
age infiltration capacity for several other flood periods was made and, would 
expected, showed rather wide range values which varied not only with 
season but also with antecedent rainfall and condition the soil. 

Plot studies have shown that, for constant slope, the main factors governing 
erosion are rainfall intensity, rainfall duration, initial moisture content the 
soil, and the effect the condition the soil surface. The most important 
factor affecting runoff and soil erosion the rainfall intensity, with the soil loss 
being affected greater degree. Describing the phenomenon briefly (7), 
when rain falls dry soil, the soil slaked and thrown into suspension. 
the rainfall intensity high enough cause runoff within few minutes after 
the rain starts, this suspended material carried off, and upon entering stream 
becomes the major portion the wash load. Rain falling moist soil only 
packs down and creates pavement effect that sheds the water, but erosion 
soil that originally was dry, most the material readily available for 
transportion suspension soon washed away, and the remaining wet soil 
packs down into smooth pavement-like surface. Although the runoff in- 
creases the rain continues, the soil losses decrease. 

That these general principles governing loss from small agricultural plots 
are also reflected the suspended-matter content drainage stream has been 
confirmed data obtained the East Fork Deep River July 29-31, 1936 
(Fig. first rise resulted from high-intensity rain that fell soil ina 
dry condition time the year when the soil surface probably was 
pulverized condition. relatively high suspended-matter concentration was 
attained result. Upon subsidence this initial rise, additional rainfall 
occurred, with the resulting runoff and suspended-matter concentration ap- 
parently being decidedly affected the first rise. result rain packing, 
and the flushing away most the loose material during this initial rise, the 
suspended-matter concentration during the subsequent rises was considerably 
reduced. fact, during the fourth and last rise which the highest peak 
discharge the period was attained, the maximum concentration was but little 
greater than during the lowest the peak discharges. This phenomenon 
apparent exhaustion the supply fine material, characterized the sus- 
pended-matter concentration reaching early peak, and then decreasing even 
though the discharge may continue increase, also has been observed the 
Enoree River investigations. Obviously, when the conditions governing the 
supply fine sediment stream are changing almost continuously with time, 
simple relationship between discharge and suspended load cannot 
exist. 
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The eventual solution the problem erosion will require considerable 
information the rainfall characteristics during each storm. reporting the 
results 1940 studies raindrops and erosion, Laws (5) made the 
following statement: 


intensity and duration records themselves not constitute 
adequate description rain. adequately describe rain far 
erosion and infiltration are concerned, necessary addition know 
the size and velocity its drop. The relative significance the drop 
characteristics depends upon multitude factors, such cover, slope, 
and soil type, and has yet evaluated.” 


CONCLUSIONS 


using distribution graph suspended-matter concentration con- 
nection with the observed hydrograph and single observed suspended-matter 
concentration, the total suspended load fine material transported during 
particular stream rise can calculated with reasonable accuracy. Previous 
records are necessary, course, developing set distribution graphs for 
particular location. 

Although the shape the distribution graph depends the period 
rise and, some extent, the intensity rainfall, unless considerable basic 
data are available define the effect intensity adequately, average distri- 
bution graph for various periods the rise should used. 
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DISCUSSIONS 


STABILITY GRANULAR MATERIALS 


Discussion 


objects the use law for the material the elastic zones, and 
quotes test results support his objection. The type stress-strain rela- 
tion expressed Eq. obtained triaxial test when the second and third 
principal stresses are held constant and the first principal stress raised 
progressively until failure occurs. The elastic modulus obtained dividing 
the first principal stress the total strain variable, beginning finite 
value and descending toward zero the test proceeds. The finite value 
the beginning the test due the initial pressure which must placed 
when the test set up. the test performed starting all pressures 
zero and thereafter bringing them all together, that definite and 
stable relation maintained among them, different result will obtained. 
this case the modulus, obtained the manner described, will start zero 
and increase the first principal stress increases. 

third result can obtained holding the second and third principal 
stresses constant, bringing the first principal stress about halfway the failure 
point, and then determining the elastic modulus observing the strains ac- 
companying partial release and restoration the first principal stress. 
change density the specimen due different packing will also modify 
the results. Some test values are given Table 

view these differences would appear that closer approximation 
the actual stress-strain relation than the simple one used the authors 
could chosen only after careful study the material, the proposed method 
placing the dam, and the loads placed thereafter. 

connection with Figs. and Mr. Benscoter presents discussion that 
appears reasonable for granular mass with cohesion. The horizontal direct 
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stress point Fig. corresponding point Fig. although not 
principal stress, only slightly greater than fourth the vertical stress. 
Similarly, for point slightly below the surface, the direct stress parallel 
the downstream slope almost twice that perpendicular the slope, although 


(Grain Size, 0.04 0.60 mm) 


CasE FROM VOLUME CHANGES 2.—CoMPUTED FROM TRIAXIAL TEST; 


Test pressure Volume strain Elastic Deviator stress Axial strain Elastic modulus 


(Ib per in.) (in. per in.) per ft) (Ib per in.) (in. per in.) per ft) 
10 0.02388 181,000 118 0.0105 1,618,000 
25 0.02453 440,000 184 0.0215 1,232,000 
50 0.02541 850,000 223 0.0323 994,000 
100 0.02641 1,636,000 249 0.0435 824,000 


Computed dividing the stress the total The linear strain computed dividing 
strain Computed dividing the deviator stress the total axial strain after the start 
of the test. 


these stresses are not the principal ones. states, Fig. corresponds 
non-linear solutions, whereas Fig. was for linear ones. 

Good agreement was found between predicted failure the reentrant 
corner common regions VI, and VII and actual performance several 
models, when instability was caused raising percolation gradients region 
VII. These tests are discussed further herein under remarks Professor 
discussion. 

The writers agree with Mr. Benscoter that the application mathematical 
theory soil analysis problems behind similar practice other fields, and 
they look forward further development. They also agree that new methods 
for measuring pressures embankments are needed urgently. 

Professor Hrennikoff comments some length the uniqueness the 
solutions and concludes that the proposed solutions are unacceptable unless 
some additional proof their reasonableness can supplied. support 
his conclusion describes procedures which attempts show that 
wide variety stress states could installed the dam some artificial 
conditioning the soil mass. The writers object his demonstration the 
following grounds: The passive, plastic-state, stress distribution requires the 
presence shearing stresses horizontal and vertical planes, and therefore 
would impossible set such state using smooth wall the 
manner describes. This statement true also for the active, plastic-state, 
stress distribution and all the possible intervening elastic states long the 
boundary slope not horizontal. similar difficulty will encountered 
the wall made rough. this case, any the foregoing stress states 
could set because the shearing stresses required the vertical planes 
could now maintained; but would not possible eliminate the central 
gap without altering the shearing stresses because removal the walls 
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would leave discontinuity the shear stresses along the plane symmetry. 
This discontinuity represents major violation the requirement for equi- 
librium. elastic zone introduced, these difficulties can removed and 
the shear stress can made pass through zero the plane symmetry, 
should. 

These difficulties with Professor Hrennikoff’s demonstration are considered 
sufficiently serious warrant reconsideration the adverse criticisms 
his section the latter part this section objects the use plastic 
states zones III and Fig. and expresses preference for elastic 
state throughout for slopes low those shown. reply this objection 
can stated that solution satisfying the elastic equations, and free from 
the aforementioned difficulties, will require the existence stresses along the 
boundaries and these cannot sustained cohesionless material. The 
active state was chosen for regions and VII for the examples Figs. and 
because was considered the state which probably would exist 
granular mass that had been deposited running water, would the case 
with many stream beds. 

Having disposed these details, consideration the question unique- 
ness may resumed. Professor Hrennikoff correct when states that the 
solutions are not unique. The writers would like urge, however, that the 
difficulty not due the form their proposed solutions but inherent 
the nature the problem. This may illustrated the retaining wall 
problem. Here the active and passive states represent limiting states beyond 
which stress distributions the linear type cannot exist. Between these 
limits, however, there are multitude elastic states possible, and should 
noted that all them, including the active and passive cases, are solutions 
the differential equations (Eqs. and 12) which meet the boundary condi- 
tions. This does not include all the possibilities, since additional solutions 
non-linear type probably can exist also. evident, then, that Rankine’s 
active state solution for the case surcharge, for example, not unique 
the sense that this distribution must exist the exclusion all others. 
remains physical possibility, however, spite the fact that tests and ex- 
perience show that greater horizontal pressures can exist. practice, other 


factors besides the satisfaction the differential equation and the boundary 


conditions must taken into account when decision must made 
the design pressure retaining wall. These factors will include considera- 
tion the manner placing the backfill and the conditions service after 
completion. 

necessary, also, exercise intelligent choice when selecting the 
stress states for the plastic active state for the downstream 
slope might reasonable for dam that constructed rolled-fill 
methods, whereas hydraulic fill dam might require different choice. These 
remarks will apply the comments Professor sections and 
The questions raised section regarding the use Hooke’s law, have 
already been dealt with connection with Professor Hennes’ discussion. 

regard the question raised Hrennikoff section con- 
cerning continuity stress, may stated that the choice here was matter 


4 
t 
n 
a 


1402 GLOVER AND CORNWELL GRANULAR MATERIALS Discussions 


judgment. the case dam subjected changing load, where the 
boundaries between elastic and plastic zones must move, was considered 
that discontinuity the normal stress parallel the boundary would lead 
difficulties regard compatibility strains the elastic side the 
boundary the case where this zone widening. other cases, choice 
equality for shear stresses along the boundary and for direct stresses normal 
zone VII are raised certain amount, the elastic zone will collapse 
simple boundary, subjected such stress condition. Example the 
percolation pressure portion incomplete, was stated the original de- 
scription. This due the limitations imposed the available solutions. 
this case, the choice percolation pressure was made that the actual 
body forces would represented well possible. determine how 
seriously the validity the solutions would affected these limitations, 
some tests were conducted model dams. These were built special 
tank which was about long, deep, and wide. The model dams 
were approximately high and rested foundation about deep. 
The case shown Fig. 6(c) was tested several occasions with dam whose 
upstream and downstream slopes were and respectively. The 
computations indicated that sliding failure should occur the downstream 
toe, with upward percolation pressure gradient zone VII, which would 
slightly less than that required cause piping. The exact figures for one 
the models tested were 113.25 and 123.10 per per cause sliding and 
piping, respectively. The sliding failure was observed take place the 
expected gradient and before any material had been removed piping. The 
estimated gradient cause piping was later confirmed very closely. was 
concluded from this test that the neglected portions the percolation pressure 
pattern did not influence the results the computations greatly. 

only necessary test one point inside elastic region the equality 
satisfied the two boundaries. The reason for this may stated fol- 
lows: The equation for second degree for any constant value 
(say when linear stresses and uniform friction are assumed. 
plot this equation will show intercepts the z-axis both 
and when plastic regions bound the elastic zone along and 
Therefore, check this inequality, the condition only 
one point along the line between and sufficient. The same 
reasoning would apply for any other constant value such 
therefore, check the inequality only one point within the entire region 
sufficient. 

Tables and Appendix were worked out restricting solutions 
Eq. real values that stresses given Eqs. would real. The 
discriminant was set equal zero, and values were determined for the 
various weights, coefficients friction, etc., shown the tables. The two 
systems percolation pressure contours used are shown Figs. 6(a) and 
Table was computed determining the percolation pressure gradient 
zone VII which would cause this zone assume the passive state. For some 
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what lesser gradients, zone absent and zone VII the elastic state, 
explained previously. 

Professor Krynine presents summary present methods analysis 
earth structures, which follows with brief geometrical interpretation 
conditions the plastic regions. This latter interpretation, together with his 
explanations some the mathematics, should prove helpful any one going 
through the method for the first time. His section showing that differential 
equations are not necessary plastic regions interesting presentation 
but would not shorten the computations for such regions. was postulated 
that the material the elastic regions obeys Hooke’s law, but this assumption 
need not made more reasonable behavior can stated. Whether 
law postulated, rectilinear stress the elastic region assumed, 
makes difference, course, when the designer limiting himself linear 
solutions, but stated the conclusion the paper general method 
approach was being developed. The procedure used the paper would 
hold even neither Hooke’s law nor rectilinear stress were assumed. 

closing, the writers wish take this opportunity express their sincere 
thanks those who have contributed the discussions. 
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DISCUSSIONS 


LATERAL STABILITY UNSYMMETRICAL 
I-BEAMS AND TRUSSES BENDING 


Discussion 


GEORGE WINTER, 


GEORGE WINTER,” connection with the general tendency to- 
ward higher working stresses and more slender members engineering struc- 
tures, problems the stability such members are becoming increasingly 
important and are given wide attention engineering research. Since the 
first publication the writer’s paper December, 1941, three articles re- 
lated topics have been published. Prof. has given com- 
prehensive and, many respects, fundamentally new treatment stability 
problems involving the simultaneous occurrence compression, bending, and 
twist. Mr. has investigated the buckling unsymmetrical I-beams 
pure bending—that is, the same problem which forms part the 
paper. Whereas the writer derived his equations means energy prin- 
ciples, Mr. Hill used differential equations for this purpose. His discussion 
the writer’s paper represents abstract his aforementioned 

For the case pure bending, Mr. Hill derives Eq. which, form, 
differs from Eq. 14a for the same case. Mr. Hill rightly points out, his 
equations are rigorously exact only the centroid the section midway 
between the flanges. The same holds true for the writer’s equation. This 
does not mean that these equations apply symmetrical beams 
though the areas the flanges must equal this case, the -ratios for the 
two flanges may different width, thickness). many instances 
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beams are designed structural practice, especially avoid local crippling 
the compression flange due large -ratios. 


The writer gratified find that, for the case equal flange areas (the 
only one for which both equations are strictly exact), Mr. Hill’s and his own 
expressions are identical. demonstrate this fact 
let Fig. represent I-beam with centroid 
flexural center and midpoint Then the 
distance introduced Mr. Hill, the distance 
from For equal flange areas the points 
and coincide. this case, with the notation 
the paper and Mr. Hill’s rule signs, 


which identical with the second term Eq. 14a. Likewise, substituting the 
value the fourth term Eq. 71, and combining with the third term, 
results in: 


2 


which, turn, identical with the corresponding term Eq. 14a. For this 
case, then, Mr. Hill’s and the writer’s expressions are the same. 

and not coincide, neither Mr. nor the equations 
are strictly exact. this case the difference between the two approaches 
that Mr. Hill implicitly takes account the distance from the centroid 
the midpoint However, neglects the quantity (see Eq. 

magnitude found from 


which and are, respectively, the areas the two flanges and 
the easily seen that the quantity the same order magni- 
tude Therefore, little gained accuracy taking account but 


neglecting Except for the rare case radically different values and 


A», the quantity seen rather small compared with orc. There- 


fore, the writer’s simple Eq. 14a represents close approximation even for the 
case which the flange areas are different. 
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The practical advantage the writer’s work-energy method that pro- 
vides rather simple way analyzing the stability beams under transverse 
loading rather than only pure bending. writer’s equations for pure 
bending and for center load are supplemented Lieutenant Schrader’s con- 
tribution those for uniformly distributed loading and for two equal loads 
symmetrically located the beam, derived the writer’s method. 
investigation these cases means differential equations would have met 
with practically insurmountable mathematical complexities. believed 
that the four types loading analyzed cover the most important practical cases. 

closing the author wishes draw attention the fact that the results 
his and other related investigations are directly applicable only the 
maximum stress due the critical load smaller than the proportional 
limit the material. This limitation, the same that the Euler formula 
for columns, common all calculations concerning the elastic stability 
structures. 
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TIMBER FRICTION PILE FOUNDATIONS 


Discussion 


Frank Am. Soc. was the author’s idea 
present, this paper, test data which seemed demonstrate some facts about 
group action timber friction piles that previously had been largely matter 
conjecture. The tentative theory for the calculation failure loads 
groups proposed the paper appeared check the test results and therefore 
seemed have some merit. Such theory can value only meets all 
conditions and improved criticism. The writer, therefore, sincerely 
appreciates the resulting comment and criticism. 

Mr. Connor offers number practical suggestions pile spacing that 
are too often neglected the designer. also gives data heavy loads 
actually carried friction piles driven into soft material which, during driving, 
did not appear able support such loads. Mr. Connor further presents 
interesting method comparing the strengths different groups piles 
the basis comparative perimeter areas soil around the groups. This 
method evidently gives close approximation between groups; but not 
possible use this method compare single piles with groups. 

Mr. Woodruff raises several points interest. The writer agrees that 
all predictions pile strength using this theory will not have the accuracy 
agreement between theory and tests apparently indicated Table nor 
does seem that such accuracy altogether necessary the design pile 
foundations. Further, quite possible, Mr. Woodruff suggests, that soil 
consolidation around the piles over long period time may affect their settle- 
ment somewhat; but does not seem likely that such long-time consolidation 
should affect the shear strength the soil adversely and lead pile failure. 
anticipated that settlement readings will made the Morganza Flood- 
way structures over period time. noted connection with 

paper Frank Masters, Am. Soc. E., was published November, 1941, 
Proceedings. Discussion this paper has appeared Proceedings, follows: February, 1942, 
Connor, Am. Soc. E.; April, 1942, Messrs. Glenn Woodruff, Jacob Feld, Greulich, 


and Paxson; May, 1942, Raymond Mindlin, Assoc. Am. Soc. E.; June, 1942, John 
Solomon, Esq.; and September, 1942, Cummings, Am. Soc. 


Cons. Engr. (Modjeski Masters), Harrisburg, Pa. 
Received the Secretary September 17, 1942. 


1407 


| 
| 
| 


1408 MASTERS FRICTION PILE FOUNDATIONS Discussions 


long-time loadings that one the pile groups the Crossing, 
group was loaded over period 243 days with constant load 
20.3 tons per pile without apparent additional settlement after the first day. 

Mr. Feld suggests that tapered timber pile actually may carry with it, 
failing, complete cylinder soil having diameter equal the butt 
diameter the timber pile. points out that under this assumption the 
shearing value the soil failure would 600 per and would 
more nearly correspondence with the value 558 per indicated 
Table for parallel-sided concrete piles. also describes approximate 
rule for determining group strengths subtracting progressive one sixteenths 
from the value all piles after the first, every new pile added the group. 
This apparently gives satisfactory results for the examples cited, but does not 
take account pile spacing relation length diameter, which the writer 
believes has important bearing the strength piles group. 

Mr. Greulich states that the law averages probably contributed the 
agreement between theory and tests recorded Table and with this the 
writer agrees. Mr. Greulich cites interesting example group action 
piles under the large foundations pivot piers swing bridges, where 
earlier designer found must restrict working loads low tons 
per pile. 

Mr. Paxson gives graphic comparison pressures beneath pile founda- 
tions compared with pressures beneath spread foundations, showing the 
considerably greater pressure under the spread foundation, and giving em- 
phasis well the fact that this greater pressure point near the surface 
where soil conditions are generally poorer than near the elevation pile tips. 

Professor Mindlin states that solution the equations elasticity, 
known the Boussinesq solution, was used the writer’s paper, and that the 
integrated Boussinesq formula not applicable pile foundations. not 
inferred that, because the writer used the Boussinesq solution, accepts 
the concept elasticity soils. The solution was used because 

‘appeared fit satisfactorily the results tests failure better than did 
any other formula. Nowhere the writer’s paper was there any mention 
the word this was deliberate omission because not believed 
that soils are perfectly elastic. This especially true connection with the 
tests given the paper which were carried failure, well beyond any possible 
range elasticity the soil. 

Professor Mindlin states that not possible compute curve Mind- 
lin distribution for uniform load per linear foot from the data given Table 
10. The method used the preparation Fig. was follows: Part (a) 
Table gives the pressures various radii from the pile for various depths. 
These data permit computation the total vertical force circle 0.1 
radius, 0.2 radius, 0.3 radius, 0.5 radius, such depths. For pile 
transmitting uniform load per linear foot the soil, the total load distributed 
distributed loads, there are subtracted the total vertical forces acting the 
circles various radii already cited, the difference each case the total shear 
carried the cylindrical surface the cylinder earth from the top pile 
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the depth 0.25 0.50 For each radius curve total shear may 
drawn, starting from the top calculated value 0.25 larger 
calculated value 0.50 The slope this curve any point the shear 
transmitted circumferential band one unit high. this slope divided 
will give the shearing intensity per square unit. The curve 
Mindlin distribution Fig. was based the slope the various curves 
midheight the pile. Any error that exists result the graphical 
method used believed slight. 

Professor Mindlin questions the principle superposition applied pile 
groups. course, there are many unusual stresses created between piles 
group due both the existence the other piles and soil stresses created 
during driving. Nevertheless, there seems reasonable way de- 
termining the build-up stresses under pile loads unless loads from various 
piles are superimposed upon one another. 

Professor Mindlin states that the real problem determine the distribu- 
tion along the pile the load transfer from pile soil. There are certain 
practical ways arriving this solution. When timber pile started the 
ground and the hammer placed upon it, will often settle into the soil 
distance several feet. the case the Morganza piles, this distance 
embedment into the soil under the weight the hammer used was such 
suggest that the shearing strength the soil near the surface was reasonably 
close the average shearing strength failure the was emphasized 
the paper, when some allowance made for the comparatively larger di- 
ameter the butt compared with the tip pile, and when there reason 
for the belief that the shearing strength soil increases slightly with depth, 
reasonable arrive such value was used the paper. In- 
cidentally, such practical data these, secured during driving operations, may 
often give clue the variations material from ground tip pile. 
recent study undertaken the writer’s office, dealing with pile foundation 
another locality, was discovered that this sinking the pile under the 
weight 9,700-lb hammer amounted ft, indicating shearing strength 
soil the surface the ground only 400 per ft, whereas the entire 
pile failed under average skin friction 975 per ft. This indicated 
considerable variation strength material and suggested method ar- 
riving value for the pile group question. 

Mr. Solomon demonstrates mathematically the effect the elasticity 
the pile upon and proves that for small pile loads the pile tip actually may 
deliver load the soil. Since the shearing stresses computed the paper 
were primarily related failure loads the piles, the value used must 
considered applying particularly such failure loads. Mr. Solomon also 
suggests that some the difference skin friction value between concrete and 
timber piles may due the effect weight the concrete piles compared 
with the weight timber piles. All the larger concrete piles were hollow 
and had weight greater than the weight the earth which they displaced. 
cylinder earth visualized standing soil, there will vertical 
shear along its vertical surfaces except that created applied load top 
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the cylinder. For that reason, the weights the piles themselves were 
uniformly and properly neglected computing skin frictions under loads. 

Mr. Cummings declares that the Boussinesq solution not applicable 
piles elastic soil, and shows that combination piles and soil not 
isotropic. restates the effect contraction the pile upon the distribution 
load the soil, and questions the equality loads the various piles 
group the Morganza tests. Mr. Cummings evidently feels that not enough 
measurements were made and that the Morganza tests were not the 
experimental standard established Galileo. that connection, although 
the writer has heard many suppositions the effect grouping upon pile 
strength, knows previous tests the scope the Morganza tests 
establish measurement what this grouping effect failure. The writer 
also wishes restate fact which appeared the paper—that these tests were 
originally intended satisfy the requirements current specifications, and 
that certain failures small groups suggested that the tests extended 
larger groups. The writer indebted the Corps Engineers, Army, for 
their cooperation making possible these supplementary failure tests large 
groups. The theory was then developed, using that method distribution 
which would check the tests failure and, such practical check, the 
method served cancel the numerous variables which are cited Mr. 
Cummings. 

Mr. Cummings also objects the use average value shear strength 
over the several crossings cases where single pile tests were not close hand. 
This was done because seemed the most reasonable approach, since funds 
had not been available make single pile test close each group test. 
However, after the paper was published, the Railway, its own 
expense, made failure test single 50-ft pile immediately adjacent the 
16-pile group, and obtained failure shear value 735 per 
this location. The 16-pile group used 60-ft piles, which might expected, 
due their increased depth, develop somewhere between this value 735 
and the average value 778 per used the study. 

conclusion, the writer hopes that opportunity may arise study the 
effect grouping the methods proposed, both the discussions and 
the paper, and urges that every effort should made relate all theories 
developed actual load tests the field. 
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DESIGN ST. GEORGES TIED ARCH SPAN 


Discussion 


the St. Georges tied arch span was made assuming the 3-ft arch rib 
entirely flexible, leaving the 9-ft tie girder resist the live-load bending 
moments without assistance from the rib. Under this assumption the action 
the structure exactly analogous that the self-anchored suspension 
bridge, inverted. The flexible arch rib, acting only direct compression, 
corresponds the cable the self-anchored suspension bridge, whereas the 
tie girder, acting both direct tension and bending, corresponds the 
stiffening girder the self-anchored suspension bridge, which must also trans- 
mit the horizontal component cable stress thrust over the entire length 
the bridge. 

The action the two structures being analogous, the equations and method 
analysis the self-anchored suspension bridge can applied with only 
slight modification the analysis the tied arch with flexible rib. One 
point difference the action the two structures found the 
effect vertical curve roadway and tie girder stiffening girder. The 
vertical curve the stiffening girder will decrease stresses the self-anchored 
suspension bridge, but the stresses the tied arch are increased the vertical 
curve. Since both structures are the type, resembling simple 
trusses that extent, the change dimensions produced deflection has 
important effect the stresses determined the elastic theory using the 
dead-load dimensions, the case either the self-anchored suspension bridge 
the tied arch. This distinction the effect deflection both the 
externally anchored suspension bridge and the arch bridge. 

more complex but similar structure consisting 3-span continuous 
girder with the center span reinforced flexible arch rib tied the girder 
analogous the self-anchored suspension bridge with unloaded side spans, 


paper Garrelts, Assoc. Am. Soc. E., was published December, 1941, 
Proceedings. Discussion this paper has appeared Proceedings, follows: February, 1942, 
bett, Am. Soc. E.; April, 1942, Jacob Karol, Esq.; May, 1942, Alexander Dodge, Esq.; June, 
1942, Carlos Bejarano, Jun. Am. Soc. E.; and September, 1942, Messrs. Eremin, and 
Freudenthal. 
Associate Engr., Robinson Steinman, Cons. Engrs., New York, 
Received the Secretary September 1942. 
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except that the case the tied arch bridge the side-span girders will not 
act direct stress. Again, however, the analysis the two structures could 
made application the same methods and equations. 

Professor Garrelts gives derivation for Eq. for the final analysis 
the St. Georges tied arch span. Although this may have been assumed 
obvious, the value this excellent paper may perhaps increased clarify- 
ing somewhat the theoretical basis for the final analysis. Hanger elongation 
under live load was considered negligible, were the tangential components 
deflection. Consequently, the angle between the rib and girder the end 
the rib, well elsewhere, not changed during deflection the 
ture, and, therefore, there can appreciable moment the rib its end 
connections with the girder. The span can made determinate cutting 
the rib the girder connections eliminate the horizontal compo- 
nent direct stress the rib and girder. The structure then will consist 
two simple beams, the arched rib, and the girder, which are connected the 
inextensible hangers and which must participate proportion their effective 
moments inertia carrying the live load. 

With only the value the horizontal component direct stress 
determined, the structure singly indeterminate, for the preliminary 
design. The moment the structure for unit value this redundant 
y’, this moment also being divided between the rib and the girder 
proportion their effective moments inertia. The direct stress the 
rib for unit value sec assuming the hanger loads equal 
throughout, and the direct stress the girder course 
for then follows from the substitution the foregoing moments and direct 
stresses the conventional Eq. for single redundant. 

Final analysis the previously mentioned tied-arch continuous girder can 
made this same method. taking the basic determinate system 
the combination the continuous girder and the center span simple arched 
beam, there would remain only the value the horizontal component 
direct stress, determined Eq. The actual values this 
equation, however, could taken for the simple span condition, and therefore 
only for the center span. Similarly, the first term the denominator could 
which represents the moment due unit value 
the basic continuous structure, and represents the moment due 
the center span simple structure. This follows from the fact 
that computing the deflection indeterminate structure the method 
the dummy unit load, the value either M’, the moment for the basic 
determinate condition, the moment produced the unit load, may 
taken for statically determinate condition. applying Eq. obtain the 
final moments, M’, course, would taken for the basic continuous system. 
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DISCUSSIONS 


PROFILE CURVES FOR OPEN-CHANNEL FLOW 


Discussion 


Messrs. STEVENS, AND Boris BAKHMETEFF AND 
NICHOLAS FEODOROFF 


Am. Soc. the writer began reading Pro- 
fessor paper, thought: ‘‘Here’s just what we’ve been looking 
and, when had read far his general differential equation for the surface 
profile (Eq. 13), appeared innocent that thought: better turn 
real mathematicians loose some our other Then, turning 
the page, that innocent looking thing the Professor states not difficult 
loomed out the darkness like the iceberg did the Titanic 
and wrecked everything. suggested that the Professor sentenced 
prepare table values his K-factor, terms and yo, varying 
tenths feet ft. One wonders how large volume would make. 

However, does prove that, defining the friction factor the Manning 
relation (which purely empirical and dimensionally incorrect), the surface 
profiles, determined for unit width wide channel (which does not exist), 
contain points inflection the range where depths are less than either 
normal critical, which range the Manning formula does not apply! 
Professor Gunder did quite well with his disclaimers, but challenged 
beat that one. 

After all, appears though hydraulic engineers will just have struggle 
along with their old step-by-step method computing water-surface profiles. 


serious error that has remained unnoticed hydraulic literature 
for years has been disclosed Professor Gunder. correctly surmised 
the author, the shape the so-called surface curve, traditionally presented 
textbooks, dates back the past century, when varied flow profiles were 


paper Dwight Gunder, Esq., was published April, 1942, Proceedings. Discussion 
this paper has appeared Proceedings, follows: June, 1942, Carl Rohwer, Assoc. Am. Soc. E.; 
and September, 1942, Posey, Assoc. Am. Soc. 
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first classified considering the special case broad rectilinear channels, 
with the further simplifying premise constant average Chézy factor 
applied over the entire range depths. Under these special conditions, the 


for all the depths, 


outline the remains concave 


tending toward approaches the critical value. further cir- 
cumstance taken into account that within the customary reaches the 
happen very exhibiting (except near only slight 
curvature and, the whole, approximating straight line. This means that 
the changes curvature, disclosed Professor Gunder, were not con- 
spicuous and the observed appearances did not upset accepted notions, 
Therefore, the over-all concave M;-profile, initially deducted for con- 
ditions, was inadvertently extended and assumed apply all cases, without 
due investigation the special effects devolving from variable friction factor. 
For such unpardonable negligence, the senior writer, the author manual 
especially devoted varied flow, presents his humble apologies. should 
stated, nevertheless, that least some the tracings the 


the conventional surface profile diagrams are certain distance 
5 t 


from the bottom, indicating that the vicinity the deductions drawn 
from the general varied flow equation cannot apply. 

Professor Gunder refers the features disclosed him 
only proper qualify previous ideas plainly erroneous, and accept 
the author’s discoveries regular characteristics the varied flow curves, 
especially these properties seem have broader significance than that 
suggested the paper, and they will found obtain quite generally, 
irrespective the particular formula used for the C-factor. 

The following analysis complements Professor pioneer findings, 
which from now should replace all former ideas relating 
Also, his discoveries are interpreted the light experimental evidence. 

(1) The first deduction, drawn Professor Gunder from Eq. 13, that 
the surface curve for broad rectilinear channels with Manning 
theoretically becomes vertical matter fact, this property 
not any particular form channel any specific friction 
formula. Indeed, constitutes general characteristic which applies all 
instances. the general varied flow 


when the depth approaches zero, the unity factors both the numerator and 
denominator become small comparison the second members, and 


Open Channels,” Boris Bakhmeteff, McGraw-Hill Book Co., Inc., New York. 
Y., 1932, 76. 
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gradient 
which the critical For small depths, tends toward unity. 
Hence, 


which the dimensionless friction drag coefficient the average wall 
stress expression 


reason its physical essence, well theoretically the customary 
Bazin, Ganguillet-Kutter, Manning, etc., relations, the C-factor for tends 


Note that this feature does not depend the value the channel slope 
other words, all surface curves, both the and the class, tend 

(2) Professor Gunder infers, properties derived from the region 
represent purely theoretical conditions. The important practical effect, 
nevertheless, that the reach small depths the M;-curve must convex 


become zero, making the critical slope ana thus infinite. 


\ 


there inflection curvature, with the surface profile changing from con- 
vex concave. The paper conveys the impression that the eventual location 
particularly affected friction coefficients expressed through the 
Manning formula. fact, the paper suggests (line following Eq. 14) that 
the range type curves, the Chézy coefficient should deter- 
mined either the Bazin Ganguillet-Kutter relation and not the Man- 
ning relation.” The recommendation apparently actuated the desire 
avoid the noted 

One should not try escape facts. What really matters whether the 
property question actual not; and under all circumstances the selected 
formula should the one that known comply best with established 
reality. Unfortunately, very little known about resistances rapid (shoot- 
ing) motion. such region investigators use formulas and coefficients 
deduced from observations tranquil motion without any direct 
assurance that they apply equally supercritical V.) flow. More- 
over, for “slowly varied (parallel) they assume that the friction rate 


Open Channels,” Boris Bakhmeteff, McGraw-Hill Book Co., Inc., New York, 
N. Y., 1932, p. 48. 
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for certain discharge and depth substantially equal what the resistance 
would have been under uniform conditions. the flow divergent 
and the divergence times becomes quite pronounced. Analogous closed 
conduits, one may anticipate complementary loss connected with the restora- 
tion kinetic energy into potential energy. Again, there material 
available present that would evidence the magnitude such restoration 
losses open-channel flow. 

Accordingly, many respects the precise answers the questions raised 
the paper will have await evidence from future research. the mean- 
time, one forced seek approximate solutions recourse one another 
the existing friction formulas, and here that the writers cannot quite 
follow Professor conclusions, for their investigations seem show 
that the choice any relation for does not materially affect the position 
the inflection point. 

(3) The location the “inflection the where the 
curvature changes from convex concave, can expressed conveniently 
dimensionless terms the ratio 


Generally, the will influenced the shape the channel, the dis- 
charge, the slope and the wall roughness. For two-dimensional motion 
(broad, rectangular channel), the shaping agencies are reduced the para- 
metric ratio 


and the form the function 


which depicts the change the Chézy factor with the varying depth flow. 

apprehend the possible effect the different formulas for neces- 
sary make clear the physics that actually underlie the structure all 
empirical relations. all formulas, “surface roughness charac- 
terized proper numerical term, such the n-factor the Ganguillet- 
Kutter and the Manning expressions. Such friction term held 
invariant for all conduits involving particular type surface roughness. 
Each formula offers the friction term with mathematical expression deemed 
represent best the factual change with the conduit all instances 
increases with the latter, the implied physical meaning being that the friction 
term exemplifies certain roughness size, and that the increase with the 
hydraulic radius reflects the eventual reduction the friction drag with the 
decrease the relative roughness. fact, relative roughness appears the 
empirical relations the form ratio the friction term appropriate 
power the hydraulic radius. Thus, the Bazin and Kutter expressions are 
built around the ratio the friction term VR, whereas the Manning 


formula the shaping ratio general, the aim research 
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the past was evolve mathematical form between and which for 
certain would best average and fit the maze observational data 
assumed respond the same roughness type. essence, the purpose was 


essentially the same, answered more recently the semi- 
rational theory, which expressed the interrelation between and 


the relative roughness logarithmic terms. 


180 


160 


Focal Point, 


140 


for 


100 


Chézy Coefficient Feet per 


Values Hydraulic Radius Feet 


comparison the different suggestions given Fig. which shows 
group f(R)-curves computed the different formulas with the 
Ganguillet-Kutter term 0.013. should remembered that, the 
otherwise rather complex and obviously artificial Ganguillet-Kutter struc- 
ture, the value C,, for 1m, constitutes focal point 
which the Chézy factor assumed unaffected the slope and 


metric measures numerically equal For feet, the focal value 


3.28 the same focal point that conditions the Man- 


ning curve, the metric form for which 


1.811 
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the other hand, the logarithmic expression 


represents, foot units, the curve, customarily referred 
the Nikuradze formula, adapted pass through the Ganguillet-Kutter 
focal point for the roughness term (In the well-known expression for rough 


pipes, 


or, foot units, 
1.811 


fitted pass through the focal point making the friction term 
— 1.811. 

Fig. the C-values are plotted logarithmically. seen, the curvature 
the outlines becomes rather unpronounced, and one naturally tempted, 


Values 


for example, replace the complex Ganguillet-Kutter relation suitable 
straight line answering the exponential relation passing through 
the focal point C,, and adjusted, far the incline concerned, approxi- 
mate best the the desired region. This explains the 
and the significance the Manning formula, which fitted conform 


with 


200 us 
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More generally, the light Fig. one may state that the various rela- 
tions versus differ the relative steepness their incline. Therefore, 
exponential expression were used for approximating the f(R)- 
function that any other reach the curve, the respective formulas could 
said differ larger smaller value the exponent. The particular 
effect Eq. would the substitution, the numerator, the ratio 


(4) the particular instance broad rectangular channel 11a), 
one obtains, with such substitution, 


The possible influence friction, appraised different formulas, de- 


pends determining the value which Eq. would corre- 
c 


ject the parametric values One could proceed directly 


o 


using the third quantity Eq. 34; but the numerical computations prove 


easier one first finds the respective and then determines 
o 


c c 


2 
the second quantity Eq. 34, and making one obtains, after 


reduction, 


0.4, thus embracing all the possible slopes Fig. presented Fig. 


c 
particular, the values for the vertical which corresponds 
flow over horizontal bed, are found from Eq. 
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dx? 


seen, the versus outline rather flat, except when approaches 
unity—that is, when the slope becomes close the critical, and the normal 
depth tends become critical. Therefore, under ordinary circumstances, for 


which, when differentiated and solved for gives 


i 


Values 


0.6 yo, the simpler expression (Eq. 37) can conveniently used for quick 
appraisals. the other hand, for all possible values the exponent, the 
inflection depth remains within 0.4 0.6 the critical, showing remarkable 
steadiness its position. Thus, rule thumb, assuming inflection about 
half the critical depth, seems justifiable approximation. 

(5) The purpose now disclose more direct fashion the exact 
possible effect that some other friction formula. Also, view the 
the desirable know more about their eventual 
outline the regions adjoining the inflection point. The general procedure 
would cumbersome and ineffective, recourse must had 
examples. Calculations were performed for wide range 
assuming two-dimensional flow ft, and making successively 


fri 


(a) 
A 
re’ 
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ft, and Using Eq. the form 


(38) 
one determined the dimensionless ratio versus with the 


values for the Chézy factor computed successively from the Manning, Gan- 
guillet-Kutter, and the (Eq. 31) formulas, and with the 
friction term 0.013. (In applying the Ganguillet-Kutter formula, 
identical set values for corresponding 0.0004, was used all 
instances. attempt was made use the Bazin relation. The latter 
well adapted for open-channel flow higher roughness 0.020 
low roughness, seen from Fig. the much too flat.) 

Computations for ft, ft, and also were made for flow 
horizontal channel, using Eq. 116 form for the latter, which would 
correspond Kq. 38, 


Prandtl Karman 


0.2 0.3 0.4 0.5 0.6 0.7 


1. dy 


Curves for the case presented Fig. The outlines 
shown are typical the plottings obtained all instances. All the 


exhibit only slight change over comparatively wide middle reach, 


revealing the previously emphasized flatness the M;-curve. The curvature 


with once more the 
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derived from Eqs. and becomes pronounced only near and 
close the bottom, where the computations based the Bélanger premises 
parallel motion with negligible curvature cease applicable. Again the 
curves corresponding the different formulas vary only detail, without 
changing the major features the picture. 

(6) behooves one now compare the results the theoretical computa- 
tions with actual observations. practice, surface curves occur the 
foot spillways and chutes, where the steep incline changes mild slope 
channel, the point where water flows from beneath sluice into hori- 
zontal mildly sloped bed. either case the initial surface, the transition 
reach, ipso facto concave. Accordingly, the initial depth less than y;, 
the subsequent flow must exhibit portion convex surface between two 
concave reaches. Such succession forms practice highly improbable. 


(a) 
Sluice Opening 
per Width 
=0.512 
Critical Depth Level 


Sluice Opening 
q=0.961 per Depth Level 
per Width 


=0.306 
= 
Distance, Feet Flume Bed Slope 


Moreover, whenever, hydraulics, certain equation results two distinct 
respective branches reflect altogether different physical regimen, that 
one the branches unstable and reality presents wavy features. 
example the double-branched which features the height the jump 
sustained flow from under sluice, and which the branch for lower 
kinetic energies corresponds jumps the type. Thus, 


Open Boris Bakhmeteff, McGraw-Hill Book Co., Inc., New York 
Y., 1932, 248. 
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respect the M;-curves (Fig. 6), one could reasonably expect that the lower 
and theoretically convex part the profile will prove reality more 
less undulated. Special experiments the 6-in. tilting flume the Fluid 
Mechanics Laboratory Columbia University New York, Y., confirm 
such anticipation. The profiles shown Fig. are typical the results 
obtained for wider range sluice openings and initial heads. each in- 
stance, the surface profile the reach succeeding the initial 
section exhibits distinct undular outlines. The waves are not large, however, 
and furthermore their size decreases with the increase the depth, the 
curve evolves into the final, practically linear, form. Most characteristically, 
the effect waviness was particularly accentuated the bottom pressures. 
Apparently, due the high prevailing velocities, substantial centrifugal effects 
were exercised even when the curvature, caused the waves, was compara- 
tively small. 

The question arises: Why were not all these facts noticed and recorded 
former research? The answer may possibly found the general flatness 
the M;-curves, and especially the fact that profile measurements were 
taken mostly considerable intervals. Under such conditions, the relatively 
small deflections either one the other direction from the approximately 
rectilinear profile could mistaken easily for erratic experimental deviations. 
shown the curves Fig. the actual features the can 
revealed only closely spaced observations. would most interesting 
have experiments performed larger scale and with equipment suffi- 
cient capacity. 

practical computations, the relatively small waviness may disre- 
garded, and the convex surface which evolves from Eq. for the range 
small depths may considered profile. 
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DISCUSSIONS 


HYDRAULIC DESIGN DROP STRUCTURES 
FOR GULLY CONTROL 


Discussion 


basins for large dams have been published, but the results have not been 
generally applicable. This paper should particular interest the 
designer small structures, the costs which are not sufficient support 
independent model studies, and the authors are commended its 
presentation, and the studies that preceded it. 

believed that word caution necessary, however, with respect 
the general applicability the results. The writer does not share with the 
authors the optimism implied the following statements (see heading 
cation Drop-Structure Rules Other Types Hydraulic 


“The principles used developing the rules for the design rec- 
tangular drop structures for gully are sufficiently general for 
application many problems. There great difference between the 
methods that have been used insuring energy dissipation and pro- 
tecting the structure from undermining and those commonly used the 
design spillways for major engineering structures.”’ 


casual reading the foregoing statements might lead the reader assume 
that the criteria developed the paper were directly applicable spillways 
for major engineering structures, although quite probable that the authors 
did not mean imply this. 

There are least two important differences between the gully-control 
structures described the paper and spillways for large dams. The first 
these differences contained the statement the authors (see heading 
“Development Design Formulas from Experimental Data: Height 
Transverse End Sill’’) 


paper Morris, Jun. Am. Soc. E., and Johnson, Assoc. Am. 
E., was published January, 1942, Proceedings. Discussion this paper has appeared Proceedings, 
follows: April, 1942, John Hedberg, Assoc. Am. Soc. E.; May, 1942, Messrs. Standish 
Hall, and Christiansen; June, 1942, Walter Wilson, Assoc. Am. Soc. E.; and September, 
1942, Proceedings, Messrs. Christensen and Dwight Gunder, and Boris Bakhmeteff and Nicholas 
Feodoroff. 

Senior Hydr. Engr., TVA, Hydraulic Laboratory, Norris, Tenn. 
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that the loss kinetic energy incident the impingement 
the nappe the apron large quantity and that this loss increases 
rapidly with increasing fall height offset, appreciable extent, 
but not completely, the increased energy the higher 


This may true for the type structure tested, which the water falls 
directly horizontal apron, but not generally true for major structures. 
There, effort usually made prevent this type energy loss because 
the tremendous impact forces that would exist immediately the toe the 
structure. Instead striking horizontal floor, the flow usually turned 
horizontal direction the toe means vertical curve bucket 
having radius, depending the depth flow, from 100 ft. This 
gradual change direction involves less loss energy and means that 
greater must dissipated the stilling basin. This difference 
not the side safety. 

The second important difference the assumption made with respect 
tailwater elevations. The tests apparently were made the assumption 
that flow below the structure would occur relatively high velocities and small 
depths. Thus, the ground roller below the end sill was always formed under 
similar conditions. Dams large rivers are usually relatively low, and the 
elevation tailwater becomes important factor. many cases, tailwater 
for maximum discharge may above the spillway crest. This true 
Dams Nos. D1, D2, and D3, Table and two other main dams the 


TABLE 5.—Comparison DIMENSIONS WITH Eas. AND 


(Fr) (Fr) 
Key designation 


= 
ee 


Tennessee River, and doubtless many other instances. When the tailwater 
sufficiently high above the spillway crest, the type flow changes and, 
instead plunging beneath the tailwater the apron the river bed, con- 
tinues near the surface for time. The conditions for energy dissipation are 
thus completely changed and the criteria proposed are longer applicable. 
The criterion for the length structure based the assumption that 
the discharge falls freely from the crest the apron. This not true for 
overflow spillways gravity dams because the requirements stability place 
the toe the dam downstream from the normal trajectory the overfalling 
sheet. result, the distance from crest end sill necessarily increased. 
and adopted the basis model tests for number dams built the 
Tennessee Valley Authority. These dams were all built with horizontal aprons 
and plain end provide typical examples both low-head and high- 
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head dams which general relationship might expected apply. The 
lack satisfactory agreement evident. 

Another factor importance many spillway designs the formation 
waves below the spillway and their effect bank erosion and navigation. The 
slope the upstream face the end sill has important bearing the pro- 
duction waves well the formation the ground roller. none 
the dams listed Table the upstream face the end sill vertical. Until 
many more experiments have been made under wide range conditions, 
special investigations should made for all major structures. 

The writer puzzled criterion (1), which states (see heading 
for Satisfactory Drop-Structure Performance’’): 


“The kinetic energy the effluent stream, measured the excess 
local velocity over the velocity for absolute minimum specific energy, 
shall minimum.”’ 


The reason for the choice the kinetic energy the critical velocity 
reference base not clear, since the total energy, and not the kinetic energy, 
that minimum that point. would also seem necessary specify the 
point which the energy (or velocity) was measured, since there consider- 
able dissipation energy short distance below the structure, illustrated 
Fig. the purpose the structure the prevention erosion, not 
lew bottom velocity the end desired? 

noted that extensive use was made photography for taking data. 
The possibilities photography for this purpose are very great but, far, 
comparatively little progress its use laboratory tool has been reported. 
The authors are commended for their use this versatile tool. 
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DISCUSSIONS 


EFFECT VARIATION ELASTIC 
CHARACTERISTICS STATIC 
UNKNOWNS 


Discussion 


effect the variation elastic characteristics and their application specific 
structures are well presented this excellent paper. the analysis given 
the author, statically determinate basic structure taken, and the redun- 
dant stresses are superimposed upon it. If, instead, statically indeterminate 
basic structure used (one having redundants, such used ob- 
taining influence line for the nth redundant), the computations are simplified 
and some new and interesting relationships appear. effect, this reduces all 
structures singly redundant types, assuming that the remaining redundants 
have been determined. avoid introducing new and possibly confusing 
nomenclature, the symbols used Professor paper will adhered 
to, but with the understanding that the moments, refer moments due 
the given loading the 1)-degree statically indeterminate base structure, 
and the moments are those the same base structure due unit load 
applied place the nth redundant. Thus, the only difference from singly 
redundant structure that the moments and must calculated for 
1)-degree statically indeterminate structure instead being determined 
alone. 


The analysis first will made for structures subjected bending alone. 
Then, general: 


— baa ( a. baa dbaa ) (20a) 


paper Hrennikoff, Assoc. Am. Soc. E., was published January, 1942, 
Proceedings. Discussion this paper has appeared Proceedings, follows: April, 1942, 
Finlay, Assoc. Am. Soc. E.; and May, 1942, Stefan Fraenkel, Jun. Am. Soc. 


(CEC), U.S.N.R., Electrical Eng. Dept., Naval Academy, Annapolis, Md.; and 
Assoc. Prof. Structural Eng., Univ. Minnesota, Minneapolis, Minn. (on leave). 
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and 
Eq. also can written: 


which the final true moment any segment, ds, and equal 
that dX, will have zero value (that is, change 


-grating along the members the structure, whenever either passes 


through zero the integrand changes sign; hence, the segments between these null 
points are alternately increased and decreased EJ. and should 
happen have common null point, this null point omitted, since both 
change sign passing through it. Thus, general rule could formulated 


follows: 


desired find the maximum variation dX, the moment 
point structure, determine the actual moments all points. Then 
consider the redundant replaced pair unit couples acting 
the same sense, and determine the moments the otherwise unloaded 
statically indeterminate base structure. Determine the null points; then 
compute dX, from 45, integrating between singular null points. 


N 


uo 


(b) CARRY-OVER AND 
(a) MOMENT DIAGRAM DISTRIBUTION COEFFICIENTS 


This analysis can readily extended structures with axially stressed mem- 
bers, indicated Professor Hrennikoff, that Eq. becomes: 


obvious that members having and like sign must increased 
area, and the other members must decreased area (or obtain the 
greatest increase, 


4 
0.8727 


em- 
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illustrate the application the foregoing principle, the frame Fig. 
will analyzed. The writer has used precise moment-distribution method® 
the analysis, gives accurate results with rapidity, but any method could 
used obtain the same results. Fig. shows the actual moments 
ends all members, and the moment diagram. The null points determined 
this method are shown heavy dots. The moment point expressed 


12 1-2? 18 


Fig. shows the analysis for Here the moment point taken 
the redundant and the frame analyzed for unit couples applied shown. 


2 12 


> 


(a) (b) 


(d) 
Fie. 


Fig. 11(a) gives the first moment distribution, assuming sidesway; Fig. 
gives the precise carry-over and distribution coefficients; Fig. 11(c) gives mo- 
ments due assumed sidesway; and Fig. 11(d) gives moments corrected for 


Precise Moment Distribution Method,” Joseph Wise, Journal, A.C.I., November, 1938, 
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Fig. locates another null point. These null points agree 
those determined the paper. 
the equations for and are written for the three parts the frame: 


Left leg— 


= 9 
2 2 
Beam— 
2 
Right leg— 
2 2 


and 


Substituting Eq. 45, 


2252 


f EI 2+ 215 dy 
| 
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This agrees exactly with Eq. 30a. However, the null points have already been 
found, and only necessary evaluate Eq. between these null points. 
The result will the same the paper. 

From study Eq. will seen that dX, will most affected 
changes where the product maximum and, obtaining the 
maximum each the three terms Eq. 52, will found that the most 


sensitive points the frame are: Left leg, 0.386; and right 


leg, Thus, curve could reveal the parts the structure 


beam 


which changes have the greatest effect the demolition 
structures might quite useful have this information. Also, the 
plotted, planimeter could used obtain the mechanical 
integration Eq. 45. 

The approximate integration finite summation can easily applied 
this analysis. The frame Fig. analyzed, the paper, substituting 
the numerical values from Table 


5,004 
1,201 


comparing the moment the crown and the moment the knee this 
frame, noted that the two are the same and, since obvi- 
ously the same for both cases, 


the variation dX, will TABLE 2.—ALTERNATE SOLUTION FOR 

obtaining 

Therefore, the K-values for 180 
inversely proportional 5.004 1.201 


the respective moments, 
248.64. The respective values for K—0.263 and 1.024—given 
the paper agree with this ratio. From this can concluded that tends 
largest where stresses are smallest. fact, one were investigate point 
inflection the structure, where zero, one would find becoming 
infinite, obviously should, since any change elsewhere the struc- 
ture would result moment this point where moment had existed 
previously. Thus, appears that not entirely satisfactory measure 
the sensitivity the structure changes the elastic constants. The 
designer not much concerned about changes moment where the moments 
are small and the stresses are far below the allowable working stresses. 
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suggested that factor 


used measure sensitivity, which the actual maximum fiber 
stress the structure the location the redundant and the allowable 
working stress. points inflection would indeterminate, being the 
form and would not considered significant. 

The suggestion made Professor Hrennikoff that sensitivity appears 
increase the number static unknowns increase can explained least 
partly the fact that, general, the redundancy the structure increases, 
its stiffness increases, that the deflections due unit couples, decrease. 
Since this term appears the denominator Eq. 45, decrease will cause 

The opinions assertions contained herein are the private ones the 
writer and are not construed official reflecting the views the Navy 
Department the naval service large. 


Corrections for Transactions: January, 1942, Proceedings, page 85, change 
read “dH (max) 4.18 Eq. 37d change 
and respectively. See also erratum page 664 
April, 1942, Proceedings. 
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STRESS CONCENTRATIONS PLATES 
LOADED OVER SMALL AREAS 


Discussion 


Messrs. MACCOLL, AND TELLER AND 
EARL SUTHERLAND 


able contribution the practical solution certain problems which are 
great importance. undoubtedly will welcomed all engineers who have 
work with the subject. 

Briefly, the idea that expressing approximately such quantities 
bending moments, twisting moments, etc., resulting from the application 
moderately concentrated load elastic plate, sums small numbers 
terms, which severally have simple interpretations terms the location and 
distribution the load. 

the writer the names ‘‘place coefficients’’ and “area coefficients’’ seem 
unfortunately chosen, for inspection the basic formulas shows that 
they are not coefficients any very usual sense the word. Eq. 13, 


are twisting moments per unit load. Hence, the 


and 


—are bending moments per unit load. Likewise, Eq. 14, 


writer would prefer some such names and moments,”’ 

suggested the paper that the so-called area coefficients, and 
may determined more easily from the partial differential equations 
and 25, than from the explicit formulas, Eqs. and The writer doubts 
that this true under any very general conditions. Whatever simplification 
the partial differential equations have afforded the particular cases con- 


paper Westergaard, Am. Soc. E., was published April, 1942, Pro- 
ceedings. Discussion this paper has appeared Proceedings, follows: September, 1942, Messrs. 
George Rich, Dana Young, and 


Bell Telephone Laboratories, Inc., New York, 
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sidered the paper has resulted, thinks, from the fact that those cases are 
very simple and special, and are not all typical the general situation. 
Certainly, set explicit formulas, comparable Eqs. and 12, the goal 
any method for solving differential equations; and should considered 
fortunate that the present case the formulas are simple and easily 
interpreted. 


portion this very useful paper which deals with the action 
large slab elastic subgrade particular interest the highway 
engineer and this particular part the paper that this discussion will 
limited. 

his original analysis the load-stress relations concrete pavement 
Dean Westergaard assumed that the load was applied the surface 
the pavement over uniformly loaded small circular area the cases 
(1) corner loading and (2) load applied interior point some distance 
from the slab edges and over uniformly loaded semicircular area for the case 
load applied edge some distance from corner. For these con- 
ditions loading, formulas for the calculation the critical load stresses 
were presented. 

the paper under discussion, the early analysis extended permit the 
computation not only the critical load stresses but also other load stresses 
the vicinity the loaded area both for areas that are circular shape and 
for areas other than circular shape. the latter, the elliptical area, such 
that used Dean Westergaard the numerical example (Table 2), 
particular interest since shape which approximates quite closely that 
the area contact large pneumatic tire normal inflation and capacity load. 

number years ago the Public Roads Administration began study 
the structural behavior concrete pavement slabs the Arlington Experi- 
ment Farm, Virginia. part this investigation was experimental study 
the load-stress relations pavement slabs uniform thickness, being, 
effect, comparison the behavior predicted the Westergaard theory 
with the behavior actually observed pavement slabs ordinary size. 
Four reports this work have been and the fifth (and final) report, 
which concerns the Westergaard analysis most directly, nearing completion 
(August, 1942). The reports already published describe the scope the 
investigation and the general procedures and apparatus used. 

One the questions left unanswered the original analysis the author 
was that the effect the shape the loaded area the magnitude the 
Senior Engr. Tests, Public Roads Administration, Federal Works Agency, Washington, 

Highway Engr., Public Roads Administration, Federal Works Agency, Washington, 


38a Received the Secretary August 29, 1942. 


“Computation Stresses Concrete Westergaard, Proceedings, Fifth Annual 
Meeting the Highway Research Board, 1926; also Public Roads, Vol. No. April, 1926. 

Structural Design Concrete reported Teller and Earl Sutherland. 

Part Description the Public Roads, Vol. 16, No. October, 1935. 

Part Effects Variations Temperature and Moisture the Size, Shape and Stress 
Resistance Concrete Pavement Slabs,”’ ibid., No. November, 1935. 

Part Study Concrete Pavement Cross Sections,” ibid., No. 10, December, 1935. 

Part Study the Structural Action Several Types Transverse and Longitudinal Joint 
ibid., Vol. 17, Nos. and September and October, 1936. 
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critical stress from given wheel load. long was unanswered has 
been the subject considerable discussion. effort obtain some data 
this particular question, limited program load tests areas simulating 
tire contact areas was made connection with the aforementioned general 
Since these tests concern one the problems analyzed the 
present paper, thought that the experimental data may interest. 
The tests were made concrete pavement slab normal size and 
8-in. uniform thickness. This pavement was laid stabilized subgrade 
formed incorporating sand with the original subgrade material and com- 
pacting maximum density the optimum moisture content. Loads were 
applied bearing plates cut shape which simulated the area contact 
low pressure pneumatic tires. This shape shown Fig. 
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Theoretical Relation 
Data 
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Fic. STRESSES WITH THOSE STRAIN MEASUREMENT 


conventionalized average fairly large number actual tire impressions 
and, even when extended the larger sizes, seems conform rather well 
available tire imprints. these comparisons the impressions airplane tires 
recommended capacity load were used. 


the loading tests four bearing plates this conventionalized shape 
the sizes were used, follows: 


Length 
(minor axis, (major axis) Area, 
inches inches square inches 
6.86 12.00 
13.71 24.00 289 
20.57 36.00 649 


27.43 48.00 1,155 
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The loading apparatus was arranged distribute the applied load 
uniformly over the area contact the bearing plate. Strains were mea- 
sured the upper surface the pavement the center the bearing area 
along both the major and minor axes the area. The gage used was the 
recording gage described the first the published reports mentioned 

The modulus elasticity the concrete equals 4,800,000 per in., 
and the modulus subgrade reaction equals 250 per in., determined 
deflection tests the pavement slab the manner suggested Dean 
Westergaard one his earlier 

Fig. 3(b) shows comparison between the observed stresses the afore- 
mentioned 8-in. pavement slab and corresponding stresses computed the 
methods described the paper. The comparison made for both the major 
and minor axes for each the four sizes loaded area that were used the 
experimental work. The wheel-load 20,000 lb; ratio 0.15; 
and the corrected combined K)’ was used. 

will noted that, general, there remarkable agreement between 
the theoretical and experimental values. The relative differences observed 
between the stresses along the major and minor axes given area are almost 
exactly those which the theory indicates should exist. The magnitudes 
the stresses for the conditions selected for comparison are close agreement 
except for the smallest bearing area where the values determined measure- 
ment are about 10% greater than the corresponding computed values. This 
affects the shape the curves showing the relation between the size loaded 
area and stress evident Fig. 3(b). However, even with this exception 
included, the observed effect size loaded area stress agrees very well 
with that expressed the theoretical treatment. 

Another comparison which interest concerns the numerical example 
given Table the paper. For the conditions assumed this example, 
indicated that dual tire with load 5,000 per tire (or wheel load 
10,000 will cause maximum pavement stress which only 1.5 times 
that caused single tire with load 5,000 indicated that 
using dual tires suitable size the wheel load can doubled with only 
50% increase the corresponding critical pavement stress. 

Another interesting point shown the computed stress values Table 
that, for the case the dual tire this comparison, the critical stress 
has practically constant value from the center the area contact one 
tire the center the area contact the other tire. 

Several years ago the Public Roads Administration made limited series 
load tests which effort was made compare the stress produced 
given load when the load was applied through single circular bearing plate 
with the stress developed when twice that load was applied through two plates 
the same size, separated distance comparable the dual tire spacings 
large tires. Although the program was quite limited and the data apply 


“Analytical Tools for Judging Results Structural Tests Concrete 
gaard, Public Roads, Vol. 14, No. 10, December, 1933. 
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only circular areas contact, the tests afford interesting comparison 
with the indications Table 

the tests referred to, loads were applied 8-in. circular plate the 
interior region slab uniform thickness. Strains were measured the 
direction the longitudinal axis the pavement the center the loaded 
area and various distances from the center line perpendicular the 
longitudinal axis the pavement through the center the loaded area. From 
these strain data curve was developed showing how the magnitude the 
longitudinal stress varied across the pavement slab the distance from the 
loaded area was changed. 

With this stress-variation curve basis, the effect two loads was 
examined combining two such curves displaced horizontally in. with 
respect each other shown Fig. and obtaining combined stress 


P=8000 


Circular Plates, 
8-In. Diameter 


Ratio: =1.45 


(a) STRESS DIAGRAM DEVELOPED 
FROM MEASURED STRAINS 


Stress, Pounds per Square Inch 


values addition. This study indicated that load 8,000 lb, applied 
through dual 8-in. diameter areas, would cause stress 1.45 times great 
would caused 4,000 applied through single area the same size. 

Fig. shows some additional data obtained the same experiments 
tests which strains were measured under and between the loaded areas. 
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will noted that both the combined stress curve Fig. 4(a) and the data 
for the 11-in. spacing Fig. 4(b) indicate that the stress between the loaded 
areas practically equal that within the loaded areas. 

Although the data presented this discussion are quite limited scope 
and apply but one part the paper under discussion, the close agreement 
between the behavior predicted the theory and that observed during the 
loading tests tends inspire confidence the general soundness the the- 
oretical treatment. 


— 
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RELATIVE ANGULAR, LINEAR, AND 
TRAVERSE ACCURACIES 
CITY SURVEYS 


Discussion 


PHILIP KISSAM, AND GEORGE WHITMORE, 


CLEVELAND Am. Soc. author this paper joins 
small and select group—those who offer something for the advancement 
the professional knowledge and ability the surveyor. any issue 
almost any technical publication, structural, hydraulic, and sanitary engineers 
may find results investigation and analysis which will bear more less upon 
their daily problems, but the land surveyor the the pro- 
fession, although not the general public, whom the 

Truly, the phrase used Mr. Bauer masterly one (see 
“Under the conditions and with the equipment that the economics the pro- 
fession force upon The economics the profession have forced con- 
siderable shoddy practice many men who knew what should done give 
reliable results their clients. The clients neither know nor care how the 
surveyor arrives his results; nor are they interested the accuracy attained, 
many cases. 

The writer surveyed farm lands for many years, competition with day 
labor the crudest type. Very frequently went job alone, and assis- 
tance was provided the client, the form hissonsand nephews. This was, 
and is, condition very common occurrence and will continue long 
registers deeds are permitted law accept the descriptions and plats 
made unlicensed surveyors. Much has been done correct this condition, 
but much educational work remains before the landowner becomes willing 
pay the price which efficient help and accurate measurement land should 
reasonably cost. 

The paper deals entirely with city surveys and concerns itself with property 
considerable value. The writer’s experience was with mountain and farm 


paper Bauer, Assoc. Am. Soc. E., was published May, 1942, Proceedings. 
Lt.-Col., Corps Engrs. Army, Knoxville, Tenn. 
Received the Secretary July 17, 1942. 
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land, worth from $10 $300 per acre, and the conditions and equipment forced 
upon him the economics the profession were certainly inadequate, even 
perhaps someone will present the rural point view. 

The survey farm lands the private surveyor, distinguished from 
cadastral surveys large areas government corporation surveyors, 
generally one the following categories: 


(1) establish dividing line disputed corner; 

(2) ascertain true acreage; 

(3) subdivide property for inheritance sales purposes, either large 
tracts small lots. 


Case (1), the experience, comes more under the head art 
Tennessee, there are applications the law township 
and section, but all property traces back original grants, mixed the 
wildest confusion. case has the procedure been clear-cut and capable 
solution entirely scientific methods and careful application angular and 
linear measurement. Diplomacy important tool the surveyor, and 
practically all final results are the result compromise, either known 
unknown the parties the dispute. Under such conditions, repeating 
angles and tension tape play part. 

Case (2) requires closer adherence the rules and the one which old 
notes can reviewed and rerun. Unfortunately for such studies, very 
seldom that corner monument can exactly occupied. Nearly all lines are 
lines, with temporary monuments, tied the corners side shots. 
Here the parsimony the owner plays great part, not only preventing the 
employment skilled help but making necessary run the offsets 
places requiring the least clearing, not following the boundary closely 
should be. Here, also, the rural objection measurement horizontal plane 
manifests itself, not infrequently making necessary measure vertical angles 
any means hand and correct slope measurements the office without 
the knowledge the client. 

Case (3) extension case (2), with the boundary survey checked 
numerous crosslines. this case, much closer adjustment possible. 

Unfortunately, the writer’s notes are stored the would 
interest for some country surveyor produce tabulation his work, 
showing the accuracies attained under fire (sometimes literally), set forth 
against the value the land surveyed. All credit due the many honest and 
efficient men all parts the United States who give reasonable results under 
very unfavorable conditions and thereby aid the upbuilding the profession, 
which has been sadly disgraced different times and places men who are 
neither efficient nor honest. 


procedures, and analyses the results obtained, are meager and gratifying 
have such report Mr. Bauer. 


Associate Prof., Civ. Eng., Univ. Illinois, Urbana, 
Received the Secretary July 29, 1942. 
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well exemplifies the fact that, proper attention the more important 
sources error, the accuracy surveys, using ordinary equipment, and 
small additional expense, can greatly increased over the crude results 
often obtained. 

There are very definite criteria for comparing the merits the different 
surveying instruments, and has become customary use the least count 
the vernier designate the character the instrument used. This unfortu- 
nate because there are wide differences between transits having the same vernier 
count; also, many transits having verniers are superior others with 30” 
The better instruments yield better results with such saving 
time that poor economy use any but the best quality surveys any 
considerable extent. 

The use the for two observations computed Eqs. 
and and described under the heading Linear Readings,”’ 
justified. This quantity based the assumption that normal 
distribution the accidental errors (residuals) series observations 
under investigation, and such distribution can obtained only where con- 
siderable number observations are used. certainly dubious practice 
apply the formulas series less than ten observations. 

25” (see heading ‘‘Findings Fact from Item (f)), but 
not clear under just what conditions this value was found. particular, 
important factor the reading lens that was used. The writer found his 
investigations‘ that, comparing the results reading angles with and with- 
out reading lens, the probable error single readings with the naked eye was 
reduced about one half when three-diameter lens was used. was also 
found that large part the error single readings angles caused the 
inability the observer set the vernier precisely zero. This condition 
explains, part, why the errors reading angles with ordinary equipment 
are greatly reduced when the method repetitions used. 

The analysis results evidently has not differentiated sufficiently between 
the effects accidental and systematic errors. For example, all traverses were 
run loops, wherein the effects systematic errors not become evident. 
Thus, tape were incorrect length, all sides would affected proportion- 
ately, and the calculated error closure would way affected this 
condition. Also, constant error the measurement the temperature 
would produce like effect. Further, the duplicate measurement line 
under the same conditions will not reveal the presence any systematic errors. 
This fact probably accounts for the fact that the average accuracy lines 
measured different years was less than that for lines measured forward and 
back the same day (see Fig. 3). 

Accidental errors are compensative their effects, and hence the errors 
expressed this report should smaller for large number observations 
than for few. This exemplified the fact that the accuracy found 
taping long lines was greater than for short lines (curves and Fig. 
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latter consideration not verified, however, the case traverse closures, 
for significant difference was found for traverses different lengths. This 
unusual result seems call for some explanation. very probable that 
the computed errors closure would larger the traverses had closed upon 
each other, instead being run closed loops. 

would seem that, where much care was used measuring the distances, 
reference would made the Bureau Standards determine the true 
lengths the tapes used. The fact that two tapes are the same length 
not complete evidence that they are correct. Their true length never known 
until they have been compared with the United States standard, with 
tape which has thus been compared. 

The effects the combined errors distance and angle measurements 
traversing are scarcely simple that indicated Eq. The writer 
attempted analysis® these effects 1928. The relations between the 
linear and angular errors are complex and, general, their effects must aug- 
ment each other, that the total error closure larger than the effect 
either one taken separately. 

Considering the foregoing remarks, the following comments may made 
about the conclusions offered the author: 


(a) The analysis the sixth paragraph preceding Eq. leads conclusion 
(2)—namely, that line, forward and back, under average condi- 
tions, adds nothing actual This analysis faulty that 
disregards the effects systematic errors that would appear measurements 
made different times but not repeated measurements made under the 
same conditions. other words, probable that the true errors the 
measurements made the same day are larger than the errors indicated 
the duplicate measurements. so, conclusion (2) cannot justified the 
evidence and analysis used. Clearly, regards the accidental errors involved, 
the error the mean two measurements will less than that for single 
measurement. 

(6) Conclusion (3) based the oversimplified analysis the relations 
the various errors traversing. conceivable that, given traverse, 
the effects angular errors will reduce the effects the linear errors, vice 
versa, but general, where the laws accidental errors are operating, the 
total error closure will comprised both angular and linear errors 
will greater than the error produced either class alone. 


Am. Soc. definitely, the paper Mr. 
Bauer contribution the art surveying. the author states, too 
little known the actual accuracies obtained the field the usual property 
just what can done obtain the desired accuracies the 
minimum cost. more surveyors would follow Mr. lead compiling 


“Specifications for Transit Traversing and Stadia Leveling,’’ Rayner, Transactions, 
Soc. E., Vol. (1930), 679. 


Associate Prof., Civ. Eng., Princeton Univ., Princeton, 
Received the Secretary August 1942. 
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data has, there would question whatever that accuracies, methods, 
and instruments would improved and costs reduced. 

Errors Transit analyze Mr. conclusions, the writer 
has classified the errors the transit follows: 


Major errors scale reading— 

(1) Errors determining what the scale reads (probable error r). 
Minor errors scale reading— 

(2) Errors graduation. 

(3) Errors due play between the inner and outer spindles. 
Major errors 

(4) Errors due lateral refraction. 

(5) Errors due the positions the instrument and signals (including 
phase). The probable error Items and The two 
types cannot separated using the data given. 

(6) Errors due not pointing the cross hair the apparent signal 
(probable error p). 

Minor errors pointing— 

(7) Errors due temperature changes the instrument and tripod. 

(8) Errors due lack adjustment the instrument. 

(9) Errors due play the slide and the trunnions and errors 
optics. 

(10) Errors due lack verticality vertical axis setup. 
(11) Errors due play between the outer spindle and the leveling head. 


This classification depends the assumption that the transit good working 
order and that proper operational procedure used. 

this discussion only the major errors will considered the minor 
errors are negligible under the conditions stated. 

Mr. Bauer has determined the accuracy single reading transit 
tabulating 176 sets angles. each case the first reading and the mean 
six turns were recorded, together with the difference between the first reading 
and the mean six. The mean difference the entire 176 angles tabulated 
was found According the theory errors the probable error 
would be: 0.848 21” 

will seen that this error made follows: 

was made one time and with the same positions for the signals and the 
instrument. Since each pair observations was made the same time 
assumed that the lateral refraction was also constant; thus: 


which the actual error the mean six repetitions, and the ratio 
between the error setting the vernier zero and the error reading 
transit the actual error the first reading: 
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and 
The probable value Eq. is: 
Mr. Bauer finds the probable error set six turns 3.8; hence: 


Solving Eqs. and simultaneously: 


0.3 16.2 4.4 
0.5 15.5 4.3 
0.8 14.0 4.0 
1.0 13.0 3.9 


Assume the third set values correct 0.8). 

possible determine Mr. Bauer has determined the probable 
error angle measured six repetitions when the measurements were 
made different times. Such error includes The value determined 
fraction and the observer can assume that the error position the instru- 
ment one half that the signal, the angular error due setup 

Since the average length line 700 ft, the probable error position 
expressed follows: sin 700 0.0066 ft, which the error 
setting the transit. 

This computation results the following probable errors angle determi- 
nation: 


(reading scale) 14.0” 
(setting, zero) 11.0” 
(position transit and signals 700’) 5.0” 
(pointing apparent signal) 4,0” 

(transit position) 0.0066 


These depend two assumptions. The first that the vernier can set 
zero with 0.8 the error reading the scale. The second that the error 
the signal position twice the error the transit position. 

There are data Mr. Bauer’s paper from which can deduced the 
decrease the probable error that would result from reading the minute transit 
15” the effect finer graduations. the opinion the writer that: 


There considerable reduction probable error when the vernier 
transit read This procedure has been followed since 1933 
the New Jersey Geodetic Control Survey with considerable reduction 
rejected angles. 
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The fineness graduation has little effect the accuracy, the accuracy 
being chiefly function the clarity the graduation presented the 


Accuracy Linear Bauer has given the probable errors 
linear measurement for immediate retaping and for retaping later date. 
states, this not based comparisons with precisely measured high- 
accuracy surveys. the accuracy taping for given field procedure de- 
pends chiefly tension, temperature, slope, and tape calibration, retaping 
the same party the same day with the same spring balance gives poor 
determination the accuracy, and, Mr. Bauer points out, little practical 
value except pick blunders. The probable error for this condition 

From the theory errors (which Mr. Bauer quotes) the following evident: 

error 
Probable error 0.954 Hence, the total error accuracy for 
linear measurement 10,880. 

Traverse would seem that when traverses are measured the 
field procedure should such that the probable error the transit would 
about equal that the linear measurements, Mr. Bauer suggests. the 
angles were measured once with the telescope direct and the measurement 
repeated once with the telescope reversed (indicated herein DR); and 
both verniers were read each time, the expression for the probable error 
angular measurement would Substituting the 
values computed, the probable error 

Introducing the errors position signals and transit, the resulting 


accuracy This probable error corresponds very 
closely with the probable error linear measurement found Mr. Bauer, 
which the equivalent angular accuracy Thus, where average 
accuracy 10,000 for traverse required, the angles should measured 
DR, both verniers should read, and the usual methods distance measure- 
ment used. Such traverses would have accuracy between 5,000 and 
15,000. 

higher accuracies are desired, the method taping should improved. 
probable error 52,000 would necessary. The use the plumb bob 
would have eliminated. Probably this could obtained commercial 
work using 300-ft tape fully supported some surface, like sidewalk 
street pavement, using thermometer and spring balance, and running rough 
profile over the tape for slope corrections. When such procedure used the 
accuracy measurement will usually attain better than 25,000, having 
probable error 52,000. 

probable error the transit work would require considerable improve- 
ment the centering signals and possibly the transit itself. 


Transit Work Depends Clearer Circle Graduations,” Philip Kissam, Engineering 
News-Record, April 10, 1941, 51. 
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The desired results could obtained the signals were designed have 
phase and the work could performed any time when there was least 
some haze cloudiness that the lateral refraction would reduced, and 
the signals and the transit were set that the probable error their positions 


was 0.005 ft. Under these conditions 2s; but, 700 ft, 

Transit accuracy also would have maintained that lower probable 
error would result than Mr. Bauer found the rule. 

obtain this accuracy suggested that one set used, reading 
both verniers the beginning and end the set. Then, the probable error 
operation and the error the positions the transit and the signals, the 
probable error angle determination would (2.8)? 4.2”. 
This practically low enough. 

This error would reduced, course, using lines more than 700 
that reasonably good balance between transit accuracy and taping accuracy 
would maintained. Assuming that transit errors and taping accuracy 
affect the traverse equally, the probable error such traverse would 
about 50,000, and the accuracy about 24,000. 

Assuming that Mr. traverses represent the usual 
traverses run good land surveyors, the accuracy expected between 

The transit work much higher degree accuracy than the taping 
and long the present method taping persists there advantage 
turning the angle more than DR, using both verniers. 

high-grade traverse work required, accuracy 24,000 can 
obtained with little, any, increase cost 300-ft tape used and the 
angles turned reading both verniers. 

These recommendations are based the assumption that better target 
used than ordinarily found practice and that the noon hours clear days 
avoided for angle measurement. 

will noted that determining the order accuracies traverses, 
Mr. Bauer has assumed that Manual No. uses error speci- 
fying the order accuracy. the belief the writer that “error 
meant the total error divided the total length, value which 
much more difficult obtain. For example, the average position closure 
required for first-order traverse 35,000. obtain this accuracy the 
average probable error would have 73,500. 


Soc. E., anp interesting approach the complex 
subject traverse accuracies contained this paper. The writers have long 


Chf. Surveys, Maps and Surveys Div., TVA, Chattanooga, Tenn. 

Civ. Engr., Maps and Surveys Div., TVA, Chattanooga, Tenn. 

Associate Mathematician, Maps and Surveys Div., TVA, Chattanooga, Tenn. 
10a Received the Secretary September 24, 1942. 
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felt that too much emphasis has been placed the closure 
the entire loop, and not enough designing measuring instruments and measur- 
ing procedures that will insure uniform precision for each and every tangent 
the entire traverse. Experience shows that traverses having circuit closures 
high order, such 10,000 20,000, may conceal errors, the indi- 
vidual courses, the order 1,000. writers believe that the accuracy 
specified for traverse loop line should indication also the approxi- 
mate precision each individual course the entire line. 

Traverse does not have many rigid geometric checks available does 
triangulation, and for that reason safeguards against blunders are imperative. 
Within the limits accuracy desired the traverse, known systematic errors 
should eliminated, and measuring procedures should used which will mini- 
mize the unknown systematic and accidental errors. Measuring instruments 
should used that will provide the required accuracies most economically. 

Unless consideration given some the systematic taping errors, third- 
order traverse accuracy can attained only under exceptional circumstances. 
For instance, any one the following conditions would result systematic 
errors that would equal the limits for third-order traverse accuracy 5,000): 


(a) the actual length 100-ft tape differed much 0.02 
from the true length, systematic error would result. Such 
change length may caused stretching kinking. 

(b) Under certain conditions systematic error the tension applied 
the tape would result systematic error 5,000. 

(c) the actual temperature the tape error 30° systematic 
error would result. Such error may result from applying 
temperature correction, or, under extreme conditions, from applying the 
temperature the air rather than the tape itself. 


There are, course, other systematic taping errors, the accumulated effect 
which may appreciable. 

The need for reducing systematic taping errors felt immediately when the 
traverse runs between two previously fixed points, such high-accuracy tri- 
angulation stations, part traverse net run under varying weather con- 
ditions two more parties using different measuring equipment and pro- 
cedures. These errors, obviously, will not ordinarily felt when dealing only 
with single traverse that run loop, originating and terminating the 
same point, such described the author. 

The attempt eliminate systematic errors, both linear and angular, will 
influence, large extent, the measuring procedures and instruments used; 
study surveying accuracies incomplete until the effects such errors are 
considered. Systematic angular errors fortunately can largely eliminated 
proper observing procedures and distributing angular closures before 
computing latitudes and departures. 

The writers believe that accidental angular errors have considerable 
effect the quality traverse, but that accidental taping errors, because 
their more nearly complete compensation, have much lesser effect and account 
for only small part the total traverse closure. The effect accidental 
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angular errors can reduced, however, means any all the following 
safeguards: (1) Frequent azimuth control lines through triangulation; (2) 
frequent astronomical observations; and (3) use long azimuth control lines, 
which span several the shorter, directly taped courses. 

The writers would like add the following thoughts the three questions 
posed the end the paper: 


The limiting the closure error not adequate assure desired 
traverse accuracy. The writers believe that the specifications should include 
the measuring procedures and instruments used, that each individual 
tangent will the required accuracy. 

Smaller plate graduations transits for city work are justified primarily 
for reasons economy, because the desired angle accuracy can obtained 
with fewer number turns. 

The writers find little criticize the design the newer transits 
now the market. course, improvements will continue made; 
however, they will probably the direction greater accuracies and more 
stability. 


AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


DISCUSSIONS 


MOMENT BALANCE: 
SELF-CHECKING ANALYSIS RIGIDLY 
JOINTED FRAMES 


Discussion 


Cornish’s original adaptation the slope-deflection formulas occurs con- 
nection with moment calculations for “ordinary defini- 
tion such frames appears available but, doubt, essential features are 
that the beams shall prismatic; that the inertia moment any haunch, con- 
sidered rectangular section, shall not greater than the assumed inertia 
moment the midspan tee section; and that sidesway due vertical loading 
negligible. 

The number moment calculations for complete multi-story building 
frame usually large and magnified the alternative arrangements 
superload which give maximum moments. Since the relative K-values are 
the first place unknown, common practice proceed means intelligent 
approximation these values, based the first place consideration fixed- 
end beam moments and statical reactions, afterward adjusted the simplest 
possible end moments for limited sections the frame. For 
these preliminary estimates moments doubtful whether any simpler 
method has been found, likely found, than single-cycle distribution 
and carry-over individual fixed-end moments for dead load, summation 
such moments and further summation the addition moments for alterna- 
tive superload arrangements, obtained proportion from the relation 
superload dead load. 

the final stage calculation, which follows the preliminary estimates, 
that the characteristics the method appear most 


paper Cornish, Esq., was published May, 1942, Proceedings. Discussion 
this paper has appeared Proceedings, follows: May, 1942, Messrs, Matthews, and William 
Larsen; June, 1942, Messrs. Eremin, Grinter, and Aleck; and September, 1942, 
Messrs. Frederick Merritt, and Ralph Stewart. 

Cheadle Hulme, Cheshire, England. 
Received the Secretary August 29, 1942. 


Practice and Standard Specifications for Concrete and Reinforced Rept. 
the Joint Committee Standard Specifications for Concrete and Reinforced Concrete, Proceedings, 
Am. Soc. E., June, 1940, Pt. paragraph 805(c), 46. 
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valuable—namely, rapid convergence, simple tabulation and ease checking. 
The method somewhat lengthy several balances are required. im- 
portant, therefore, that the estimates the M’-values should reasonably 
close. 

the insertion, their appropriate columns, the M’-values and K-values 
already obtained, one operation usually sufficient give results 
accurate enough form the basis for design. 

point favor Mr. Cornish’s method, for many designers, will its 
clear derivation from the well-known slope-deflection formulas. The substi- 
tution values formulas practice well established among structural engi- 
neers, fact which doubt explains the many efforts express the processes 

example section simple building frame shown Table The 
and K-values are those obtained previous approximation. Check 
calculations made the usual slope-deflection and full moment-distribution 
methods give results practically identical with the shown Table 
The necessary summations are common any method and are not given. 


For example, see Practice and Standard Specifications for Concrete and Reinforced 
Rept. the Joint Committee Standard Specifications for Concrete and Reinforced Concrete, 
Proceedings, Am. Soc. E., June, 1940, Pt. Appendix 110. 
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CLASSIFICATION IRRIGABLE LANDS 


Discussion 


direction when the state and federal agencies concerned with the development 
new agricultural lands learned that, when the engineering problems sup- 
plying water had been solved, there was still the question the suitability 
the land taken into consideration. The sequence these two actions 
vital importance and when properly timed may avoid many disappointments. 

common knowledge that water has been made available some lands 
that could not profitably farmed even though water were free. Although 
not recent discovery the Bureau Reclamation that the inherent qual- 
ity the soil major importance, the Bureau not entirely without fault 
even some its later projects. long ago 1925 the writer was employed 
make soil classification Reclamation Service project aid 
making adjustments water charges which had been found out line 
with the producing power the soil. 

Considerable advancement has been made the land-classification con- 
cept since the passage the 1924, and course that 
should be. Mr. Johnston has presented the results some years 
Bureau Reclamation experience attempting arrive sound basis for 
land classification. 

The author intimates that the purpose land classifications deter- 
mine what lands are arable (using his definition the word and 
that this information will used the design the irrigation system. Un- 
fortunately, the timing these land-classification surveys faulty, and even 
some the recent surveys have not been completed until after construction 
well under way completed. 

For many years the University California, Berkeley, has urged 
detailed soil survey and land classification the East Imperial Mesa Cali- 
fornia, area some 225,000 acres desert, for which water might 


paper Johnston, was published May, 1942, Proceedings. Discussion 
this paper has appeared Proceedings, follows: May, 1942, Earl Storie, Esq. 

Drainage Engr., Div. Soils, Univ. California, Berkeley, Calif. 

Received the Secretary August 31, 1942. 
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made available from the All-American Canal. The classification the soils 
this mesa example poor timing, was not until after the comple- 
tion the All-American Canal and the Coachella Branch (from which branch 
much the land irrigated) that the Bureau Reclamation undertook 
land-classification survey. The writer member the Advisory Com- 
mittee Land Classification for this area and assisted setting the stand- 
ards classification. These followed substantially the system described 
Mr. Johnston’s paper. 

The Coachella Branch was constructed with capacity sufficient irrigate 
nearly 200,000 acres this mesa. This was based very rough recon- 
naissance soil survey made 1920, which would not now qualify recon- 
naissance land classification described Mr. Johnston. The detailed 
classification completed 1942 places the arable lands the mesa sub- 
stantially lower amount—in fact, much lower that, had the information been 
available sooner, the construction canal substantially different cross 
section might have resulted. Although the writer does not have the data 
hand, probable that somewhat similar situation may found the 
Columbia Basin area the State Washington. 

With regard the East Imperial Mesa survey, the Advisory Committee 
was full agreement the standards established, having mind the par- 
ticular use that was made the survey, but the Committee strongly 
advised that detailed soil survey made either before concurrently with 
the land classification. This suggestion has been followed; the University 
California cooperated with the Bureau Reclamation and the two surveys 
were made concurrently, the advantage both agencies. 

One criticism that may made the land classifications outlined 
this paper that the procedure such that can followed engineers 
whose training soils may limited. The Bureau Reclamation personnel 
consists largely men engineering training and has been exceedingly 
slow going far afield for men trained other branches science. Land 
classification not strictly engineering job and, aside from the reconnais- 
sance type survey, should not made engineers. 

The writer frequently has disagreed with the Bureau Reclamation its 
contention that detailed soil survey data are not directly applicable its needs 
and that the Storie Index for soil rating cannot applied generally soils 
outside California because the basic principles used the index are profile 
development and texture the surface soil. 

Closer correlation land classification can maintained among different 
areas the same different states through the use the Storie Index and 
Earl Storie’s Land Types” than any method with which the 
writer familiar; but they must applied men familiar with soil. Land 
classification, made for whatever purpose, only special adaptation soil 
survey information. 

Mr. Johnston’s statement that soil characteristics vary greatly that 
difficult make uniform application guide irrigability, and his 
comparison conditions the Sacramento Valley and eastern Montana, will 
hardly support contention that detailed soil surveys are not the most logical 
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basis for land classification. This contention may have some merit the 
classifiers are not soil trained, but certainly not true little closer co- 
ordination could secured between the engineer and the soil technician. 

has been found the work Mr. Storie and the writer the University 
California that, the more detailed the soil survey, the more accurately 
land classification can made. Their most severe criticism has come from 
attempts make classifications the absence up-to-date soil data. 
up-to-date soil data meant detailed soil survey which not only delineates 
and describes the soil profile terms its development and texture the 
various horizons and strata, but also contains information alkali, drainage, 
water-holding capacity, slope, and erosion. Many the published soil survey 
reports are admittedly not satisfactory for the purpose, and only the most 
recent surveys are adequate. This condition only expected, however, 
soil surveying well land classification has developed rather rapidly, 
from modest beginning, into science. 

The writer can agree with Mr. Johnston statement that future drain- 
age requirements new irrigation projects are difficult and often impossible 
determine advance. has never fully agreed with some persons (not 
engineers) who have advocated the construction drainage works the same 
time the construction irrigation works. the other hand, knowledge 
the soil, its age, profile development, mode formation, parent material, 
and other features that may furnished soil survey are material aid 
predicting what lands may subject seepage, waterlogging, salt accumu- 
lation, and their placement the proper land-classification group. 


land-classification work which Mr. Johnston 
and the writer were closely associated before and after the enactment the 
1924 amendment the Federal Reclamation Act has been developed high 
state precision and usefulness. Members the five-man board review 
class and class lands had opportunity test the accuracy and value 
the Columbia Basin project classification July, 1941. similar procedure 
may used soil survey and classification irrigation-project lands 
undertaken for the Venezuelan Department Public Works beginning ap- 
proximately September 15, 1942. 

The usable water capacity the soil perhaps its most important physical 
property land-classification work for irrigation. The water capacity the 
soil profile used plant roots seems more useful criterion than definite 
limitations soil depth. successful projects with surface irrigation 
‘Northwestern United States, this should least acre-in. Some improve- 
ment soil profile water capacity may effected the use deep-rooted, 
humus-building legumes and farm manures from livestock enterprises. 

irrigation permanent and repay those financially concerned, 
should preceded rigid land classification. This paper agrees closely 
with the writer’s experience. 


4Soil Scientist Chg., Oregon Agri. College and Experiment Station, Oregon State College 
Corvallis, Ore. 
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Discussion 


ALFRED NILEs, Assoc. AM. Soc. 


tion described this paper very ingenious, and should prove great 
timesaver the solution many types problems. Although the author 
has shown applications both beams and columns single span, has failed 
warn the reader that his method not directly applicable continuous 
members, for which the bending moments over the supports must obtained 
the use the three-moment equation, the method moment distribution, 
equivalent method. Perhaps the most awkward member this type, 
from the point view the stress analyst, continuous beam nonuniform 
section that subjected combined bending and compression. This 
handled most readily the method moment distribution, proper allowance 
being made for the effect the axial load when computing fixed-end moments 
and the carry-over and stiffness factors. The computation these quantities 
for members nonuniform section previously published methods slow 
and tedious procedure. The methods proposed the author appear more 
convenient for this work than any other that has yet been suggested, including 
that the writer and Newell.’ 

his numerical examples the author divides the beam into segments 
equal length. This greatly simplifies the work, and nearly always allowable 
when the transverse load (or ‘‘angle change’’) can completely represented 
asmooth curve. this not allowable, however, when unequally spaced 
concentrated loads are present, much the advantage lost. fact, the 
author’s method becomes practically the same that the only 


paper Newmark, Assoc. Am. Soc. E., was published May, 1942, 
Proceedings. Discussion this paper has appeared Proceedings, follows: May, 1942, Bruce 
Johnston, Assoc. Am. Soc. E.; June, 1942, Ketchum, Jr., Assoc. Am. Soc. E.; and 
September, 1942, Messrs. John Wilbur, Ralph Stewart, and Stefan Fraenkel. 

Prof., Aeronautic Eng., Leland Stanford Junior Univ., Stanford Univ., Stanford University, Calif. 

Received the Secretary August 17, 1942. 


Structures,” Niles and Newell, Ed., New York, Y., 1938, Vol. pp. 
133-136, and Vol. II, pp. 


Ibid., Vol. pp. 58-60. 
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differences being the method recording the computations, and the more 
accurate, although more time-consuming, method allowing for the curvature 
the loading (or diagram, which the writer treats composed 
straight segments. most practical problems, however, the value 
the greater accuracy the method allowing for this factor 
questionable. 

computing “equivalent concentrated angle the author uses 
formulas from Figs. and that require division 12, 24, depending 
the shape the angle-change curve. would seem simpler include the 
factor 12, the “common factor” which the moments deflec- 
tions are multiplied get the final results. This would result most 
the values the tabulated computations being 12, times large 
those obtained the method described, but would not affect the final results. 
Care would have taken, however, use the same factor throughout the 
span when part the loading curve was straight and part curved, but any 
given problem would easy multiply the formulas Figs. and 


order obtain common denominator for all the formulas used. 


some his examples (as Fig. the author computes the actual 
value the shear the left end the beam before computing the shears 
other points, whereas others (as Fig. 1(e)) starts the shear computations 
from arbitrarily assumed figure and makes final correction the bending 
moments, necessary. The writer has two objections the latter practice, 
although seems preferred the author. The first that often 
necessary know the shears various points along the span, and, the latter 
practice, would necessary remember correct the values originally 
found obtain the true ones. This could done easily, and the objection 
would unimportant stood alone. The more serious objection that 
the practice eliminates valuable internal check the computations. the 
actual shear the left end first computed, then the moment the right end, 
computed summation the shears along the span, should the same 
that stated the formulation the problem. the two values are not 
substantial agreement, error has been made. the method Fig. 
one does not know whether the necessary moment correction due solely 
the difference between the assumed and actual shears the left end, whether 
partly due numerical error computation. Since the actual shear 
the left end can computed quite easily, the check obtained justifies the little 
additional work involved using it. 

The writer notes that the author has reversed the usual convention and 
considers that loads are positive when they act upward. heartily indorses 
this practice. wonders, however, why the same change was not made 
the conventions for slope and deflection. That would have involved the elim- 


EI > but 
that sign not essential. there only reconcile some independently 
assumed conventions which proved lacking logical consistency. 


really much simpler assume upward loads and deflections positive. Then 


ination the minus sign from the definition ‘‘angle 
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one can differentiate the equation the elastic curve four times, obtaining 
successively the slope, bending moment, shear, and loading, without having 
remember reverse signs arbitrarily various steps. 

related point connection with nomenclature. his paper the 


author terms the quantity, the “angle Actually, states, 
4 
using that expression allowable approximation for what mathe- 
call the which the radius curvature. Since 


the mathematicians already have given the quantity name, why was 
“curvature curve,” although the expression should quite clear ‘‘angle- 
trated angle What the author has termed really 
“rate slope and the latter term would really preferable the 

studying the numerical examples, the writer was unable verify one 
the author’s Fig. 15(a) the equivalent concentrated angle change 
the section change moment inertia shown —404.90. This 
appears quantity obtained use the formulas Fig. 5(a), 
assuming the distributed angle-change curves produced have ordinates either 
one tenth of, ten times, those the actual curve, the adjacent segments 
the beam. this basis, the concentrated angle change question would 
value and the practical consequence, but would interest- 
ing learn whether the figure the text was computed some other method. 

Although the paper subject the foregoing minor criticisms, 


deserves much for developing valuable new tool for the use the 
structural engineer. 


Correction for Transactions: May, 1942, Proceedings, page 715, Fig. 
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STATISTICAL ANALYSIS HYDROLOGY 


Discussion 


STANDISH HALL, AND ALDEN FOSTER 


statistical analysis hydrology, this paper should prove valuable many 
respects. The accuracy the projection short-term record statistical 
methods determine the probability the occurrence event over 
longer period time has long been open question. fact, even the 
basis statistics, the probable error such projection 

The theory the frequency distribution curve, which the duration curve 
based, predicated certain assumptions indicated the author. These 
conditions are not rigidly applicable hydrological data, since tacitly 
assumed: 


(1) That there difference any essential respect between the condi- 
tions and methods from which the observations are drawn; and, the observa- 
tions have been made different epochs, that there has not been any essential 
change during that period; 

(2) That the conditions that regulate the appearance any event must 
the same for every series, and also for every observation the series; and 

(3) That the individual events must completely independent one 
another. 


these conditions, the second and third are marked limitations. 
satisfy the first condition stream-flow measurements, changes should 
have occurred from the transfer the measurements from one location 
another, the accuracy with which gage heights are observed stream 
flow measured estimated. well known that this condition has not been 
fulfilled, since there has been progressive improvement over the period 
available records with regard all these items. Recording gages have 
replaced staff gage readings; improvements have been made the technique 
measuring stream flow and estimating flood discharges surface measure- 

paper Beard, Jun. Am. Soc. E., was published September, 1942, Proceedings. 

Hydr. Engr., East Bay Municipal Utility Dist., Oakland, Calif. 


Received the Secretary August 29, 1942. 
Transactions, Am. Soc. E., Vol. (1927), 58. 
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ments, slope observations, water-surface drop contracted 
general, attempt has been made review previous flood discharge observa- 
tions the light more recent knowledge. The same situation also applies 
rainfall observations and most long records are composite observations 
several locations given locality made under varying conditions accu- 
racy, due exposure wind and height gage above ground. Also, assuming 
that the measurements stream flow were all equal accuracy, there have 
been changes many streams due the works man causing changes 
channel storage reservoirs, levees, and channel improvement, changes 
permeability the surface roads, buildings, cultivation land, soil erosion, 
and deforestation lumbering fire. 

The satisfaction the second condition requires that cyclic changes 
affect the hydrologic record. Although the climatic divisions the world 
appear, general, constant, there are known minor fluctuations having 
duration ranging from few years many centuries. the shorter cycles 
the best known are the sunspot cycle and the 35-yr Briickner cycle. 
the longer meteorological periodicities, the glacial and inter-glacial cycles are 
the most pronounced; but even within historic times there definite evidence 
wet and dry cycles. the western United States evidence climatic changes 
have been derived from the width the annual growth the Big Trees 
California, checked and adjusted the variations the levels salt lakes 
the desert area and archeological remains. result the correlation 
these data, the cycles precipitation have been traced back 1,000 
indicating the alternate wet and dry periods, each having dura- 
tions several hundreds Deyiations from the mean extending 
over cycles such length could not ascribed fortuitous chance. 

The third condition that individual events completely independent 
each other cannot realized the hydrologic cycle. Since climate con- 
tinuous process, each event influenced greater less degree the 
weather the preceding days, weeks, months, and some instances even 
years. particular, flood discharges are circumscribed this respect. 
given precipitation may result either moderate runoff flood un- 
precedented magnitude, depending the temperature and precipitation 
the preceding weeks and months, these may affect the soil infiltration rate, 
the accumulation snow the ground sufficiently shallow depth 
permit melting during warm rain. This complex set circumstances 
surrounding the occurrence major floods results marked skewness 
the data, which turn affect unfavorably the accuracy prediction from 
short records future flood frequency. some flood records years 
length, there has been included maximum greatly excess the frequency 
that its occurrence could not possibly predicted mathematical 
analysis. More accurate flood-prediction frequencies would possible 
more study were made precipitation intensities conjunction with correla- 
tion with the other factors affecting runoff. general, rainfall records not 


Measurement and Computation Flood Discharge,” Carl Paulsen, Transactions, Am. 
Geophysical Union, 1939, Pt. II, 


Encyclopedia Britannica, 14th Ed., Vol. 827. 


| 
ik 
_» 


1460 HALL STATISTICAL ANALYSIS Discussions 


exhibit the degree skewness evidenced flood-flow records, and hence can 
projected with relatively greater accuracy. For example, tropical 
semi-tropical climate lacking the influence snow accumulation, the runoff 
from identical storms can vary considerably, depending whether the storm 
occurs soil either dry, average moisture content, saturated. The last 
condition will naturally produce the maximum runoff, and assuming that this 
frequency one fourth, precipitation occurring once twenty-five years 
would produce maximum flood only once one hundred years. Such full 
use available data, although involving much statistical analysis, would 
undoubtedly yield results greater dependability. regions having accu- 
mulations snow the ground, another variable added involving such 
factors saturation the soil before the snow fell, water content the snow, 
depth accumulation, temperature, and precipitation. favorable coinci- 
dente these factors produces flood unprecedented magnitude. 

The author’s determination the location the plotting position each 
series occurrences based the theoretical expression the binomial 
expansion the form the frequency distribution." The lack complete 
mathematical treatment makes this portion the paper obscure, and would 
difficult for one compute the plotting position any number 
rences without the author’s formula. This unfortunate, the author has 
not included complete table plotting positions. 

The basis the formula the proposition that the frequencies 
successes events are given the successive terms the binomial 


which represents the frequency successes and the frequency failures. 
The plotting position determined equating the summation successive 
terms the binomial series 
There some question the writer’s mind the propriety this 
assumption. Considering, for example, 10-yr period, the first term does not 
indicate the probability equal chances failure success secure the 
median the events occurring one tenth the time, but rather the 
cesses the series. 
This illustrated the following tabulation, which shows the 
No. Probability No. Probability 
successes per 100 series successes per 100 series 


Introduction the Theory Udny Yule, Charles Griffin Co., Ltd., 
pp. 291-295. 
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Two successes would normal and would have the greatest probability, 
but cases one success are almost frequent and even nine successes may 
about once 10,000 trials. The probability obtaining more 
100 series events each 60.82 against 39.18 chances 
successes only one success. 

all cases which does not equal increase naturally raises the 
number successes most probable occurrence, but the same time 
increases the dispersion and decreases the skewness. This revealed the 
distribution for compared with that for 100 the expansion 
the series (0.9 namely: 


No. Probability No. Probability 

successes per 100 series successes per 100 series 
0.16 14. 5.13 


this case, ten successes would normal and would have the greatest 
probability, but the frequency from successes nearly equal. The 
skewness the series much less pronounced than the previous example. 
The probability obtaining more successes 100 series 100 events 
each 54.87 against 45.13 chances obtaining from zero nine successes. 
will noted that practically chance failure exists 100-yr record. 
necessary. This indicates the length record necessary for similar prob- 
ability for 0.01. 

The chances successes the series p)" where 1/p increases 
becomes smaller follows: 


Probability 


successes per 100 trials 
0.10 0.90 34.87 
0.05 0.95 35.85 
0.02 0.98 36.42 
0.01 0.99 100 36.60 


other words, the chances failure obtain during the once-in-a-hundred 
event series 100 trials are greater than the chances failure obtain 
during the once-in-ten event series trials. This offset the 
probability one more successes each series trials. Hence, any 
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given series, justified assume that each occurrence representative 
the period and can considered the median for 1/nth the period. 

selecting the data form the duration curve for maximum daily pre- 
cipitation maximum flood discharges, the maximum for each the years 
chosen maximum events can selected regardless the year 
which they occurred. Assuming that the occurrences “universe” 
were known, then, with records years’ length according the 
binomial expansion, the chance failure would occur about the 
years and the successes would grouped the remaining 65% the 
The selection the maximum events would the proper manner arrang- 
ing the data. records this type usual find wet years with two, 
even three, occurrences that exceed the maximums the dry fact, 
many western streams not unusual find wet years having 50% more 
the daily discharges excess the maximum daily flow the driest years. 
Under these circumstances the selection the maximum occurrence each 
calendar year obviously does not afford the best use the available data. 

The author indicates that the Gauss-Laplace normal distribution curve 
the only one having rigid mathematical derivation. Many rigid mathe- 
matical treatments are possible based assumptions differing 
forming the basis the normal curve error. possible that many 
distributions covering all types may follow this law, but, general, hydrologic 
data differ too greatly from make possible any forced agreement. 

The author states (see heading Normal Duration Curve’’) that. the 
transformation the variate often has been used with exceptional 
This assumes that the geometric mean the true average the series, 
and that the logarithms follow the normal law error. The geometric mean 
series always less than the arithmetic mean, and most dats 
have distribution with the mode less than the mean. This does not indicate 
that the geometric mean the most probable value. many series having 
large number zeros the geometric mean obviously does not apply, for 
example the distribution daily precipitation. very cursory examination 
such data will indicate zero the most probable value. Such data cannot 
possibly follow the normal law error. 

The method used plotting the duration curves Fig. not stated, 
the author’s method has been used the location the plotting positions 
would facilitated appending table. The selection the data (which 
24-hr rainfall, inches) perhaps unfortunate, due the fact that 
many wet years there might two more occurrences 24-hr precipitation 

excess the maximum occurrences dry years. The method used 
the statistical array not stated. discontinuous type 
rence, the ones illustrated, best treated other For 
discharges, for example, combination probability methods determine 
rainfall, combined with methods for flood runoff unit 
graph most advisable. However, this criticism does not affect the validity 
the author’s method analysis. 
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Estimating Floods” from the United States,” Water Supply Paper 
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October, 1942 FOSTER STATISTICAL ANALYSIS 1463 


The discussion the accuracy the duration curve very valuable, but 
the author does not explain his method determining the probable error and 
the largest deviation sufficient extent indicate the applicability his 
results other statistics. the paper argumentative rather 
than analytical, but this not without value, since presents the need for 
complete analysis statistical methods applied hydrology order 
evaluate the worth various probability methods brought forth over the 
past years and determine the proper application the various types 
probability studies engineering statistics. The investigation this field 
hydrology would merit further study. 


connected with the use duration curves the analysis hydrological rec- 
ords contained this paper. number these questions have been treated 
some repetition, the following comments are submitted: 

Mr. Beard devotes considerable space demonstration his conclusion 
that the generally accepted method for plotting items the duration curve 
error, although concedes that, general, the presumed error 
negligible. writer has been unable follow the demonstration Mr. 
theorem. 

the first place, states (see heading the Duration 
that largest occurrence 10-yr period considered the 
median the largest occurrences all the 10-yr periods ever 
the opinion, this statement not correct. the duration curve 
generally considered, the largest occurrence 10-yr period should 
taken the median the largest occurrences 100-yr record, the 
median the 100 largest occurrences 1,000-yr record. 

make this clear, assume duration series, each covering years, 
taken from 100-yr record: 


Also assume that these series have been arranged that the terms are 
descending order magnitude from Mr. Beard’s method, 
only one 10-yr series available, the term that series assumed equal 

Now quite possible that one more the ten terms may greater 
than the smallest term, Consequently, duration series prepared 


With Parsons, Klapp, Brinckerhoff Douglas, New York, 
Received the Secretary August 29, 1942. 
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using the entire 100-yr record, the first terms this combined series may 
include some the terms, which case will not include all the terms. 
Hence, the median the “largest all the 10-yr (a, 
will smaller than the median the largest occurrences the 
100-yr period. 

true (see heading ‘‘Construction the Duration Curve’’) that 
probability that the maximum occurrence set chosen random will ex- 
ceed the median the maximum occurrences all the sets The 
probability that the said maximum occurrence will less than the median 
also 50%. However, the probability that the ten occurrences chosen 
random will exceed the median the maximum occurrences all the 
not 50% stated the author. The 50% probability previously noted 
applies the maximum occurrences (the aforementioned terms); that 
the probability that any one the ten terms one series (such 
for example) will exceed the median the ten terms 0.50 0.10 0.05, 
5%. Hence the probability that none the ten will exceed 
the median should 95% instead 

seems obvious the writer, therefore, that Mr. method 
puting the percentage occurrence corresponding the median the maxi- 
mum occurrences from all the sets cannot justified. The same conclusion 
applies his method for computing the other plotting positions. the 
writer’s opinion, Eq. gives the theoretically correct plotting position. 
cidentally, this method has been quite generally accepted for the plotting 
dufation curves hydrological studies during the years since the first 
World War. 

The writer full agreement with most Mr. Beard’s discussion dura- 
tion curves, except the general use the normal, Gaussian, probability 
curve. has been quite well established, both actuarial studies and 
various statistical studies other sciences, such biology, that unsymmetrical 
distributions many types occur. There seems reason why similar 
distributions should not occur hydrology and meteorology. fact, 
the writer’s observation that the unsymmetrical, skew, distribution the 
rule rather than the exception, hydrological data. The mathematical form 
the skew probability duration curve has been discussed extensively else- 
and noted Mr. Beard, and would seem outside the 
scope the present discussion. 

Mr. Beard’s study the accuracy the duration curve interest. 
The writer fully concurs his statement (see heading the Dura- 
tion Curve’’) that the advantages design the duration curve 
rather than the maximum occurrence record that greater accuracy 
The discussion the use such terms probable 
also well taken. 

conclusion, Mr. Beard should commended for his description the 
duration curve and its application problems hydrology. Although the 
writer does not agree with Mr. Beard’s method locating plotting points for 
the curve, this discrepancy only minor significance the practical 
tion the curve. 
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ALLOCATION THE TENNESSEE VALLEY 
AUTHORITY PROJECTS 


Discussion 


problem allocating the cost enterprise the various services renders 
the commodities produces difficult, even when the services rendered 
the commodities produced are all the same general kind. The problem 
allocation still more complicated for multiple-purpose enterprises projects 
which the services products are very different kinds. 

Strictly speaking, the so-called theories joint-cost allocations, with the 
possible exception the use-of-facilities theory, are never actually what 
their names indicate; rather, they are mathematical methods that arbitrarily 
assign responsibility for the joint cost the various purposes project. 
These methods are just applicable single-purpose project 
with respect the various users its services multiple-purpose project 
with respect its various purposes. not difficult imagine situations 
where, for single-purpose project, the application these cost-allocation 
theories would lead inequitable and discriminative results. There are many 
cases which the cost obtaining given service one way has conceivable 
relation the cost obtaining another way; for example, the customers 
water utility may able supply their own requirements privately 
varying unit costs, depending their respective There 
basis for assuming that the cost the utility supply these various customers 
would bear any relation the costs such users supply themselves. Since 
the application the alternative cost any other method may lead 
arbitrary and discriminative results when applied such relatively simple 
case, follows that the same result may obtained when these same methods 


paper Theodore Parker, Am. Soc. E., was published December, 1941, 
roceedings. Discussion this paper has appeared Proceedings, follows: February, 1942, Malcolm 


Am. Soc. E.; March, 1942, Chandler, Am. Soc. E.; May, 1942, Eugene 
Am. E.; and June, 1942, Messrs. Sherman Woodward, and Albert Arledge. 
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are applied the more complicated problem allocating cost for multiple- 
purpose projects. stated the author, however, the more generally 
accepted methods for allocating the cost single-purpose project scarcely 
would adaptable multiple-purpose project because the absence 
common characteristics the services the various purposes such projects. 

Even for electric utilities, there one generally accepted method for 
allocating cost the various classes customers. The reason that cost alloca- 
tions are made such cases determine the reasonableness the charges 
made obtain basis for determining what charges should 
gardless the allocation cost and the resulting charges for services, the total 
annual revenue received the utility would depend only the total 
annual cost rendering service, provided, course, that such cost less 
than the cost alternative supply. 

The justification multiple-purpose project, from engineering point 
view, depends only the total annual cost the services from the project 
being less than the minimum total annual justifiable cost (that is, cost which 
less than the annual benefits) obtaining the same services otherwise. The 
total annual revenue, for single-purpose project, should depend only the 
total annual cost rendering the services, provided not greater than the 
minimum annual cost obtaining the same services otherwise. Except for this 
limitation, the total annual revenue should independent both the cost 
obtaining the same services otherwise and any cost allocation that may 
made. 

The only reason for making cost allocation obtain basis for the 
rates and charges for various services that will result total annual revenue 
nearly equal the total annual cost possible. the public policy 
not collect revenue for various services government-owned, multiple- 
purpose projects and make only nominal charges for others, the resulting 
deficiencies will have subsidized from general taxes. The amount 
such subsidies should determined directly matter public policy and 
not indirectly (as most often done), means cost allocation. 

The questions whether not project should subsidized and the 
extent the subsidy are political questions rather than problems engineering 
economics. The application engineering economics project whole 
ends when has been determined whether not the project economically 
justified; and the political decisions with respect such projects, rather than 
the engineering and economic facts, should based the declared public 
policy. The answer the question what extent, any, the public 
policy subsidize such projects depends turn the public’s knowing the 
true annual cost, including the controversial question the money.” 

The money invested belongs the public regardless whether the project 
owned; the one case the money invested belongs 
specific individuals groups and the other the public generally. The 
economics project not change merely because sponsored the 
government. important difference between the financing project 
the government private company that the government, particularly 
the federal government, has much greater resources. apparent difference 


— 2p 


‘ 
| 


October, 1942 WATCHORN AND ALLNER TVA PROJECTS 1467 


arises from the fact that, because tax regulations and the safety investors 
government securities provided its taxing powers, the government 
usually able borrow lower rate than private enterprise. This does not 
necessarily mean, often assumed, that the cost money the government 
for any particular purpose less than that responsible private company for 
the same kind undertaking. This true even though the government, 
desires, may finance project entirely with borrowed money, whereas 
private companies must usually provide substantial portions equity money 
addition that borrowed. 

frequently forgotten that the government’s “equity” provided the 
taxpayers who must bear the burden all losses suffered account invest- 
ments such projects. difficult understand logic that assumes the 
public places less value its money and willing for the money earn lower 
return simply because was invested through the government instead 
private company. The answer that the money used the same project 
worth the same, regardless who made the investment, assuming equally 
capable management each case, since the risks would then the same. 
Further, not believed that the interest rate paid the government 
alone tells the entire story the cost the public’s money used the 
government for such purposes. The present situation taxes makes this 
point more clear. The unevaluated portion the cost the public’s money 
used the government the necessary sacrifice the public make money 
available the government for its various activities. This fact has always 
existed and only made more clear present conditions. For example, 
one assumes public debt (as result pay-as-you-go program), the 
government would pay interest the money expended for any purpose. 
Should one then assume zero cost for the money spent the government for 
such projects? The fact the matter that there magic government 
ownership, and the cost money the average return demanded the 
public itself for investments reasonably safe private undertakings. Such 
basis for determination serves indirect, but nevertheless reliable, method 
evaluating previously neglected elements cost money the government. 
The economic justification the project and its true annual total cost should 
both computed the foregoing basis for determining the cost money. 
The total annual revenue should equal this total annual cost make the 
project self-supporting. However, may the public policy that revenue 
received from some the purposes the project, with only nominal revenue 
from others, and normal revenue from the remaining purposes. 

illustrate, assume multiple-purpose project built for flood control, 
navigation, power and energy production, and recreation; and assume that 
the public policy that direct charges made for flood-protection benefits, 
that only nominal charges made for use the recreation facilities, and that 
normal charges made for navigation and power and energy. The total 
annual cost the project $3,700,000, which $300,000 the annual cost 
incurred specifically for flood control, $400,000 specifically for navigation, 
$1,500,000 specifically for power and energy production, $100,000 specifically 
for recreation, and the remainder $1,400,000 for joint facilities. nominal 
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charges made for the use the recreation facilities result annual revenue 
$20,000. The problem determine the rates and charges that should 
made for the use the navigation facilities and for power and energy. 
believed, stated previously, that there correct engineering solution 
this problem, the scientific sense, and that the determination the cost 
allocation the various purposes necessarily arbitrary, regardless how 
may derived. 

The total annual revenue desired likewise matter public policy, which 
one may assume would either (1) large possible, (2) approximately 
equal the total annual cost the project even though flood control 
wholly, and recreation largely, subsidized, (3) equal some specified amount, 
(4) less than the total annual cost the project some specified amount, 
(5) equal some specified percentage the total annual cost the project. 

believed that the underlying principle the vendibility theory cost 
allocation (namely, that services and commodities jointly produced are mutual 
by-products each other and that the actual total annual revenue should 
nearly possible equal the desired total annual revenue) the most logical 
basis for arriving equitable, but admittedly arbitrary, basis for allocating 
the cost government-owned, multiple-purpose project its various pur- 
poses. The total revenue received account any the purposes the 
project immaterial, long not greater than that required defray 
the annual cost for the most economical means obtaining the same services 
otherwise. 

The charges for each the services should equal less than the 
smaller (1) the economic rate (that is, the rate that would produce the 
maximum net income), and (2) the rate required paid for the most eco- 
nomical available means obtaining the same services otherwise. 

the total annual revenue resulting from the predetermined nominal rates 
for certain services and those determined the foregoing manner for the other 
services (which may called either separately total the unadjusted 
revenue) greater than the total desired for the project (as determined 
policy) the rates should further reduced that the desired over-all 
results will obtained. Such excess revenue may result from various 
combinations different amounts earnings account the various pur- 
poses the project. Whatever the combination that causes this excess, 
suggested that the reduction based comparisons the unadjusted 
annual revenue with the annual cost the specific facilities for each the 
purposes the project, and the total unadjusted annual revenue with the 

total revenue desired. 

The simplest case results when the excess, any, unadjusted total annual 
revenue over the desired total annual revenue derived from only one 
the project. This will the case when there only one purpose for which 
the unadjusted annual revenue exceeds the annual cost specific facilities and 
this excess greater than the amount the unadjusted total annual revenue 
exceeds the desired total annual revenue. would appear illogical, under 

such circumstances, further reduce the annual revenues for those purposes 
for which the unadjusted annual revenues fail cover the annual cost 
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specific facilities; and consequently, such reduction should all applied the 
one purpose for which the unadjusted annual revenue exceeds the annual cost 
specific facilities. However, the event (even though there only one 
purpose for which the unadjusted annual revenue exceeds the annual cost 
specific facilities for that purpose) that the excess the unadjusted total 
annual revenue over the desired total annual revenue greater than such excess 
for that one purpose, the situation becomes rather complex and more com- 
plicated treatment the problem, discussed later, necessary. 

there are several purposes the project for which the unadjusted annual 
revenues exceed the annual cost specific facilities, there are, brought out 
the preceding special case, two possible conditions, one which relatively 
simple and the other more complex. The simple condition results when the 
excess the unadjusted total annual revenue over the total annual revenue 
desired less than the total the excesses the unadjusted annual revenues 
over the annual cost specific facilities for those purposes the project for 
which there are such excesses. reduction made for the simple case 
should proportional basis those purposes the project having such 
excess revenues. 

illustrate, assume for the example given previously that the unadjusted 
annual revenue $1,500,000 for navigation and $2,800,000 for power and 
energy, which for navigation $1,100,000 and for power and energy $1,300,000 
excess the annual cost the specific facilities for those purposes. These 
unadjusted annual revenues, together with those from recreation, give 
$4,320,000 the total annual revenue from the project. the desired annual 
revenue $3,000,000, the total reduction for navigation and power and energy 
would $1,320,000, whereas the sum the excess revenues for these purposes 


$2,400,000. The reduction applied navigation then would 


$715,000. The resulting adjusted revenue for navigation would $895,000, 
and for power and energy $2,085,000. 

The more complex situations result when the desired total annual revenue 
the project low, comparison with the unadjusted total annual revenue, 
that the necessary reduction greater than the total excess unadjusted 
revenue over the annual cost specific facilities for those purposes the project 
for which such excess revenues exist. (The case which there was purpose 
the project for which the unadjusted annual revenue exceeded the annual 
cost specific facilities would also included here.) order make the 
desired reduction total revenue for such cases, not only revenues from other- 
wise purposes the project, any, must reduced below the 
annual costs specific facilities for these purposes, but revenues from the 
“unprofitable” purposes the project may also further reduced. 

The solution suggested herein involves the finding ratio, less than unity, 
which when applied the specific cost facilities for any purpose will give 
the adjusted revenue for that purpose, provided the revenue reduced thereby 
(that is, the revenue for any purpose reduced must remain un- 
changed). The ratio determined must also satisfy the condition that ad- 
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justed revenues (including those that remain unchanged) must total the 
desired revenue from the project. 


illustrate, assume the same data for the project previously 


except that the total desired annual revenue only $1,000,000 (that is, the 
subsidy has been increased $2,000,000 compared with the previous example). 
The total annual cost specific facilities, excluding those for flood control, 
$2,000,000. This total multiplied 50%, the ratio the desired annual 
revenue the total annual cost specific facilities, would give the desired 
$1,000,000 total revenue. immediately apparent, however, that the 50% 
ratio applied the cost facilities for recreation ($100,000) would give 
revenue excess the $20,000 predetermined for this purpose. Consequently, 
the computation the ratio must adjusted excluding the annual cost 
recreation facilities and recreation revenue. The total annual cost specific 
facilities for navigation and power and energy $1,900,000, and the desired 
revenue, after deducting $20,000 for recreation, $980,000. The ratio de- 
sired revenue cost then 0.516. The adjusted revenue for navigation would 
then $400,000 0.516, $206,400, and for power and energy $1,500,000 
0.516, $773,600. The total adjusted revenue, including $20,000 from 
recreation, the desired $1,000,000. 

These two methods distributing the excess the unadjusted total annual 
revenue over the total annual revenue desired among the revenues from the 
various purposes multiple-purpose project have common point which 
their results necessarily agree. this point the desired reduction total 
revenue just equal the sum the excesses unadjusted annual revenues 
over the annual cost specific facilities for those purposes for which excesses 
exist. 

This suggested basis for determining the charges made for the various 
services rendered multiple-purpose projects not, strict sense, basis 
for allocating cost. merely basis for the determination the charges 
made for these various services accordance with declared public policy, 
not only what purposes the project are subsidized, but also 
what extent the project whole subsidized. The procedure used 
the author (assuming that the rates for the various services are based the 
cost allocation) results subsidy for the project whole which varies 
the cost allocated the various purposes varies. The subsidy increases 
the cost allocated the no-revenue and small-revenue producing purposes 
the project increases. the other hand, the application the foregoing 
suggested basis for determining the charges made for the various services 
the project assures the return that desired from the whole 
accordance with public policy, and the public policy what purposes shall 
subsidized no-revenue and small-revenue producing. same time 
makes available the various services, for which normal charges are made, 
the same less cost than the minimum for which they could obtained 
otherwise. Thus, the beneficiaries such services have basis for complaint, 
since they are not hurt, whereas the interests the general public are pro- 
tected definite predetermined extent. 
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